Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 834,164 A→G 100% intergenic (‑830/‑236) mcpA1 ← / → hoxH methyl‑accepting chemotaxis protein McpA/NAD‑reducing hydrogenase HoxS subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034834,1640A→G100.0% 820.0 214intergenic (‑830/‑236)mcpA1/hoxHmethyl‑accepting chemotaxis protein McpA/NAD‑reducing hydrogenase HoxS subunit beta
Reads supporting (aligned to +/- strand):  new base (114/100):  ref base (0/0):  total (114/100)

Only 100 of 217 total aligned reads displayed.

TCGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGAAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACCAC  >  NC_014034/834066‑834261
                                                                                                  |                                                                                                 
tCGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAg                                                                                                 >  1:3260157/1‑101 (MQ=255)
 cGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAgg                                                                                                >  1:4431154/1‑101 (MQ=255)
  ggACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAggg                                                                                               >  1:373213/1‑101 (MQ=255)
   gACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGc                                                                                              >  1:2672157/1‑101 (MQ=255)
   gACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGc                                                                                              >  1:1785697/1‑101 (MQ=255)
   gACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGc                                                                                              >  2:1493233/1‑101 (MQ=255)
     cGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCg                                                                                            >  1:1815349/1‑101 (MQ=255)
     cGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCg                                                                                            <  2:2698536/101‑1 (MQ=255)
      gCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGc                                                                                           >  2:4166929/1‑101 (MQ=255)
      gCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGc                                                                                           <  2:3649483/101‑1 (MQ=255)
       cAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCa                                                                                          >  2:3570043/1‑101 (MQ=255)
        aGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAg                                                                                         >  2:1136847/1‑101 (MQ=255)
        aGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAg                                                                                         <  1:4241441/101‑1 (MQ=255)
        aGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAg                                                                                         <  1:3760759/101‑1 (MQ=255)
          cGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGAt                                                                                       <  2:350776/101‑1 (MQ=255)
            ggCAATTCGTTTGACCCTCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCt                                                                                     <  2:4602827/101‑1 (MQ=255)
            ggCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCt                                                                                     <  2:2007427/101‑1 (MQ=255)
            ggCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCt                                                                                     >  1:44113/1‑101 (MQ=255)
              cAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTg                                                                                   <  2:2158937/101‑1 (MQ=255)
               aaTTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTgg                                                                                  >  1:1759485/1‑101 (MQ=255)
                aTTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGa                                                                                 <  2:3986002/101‑1 (MQ=255)
                 ttCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGaa                                                                                <  1:1676545/101‑1 (MQ=255)
                   cGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGa                                                                              <  2:4441270/101‑1 (MQ=255)
                    gTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGAt                                                                             <  1:4530747/101‑1 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            <  2:2818745/101‑1 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            >  1:1525527/1‑101 (MQ=255)
                      ttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATga                                                                           >  2:1358095/1‑101 (MQ=255)
                       tGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgac                                                                          >  1:2953048/1‑101 (MQ=255)
                        gACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgacg                                                                         <  1:1473844/101‑1 (MQ=255)
                          cccgcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgacgac                                                                       >  2:3678590/1‑101 (MQ=255)
                           ccgcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGacgacg                                                                      <  2:2631563/101‑1 (MQ=255)
                            cgcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACgg                                                                     >  2:4296684/1‑101 (MQ=255)
                             gcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACggg                                                                    >  2:581302/1‑101 (MQ=255)
                             gcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACggg                                                                    >  2:715336/1‑101 (MQ=255)
                              cccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACgggg                                                                   >  1:2093567/1‑101 (MQ=255)
                                cgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCACGGGGGGCCGCAGATCTTGGAAGATGACGACGGGGTc                                                                 >  1:3994013/1‑101 (MQ=255)
                                 ggCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTcc                                                                >  1:359525/1‑101 (MQ=255)
                                 ggCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTcc                                                                <  1:2902523/101‑1 (MQ=255)
                                 ggCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTcc                                                                <  1:4266531/101‑1 (MQ=255)
                                  gCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCg                                                               <  2:1452212/101‑1 (MQ=255)
                                   cGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCgg                                                              >  2:231745/1‑101 (MQ=255)
                                      gCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTc                                                           <  2:2952819/101‑1 (MQ=255)
                                      gCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTc                                                           <  1:2943163/101‑1 (MQ=255)
                                      gCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTc                                                           <  2:2452999/101‑1 (MQ=255)
                                        tGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCgg                                                         >  2:3761280/1‑101 (MQ=255)
                                         ggCGCAATTTTGGTCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGACCGCAGATCTTGGAAGATGACGCCGGGGTCCGGGTCGGa                                                        >  2:4145968/1‑101 (MQ=255)
                                         ggCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGa                                                        >  2:3709361/1‑101 (MQ=255)
                                          gcgcAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGac                                                       >  2:3295619/1‑101 (MQ=255)
                                          cttccgatctGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGac                                                       <  2:655122/92‑1 (MQ=255)
                                             cAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGAcaca                                                    <  1:1143295/101‑1 (MQ=255)
                                             cAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGAcaca                                                    >  1:2597524/1‑101 (MQ=255)
                                             cAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGAcaca                                                    <  1:975713/101‑1 (MQ=255)
                                              aaTTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAg                                                   <  1:414641/101‑1 (MQ=255)
                                               aTTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAgg                                                  >  1:3517495/1‑101 (MQ=255)
                                                ttttGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAggg                                                 >  1:1205109/1‑101 (MQ=255)
                                                    gATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGcccc                                             <  2:617620/101‑1 (MQ=255)
                                                    gATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGcccc                                             >  1:1350685/1‑101 (MQ=255)
                                                    gATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGCCCGGGTCGGACACAGGGcccc                                             >  1:83418/1‑101 (MQ=255)
                                                    gATCCATCCCAAGGGAAGGGTGGAGGAGGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGcccc                                             <  1:654000/101‑1 (MQ=255)
                                                     atccatccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGccccc                                            <  2:4016225/101‑1 (MQ=255)
                                                       ccatccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGa                                          >  1:130509/1‑101 (MQ=255)
                                                         atccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAAc                                        >  2:4278194/1‑101 (MQ=255)
                                                            ccAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCac                                     >  2:3189836/1‑100 (MQ=255)
                                                             cAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATc                                    <  2:1759366/101‑1 (MQ=255)
                                                              aagggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATcg                                   >  2:18282/1‑101 (MQ=255)
                                                              aagggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATcg                                   <  1:3393132/101‑1 (MQ=255)
                                                              aagggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATcg                                   >  2:2035612/1‑101 (MQ=255)
                                                                 ggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGt                                >  1:1884860/1‑101 (MQ=255)
                                                                  gaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTc                               <  1:274105/101‑1 (MQ=255)
                                                                     gggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGAc                            <  2:1105093/101‑1 (MQ=255)
                                                                      ggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACa                           <  2:2671972/101‑1 (MQ=255)
                                                                      ggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACa                           <  2:3259912/101‑1 (MQ=255)
                                                                       gTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAg                          <  2:4403338/101‑1 (MQ=255)
                                                                       gTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAg                          >  1:2463445/1‑101 (MQ=255)
                                                                       gTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAg                          >  1:3183903/1‑101 (MQ=255)
                                                                        tGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGc                         <  2:1897183/101‑1 (MQ=255)
                                                                          gaggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCgg                       >  2:1298679/1‑101 (MQ=255)
                                                                           aggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCggg                      >  1:2317125/1‑101 (MQ=255)
                                                                            ggaTGTAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGTGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGga                     <  2:876452/101‑1 (MQ=255)
                                                                            ggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGga                     >  1:581651/1‑101 (MQ=255)
                                                                             gaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgag                    <  2:3395849/101‑1 (MQ=255)
                                                                             gaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgag                    >  2:1176338/1‑101 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgaga                   >  2:1087421/1‑101 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgaga                   >  1:1971300/1‑101 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgaga                   >  2:4426667/1‑101 (MQ=255)
                                                                               tGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgagag                  <  2:209947/101‑1 (MQ=255)
                                                                               tGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgagag                  >  1:4497315/1‑101 (MQ=255)
                                                                                 gagaGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGAt                <  1:3809653/101‑1 (MQ=255)
                                                                                   gagGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATcc              >  1:3617474/1‑101 (MQ=255)
                                                                                      gTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGCTCCcgc           >  1:689763/1‑101 (MQ=255)
                                                                                      gTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCcgc           >  1:1418172/1‑101 (MQ=255)
                                                                                       ttCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCcgcg          <  1:1133744/101‑1 (MQ=255)
                                                                                         cGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGcc        <  1:202967/101‑1 (MQ=255)
                                                                                         cGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGcc        >  2:1782609/1‑101 (MQ=255)
                                                                                           ggCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCgc      >  2:1340196/1‑100 (MQ=255)
                                                                                            gCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGAc     <  2:83404/101‑1 (MQ=255)
                                                                                              cgcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACca   >  1:4663878/1‑101 (MQ=255)
                                                                                              cgcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACca   >  1:77420/1‑101 (MQ=255)
                                                                                              cgcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGACGACca   >  1:724652/1‑101 (MQ=255)
                                                                                               gcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACcac  >  2:625133/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
TCGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGAAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACCAC  >  NC_014034/834066‑834261

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG