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% 891.7 232intergenic (‑830/‑236)mcpA1/hoxHmethyl‑accepting chemotaxis protein McpA/NAD‑reducing hydrogenase HoxS subunit beta
Reads supporting (aligned to +/- strand):  new base (124/108):  ref base (0/0):  total (124/108)

Only 100 of 234 total aligned reads displayed.

TCGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGAAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACCAC  >  NC_014034/834066‑834261
                                                                                                  |                                                                                                 
tCGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAg                                                                                                 >  2:803281/1‑101 (MQ=255)
 cGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAgg                                                                                                >  2:3786064/1‑101 (MQ=255)
 cGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGGGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGGGGAGGATGGAGAGGTTCGGGGGGCGGAgg                                                                                                >  2:4788384/1‑101 (MQ=255)
  ggACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAggg                                                                                               <  1:1034555/101‑1 (MQ=255)
  ggACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAggg                                                                                               >  1:2755286/1‑101 (MQ=255)
  ggACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGCTCGGGCCGCGGAggg                                                                                               >  2:3301073/1‑101 (MQ=255)
    aCGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGcc                                                                                             <  1:2404143/101‑1 (MQ=255)
       tcttccgatcttCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCa                                                                                          <  1:3471160/91‑1 (MQ=255)
       cAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCa                                                                                          <  2:3382767/101‑1 (MQ=255)
           gggCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATc                                                                                      >  1:4808474/1‑101 (MQ=255)
            ggCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCt                                                                                     >  1:1639473/1‑101 (MQ=255)
             gCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCtt                                                                                    >  1:1803598/1‑101 (MQ=255)
             gCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCtt                                                                                    <  2:1996214/101‑1 (MQ=255)
              cAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGGGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTg                                                                                   >  2:314600/1‑101 (MQ=255)
               aaTTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTgg                                                                                  <  2:934079/101‑1 (MQ=255)
               aaTTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTgg                                                                                  >  1:4121787/1‑101 (MQ=255)
                aTTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGa                                                                                 >  1:1112726/1‑101 (MQ=255)
                aTTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGa                                                                                 <  2:1658651/101‑1 (MQ=255)
                 ttCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGaa                                                                                <  1:819637/101‑1 (MQ=255)
                 ttCGTGTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGaa                                                                                <  2:2965122/101‑1 (MQ=255)
                  tCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAg                                                                               <  1:4487531/101‑1 (MQ=255)
                    gTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGAt                                                                             <  2:2257186/101‑1 (MQ=255)
                    gTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGAt                                                                             <  1:3935077/101‑1 (MQ=255)
                     tttGGCCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            <  1:24411/101‑1 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            >  1:1928403/1‑101 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            <  1:2108846/101‑1 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            <  2:3393826/101‑1 (MQ=255)
                     tttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATg                                                                            <  1:304231/101‑1 (MQ=255)
                      ttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATga                                                                           <  1:3924491/101‑1 (MQ=255)
                      ttGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATga                                                                           <  2:36240/101‑1 (MQ=255)
                       tGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgac                                                                          <  2:636597/101‑1 (MQ=255)
                        gACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgacg                                                                         <  2:1594908/101‑1 (MQ=255)
                        gACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgacg                                                                         >  1:4635583/1‑101 (MQ=255)
                          cccgcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATgacgac                                                                       <  1:556565/101‑1 (MQ=255)
                           ccgcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGacgacg                                                                      <  1:776414/101‑1 (MQ=255)
                             gcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACggg                                                                    <  2:2099867/101‑1 (MQ=255)
                             gcccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACggg                                                                    >  2:2076327/1‑101 (MQ=255)
                               ccgGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGt                                                                  <  1:1607529/101‑1 (MQ=255)
                                 ggCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTcc                                                                <  2:1099760/101‑1 (MQ=255)
                                  gCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCg                                                               >  2:3908753/1‑101 (MQ=255)
                                   cGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCgg                                                              <  2:2658860/101‑1 (MQ=255)
                                    gTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCggg                                                             >  1:4524737/1‑101 (MQ=255)
                                      gCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGGTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTc                                                           >  1:4724154/1‑101 (MQ=255)
                                      gCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTc                                                           >  1:4557693/1‑101 (MQ=255)
                                        tGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCgg                                                         <  1:609411/101‑1 (MQ=255)
                                         ggCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGa                                                        >  1:3587799/1‑101 (MQ=255)
                                         ggCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGa                                                        >  1:49770/1‑101 (MQ=255)
                                           cgcAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGaca                                                      <  1:3882662/101‑1 (MQ=255)
                                           cgcAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGaca                                                      <  2:4333892/101‑1 (MQ=255)
                                            gcAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGacac                                                     <  2:4490373/101‑1 (MQ=255)
                                              aaTTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAg                                                   <  1:4283560/101‑1 (MQ=255)
                                              aaTTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAg                                                   <  2:4261385/101‑1 (MQ=255)
                                                  ttGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGcc                                               >  1:1667827/1‑101 (MQ=255)
                                                  ttGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGcc                                               >  1:2477131/1‑101 (MQ=255)
                                                    gATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGcccc                                             <  2:2895041/101‑1 (MQ=255)
                                                     atccatccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGccccc                                            <  1:2189574/101‑1 (MQ=255)
                                                      tccatccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCg                                           <  1:120785/101‑1 (MQ=255)
                                                      tccatccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCg                                           >  1:438861/1‑101 (MQ=255)
                                                        catccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGaa                                         >  1:2342710/1‑101 (MQ=255)
                                                        catccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGGCGGACACAGGGCCCCCGaa                                         >  1:3091959/1‑101 (MQ=255)
                                                         atccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAAc                                        >  2:1949922/1‑101 (MQ=255)
                                                          tccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAAcc                                       <  1:3528434/101‑1 (MQ=255)
                                                          tccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAAcc                                       <  1:2780475/101‑1 (MQ=255)
                                                          tccCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAAcc                                       <  1:2096586/101‑1 (MQ=255)
                                                            ccAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCAt                                     >  1:4144464/1‑101 (MQ=255)
                                                             cAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGGCCCCGAACCATc                                    >  2:4603219/1‑101 (MQ=255)
                                                             cAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATc                                    >  2:4537841/1‑101 (MQ=255)
                                                                 ggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGt                                >  1:3002089/1‑101 (MQ=255)
                                                                 ggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGt                                <  2:3925219/101‑1 (MQ=255)
                                                                 ggaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGt                                <  1:408896/101‑1 (MQ=255)
                                                                  gaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTc                               >  1:3856232/1‑101 (MQ=255)
                                                                  gaagggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTc                               <  2:505670/101‑1 (MQ=255)
                                                                     gggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGAc                            <  1:4022591/101‑1 (MQ=255)
                                                                     gggTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGAc                            <  1:4603063/101‑1 (MQ=255)
                                                                       gTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGCCGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAg                          <  1:4271167/101‑1 (MQ=255)
                                                                       gTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAg                          <  1:89571/101‑1 (MQ=255)
                                                                       gTGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAg                          >  1:381949/1‑101 (MQ=255)
                                                                        tGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACCGc                         >  2:4236842/1‑101 (MQ=255)
                                                                        tGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGc                         >  2:4014133/1‑101 (MQ=255)
                                                                        tGGAGGATGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGc                         <  2:104371/101‑1 (MQ=255)
                                                                         ggaggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCg                        >  2:366343/1‑101 (MQ=255)
                                                                         ggaggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCg                        >  2:1533855/1‑101 (MQ=255)
                                                                         ggaggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCg                        >  1:4834007/1‑101 (MQ=255)
                                                                           aggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCggg                      >  2:4439470/1‑101 (MQ=255)
                                                                           aggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCggg                      <  1:2081923/101‑1 (MQ=255)
                                                                            ggaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGga                     >  2:4445232/1‑101 (MQ=255)
                                                                             gaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgag                    <  2:277606/101‑1 (MQ=255)
                                                                             gaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgag                    >  2:3191556/1‑101 (MQ=255)
                                                                             gaTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgag                    >  2:3504102/1‑101 (MQ=255)
                                                                              aTGGAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATAGCGTCGCCAGcacgacc                   >  2:2164618/1‑95 (MQ=255)
                                                                                ggAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgagaga                 >  2:2341434/1‑101 (MQ=255)
                                                                                ggAGAGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGgagaga                 >  2:1062609/1‑101 (MQ=255)
                                                                                 gagaGGTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGAt                <  1:2153548/101‑1 (MQ=255)
                                                                                     ggTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGGGAGATCCcg            >  2:4968940/1‑101 (MQ=255)
                                                                                     ggTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCcg            >  2:2744304/1‑101 (MQ=255)
                                                                                      gTTCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCcgc           >  1:2813267/1‑101 (MQ=255)
                                                                                       ttCGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCcgcg          <  2:303720/101‑1 (MQ=255)
                                                                                         cGGGCCGCGGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGcc        >  2:1352891/1‑101 (MQ=255)
                                                                                               gcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGCGATCCCGCGCCGACcaa  >  2:1720460/1‑100 (MQ=255)
                                                                                               gcgGAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACcac  <  1:1406602/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TCGGACGCAGCGGGCAATTCGTTTGACCCGCCCGGCGTGCTGGCGCAATTTTGATCCATCCCAAGGGAAGGGTGGAGGATGGAGAGGTTCGGGCCGCGAAGGGCCGCAGATCTTGGAAGATGACGACGGGGTCCGGGTCGGACACAGGGCCCCCGAACCATCGCGTCGACAGCGGGAGAGATCCCGCGCCGACCAC  >  NC_014034/834066‑834261

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG