Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,000,979 +G 3.8% intergenic (+23/+100) bhbA → / ← bfr methylmalonyl‑CoA mutase/bacterioferritin

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,000,9791.→G 3.8% 7.9 ‑4.6 ‑0.0132intergenic (+23/+100)bhbA/bfrmethylmalonyl‑CoA mutase/bacterioferritin
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (49/78):  total (54/78)
Fisher's exact test strand distribution p-value = 1.02e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.42e-01

Only 100 of 220 total aligned reads displayed.

CGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCT  >  NC_014034/1000880‑1001079
                                                                                                    |                                                                                                         
cGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGctgacacggggcagctgacacggggca.g                                                                                                          >  2:1174231/1‑101 (MQ=255)
cGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGctgacacggggcagctgacacggggca.g                                                                                                          <  1:2189745/101‑1 (MQ=255)
  aTCTTCGGGCCGGGCACCAACATTCCGGCGGGGGGGCGCGAGATCCTCGATCTGTTCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.g.gg                                                                                                       >  1:2740891/1‑101 (MQ=255)
   tCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.g.ggg                                                                                                      <  1:3744235/101‑1 (MQ=255)
    cTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.g.gggg                                                                                                     <  1:898465/101‑1 (MQ=255)
    cTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.g.gggg                                                                                                     >  2:3404743/1‑101 (MQ=255)
    cTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGAGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.g.gggg                                                                                                     <  2:4490085/101‑1 (MQ=255)
      tCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGgc                                                                                                   >  2:616598/1‑101 (MQ=255)
       cgggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGgcg                                                                                                  <  2:2741804/101‑1 (MQ=255)
       cgggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGgcg                                                                                                  <  1:26640/101‑1 (MQ=255)
       cgggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGgcg                                                                                                  <  2:1413597/101‑1 (MQ=255)
        gggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGgcgc                                                                                                 >  2:3103352/1‑101 (MQ=255)
         ggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCt                                                                                                <  1:3828475/101‑1 (MQ=255)
          gccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTg                                                                                               <  1:3163786/101‑1 (MQ=255)
          gccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTg                                                                                               >  1:4453944/1‑101 (MQ=255)
            cgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.g.gggggggggggc                                                                                             >  2:3237165/1‑94 (MQ=255)
            cgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.cc                                                                                            >  2:537985/1‑101 (MQ=255)
            cgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.cc                                                                                            <  2:1647844/101‑1 (MQ=255)
             gggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGGGCTG.ccc                                                                                           >  1:4340549/1‑101 (MQ=255)
              ggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.cccc                                                                                          >  2:258089/1‑101 (MQ=255)
               gcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.TGGGGGGCCG.cgccc                                                                                         >  2:1612755/1‑101 (MQ=255)
                cACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.cccccc                                                                                        <  2:4021177/101‑1 (MQ=255)
                   cAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGGGGCG.CCCCCCCcc                                                                                     >  1:833418/1‑101 (MQ=255)
                   cAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGcc                                                                                     >  1:299592/1‑101 (MQ=255)
                    aaCATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCGGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCt                                                                                    <  1:1446259/101‑1 (MQ=255)
                    aaCATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCAGG.GGGGGGGGTG.CCCCCCCcc                                                                                     >  1:1108209/1‑101 (MQ=255)
                    aaCATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCt                                                                                    <  1:617857/101‑1 (MQ=255)
                     aCATCCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCtt                                                                                   <  2:3311093/101‑1 (MQ=255)
                        ttCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTgg                                                                                <  1:2681668/101‑1 (MQ=255)
                          ccGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGTGTNG.CCCCCCGCCTTTGG..cc                                                                            >  2:1995159/1‑101 (MQ=255)
                           cggcggcggCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..ccc                                                                           <  1:3884053/101‑1 (MQ=255)
                           cggcggcggCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGAGCTG.CCCCCCGCCCTTGG..ccc                                                                           >  2:2112276/1‑101 (MQ=255)
                            ggggggggcgCGCGAGATCCGGGATCTGACCCGCGGGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..cccc                                                                          <  1:2376798/95‑1 (MQ=255)
                            ggcggcggcGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTC.CCCCCCGCCCTTGG..gccc                                                                          >  1:4640241/1‑101 (MQ=255)
                            ggcggcggcGCGCGAGATCCTGGATCTGACCCGCGCGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..cccc                                                                          <  1:1977466/101‑1 (MQ=255)
                             gcgggcgcccgcgAGATCCGGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..ccccc                                                                         <  2:1540692/92‑1 (MQ=255)
                              cggcggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..cccccc                                                                        <  2:799642/101‑1 (MQ=255)
                              atctggcgGGCGAGATCCTGGATCTGGTCCGCGCGGCGGGGGGCTCACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..cccccc                                                                        <  2:1883174/97‑1 (MQ=255)
                               ggcggggCGCGAGATCCTGGATCTGATCCGCGCGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCg                                                                       <  2:2025964/101‑1 (MQ=255)
                                 cggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGGGCTG.CCCCCCCCCCCTGG..CCCCCCcgg                                                                     >  1:186657/1‑101 (MQ=255)
                                 cggacctcaagaTTCTGGAGCTGATGCGCGGGTCGGGGGACTGACAAGGGGCAGTGGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCggg                                                                     <  2:1709199/92‑1 (MQ=255)
                                  ggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.GCCCCCGCCTTTGG..CCCCCCGGGa                                                                    >  1:3087133/1‑101 (MQ=255)
                                   gggcgcgAGACCCTGGATCTGATCCGCGCGGGGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGat                                                                   <  1:128191/99‑1 (MQ=255)
                                   gcgcgcGAGATCCTGGATCTGACCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGat                                                                   <  2:359256/101‑1 (MQ=255)
                                    cgcgcgAGATCCTGGATCGGACCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACCCGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGata                                                                  <  1:4072036/101‑1 (MQ=255)
                                     gcgcgAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGatat                                                                 <  2:2834541/101‑1 (MQ=255)
                                      cgcgAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATAtt                                                                <  1:151453/101‑1 (MQ=255)
                                      cgcgAGATCCTGGATCGGCTCCGGGCGTGGCGGGGAGAACGGGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATAtt                                                                <  2:3045814/101‑1 (MQ=255)
                                       gcgACATCCTGGATCTGACCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTC.CCCCCCGCCTTTGG..CCCCCCGGGATAttt                                                               <  1:1746995/101‑1 (MQ=255)
                                        cgAGATCCTGGATCTGATCCGCGCGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATAtttt                                                              <  1:2879064/101‑1 (MQ=255)
                                        cgAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGGGCTG.CCCCCCCCCCTTGGGCCCCCCCGGGATTtt                                                                >  1:3245145/1‑101 (MQ=255)
                                          agaTCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTgg                                                            <  1:2974526/101‑1 (MQ=255)
                                             tCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCa                                                         >  1:2301143/1‑101 (MQ=255)
                                              ccTGGATCTGATCCGCGCGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCaa                                                        <  1:2211894/101‑1 (MQ=255)
                                              ccTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCaa                                                        <  2:279797/101‑1 (MQ=255)
                                               cTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GAGGGGGGGGGGCCCCCCCCCCCGTGC..CCCCCCGGGCTCTTTTgcgtg                                                         >  2:1742776/1‑97 (MQ=255)
                                                tGGATCTGATCCGCGGGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGa                                                      <  1:2041838/101‑1 (MQ=255)
                                                 ggATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCCGGATATTTTGGGCAAGAt                                                     >  2:3360271/1‑101 (MQ=255)
                                                  gATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGGGCGG.CCCCCCCCCCTTGG..GCCCCCCGGGTATTTTTGGGAAGaat                                                    >  1:2021450/1‑99 (MQ=255)
                                                  gATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGAtt                                                    <  2:2343530/101‑1 (MQ=255)
                                                  gATCGGCTCCGCGCGGGGCGGGGCTTACGCGGGGCACCTGCCCCGGGGCA.G.GGGGGCCCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGAtt                                                    <  2:2733389/101‑1 (MQ=255)
                                                   tccTGACCCGGGCGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTa                                                   <  1:3171926/99‑1 (MQ=255)
                                                   aTCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTa                                                   >  1:680673/1‑101 (MQ=255)
                                                    tctgtcccGCGCGGCGCGTGCGTGACGCGGGGTATCTTACCCGGGGCA.G.GGGGGCCCTC.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAc                                                  <  1:4111170/95‑1 (MQ=255)
                                                    tCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGCTATTTTGGGCAAGATTAc                                                  >  2:4355454/1‑101 (MQ=255)
                                                     cTGACCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.t                                                <  1:3515867/101‑1 (MQ=255)
                                                           cgcgcgGCGGGGGGCTGACACGGGGCAGCTGACGCGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGaaaaa                                          <  1:1139018/101‑1 (MQ=255)
                                                                  gggggggCGGACGCGGGGCAGCTGACGCGGGGCG.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTc                                   <  2:3911871/98‑1 (MQ=38)
                                                                  ggcggGGCTGACCCGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTc                                   <  2:1283550/100‑1 (MQ=255)
                                                                   gcggGGCTGACACGGGGCAGCTGACACGGGGCAGG.GGGGGGGCGC.CCCCCCCCCCTTGC..CCCCCCGGGGATTTTTGGGAAAGATAAT.TAAAAAAgggcgcc                                   >  2:824459/1‑96 (MQ=255)
                                                                   gcggGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTcc                                  <  1:2657880/101‑1 (MQ=255)
                                                                    gggggCGGACAGGGGGCAGCTAACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCcg                                 <  1:1612860/100‑1 (MQ=255)
                                                                    cggGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCCG.CCCCCCCCCCTTGG..CCCCCCCGGGTATTTTTGGGCAGAATAC.TTAAAAAAGCGCtacc                                 >  2:4277369/1‑98 (MQ=255)
                                                                      gggCGNACGCGGGGCAGCGGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCcgcg                               <  1:1220200/101‑1 (MQ=255)
                                                                         ctgacacggggcagctgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAgg                            >  2:1502858/1‑101 (MQ=255)
                                                                         ctgacacggggcagctgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAgg                            <  1:2577456/101‑1 (MQ=255)
                                                                          tgacagggggcagcggacgcggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAggg                           <  2:2339680/101‑1 (MQ=255)
                                                                          tgacacggggcagctgacagggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAggg                           <  2:4376893/101‑1 (MQ=255)
                                                                           gacacggggcagctgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAgggg                          <  1:331242/101‑1 (MQ=255)
                                                                           gacacggggcagctgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAgggg                          <  2:910388/101‑1 (MQ=255)
                                                                            acgcggggcagctgacacggggca.g.GGGGGCGCTC.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGgc                         <  1:4277657/98‑1 (MQ=255)
                                                                            acacggggcagctgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGgc                         >  1:3054598/1‑101 (MQ=255)
                                                                                ggggcagcggccacggggcg.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCt                     <  2:973602/101‑1 (MQ=255)
                                                                                 gggcagctgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCtt                    <  2:386098/101‑1 (MQ=255)
                                                                                   gcagctgacacggggga.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTg                  <  2:56310/101‑1 (MQ=255)
                                                                                   gcagctgacacggggca.g.GGGGGGGGTG.CCCCCCCCCCTTGG..CCCCCCCGGGGATTTTGGGGAAGAATACTTAAAAAAGGGGCCCCCCGGGGGGGGCCttg                   >  1:1626948/1‑100 (MQ=255)
                                                                                    cagctgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGa                 >  2:2919445/1‑101 (MQ=255)
                                                                                       ctgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCt              >  2:2282276/1‑101 (MQ=255)
                                                                                        tgacacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCtt             >  2:4181715/1‑101 (MQ=255)
                                                                                        ggacacggggcg.g.GGGGGCCCTC.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCtt             <  2:2324097/100‑1 (MQ=255)
                                                                                         gacggggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCttt            <  1:4467682/96‑1 (MQ=255)
                                                                                           cacggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgt          >  2:1493634/1‑101 (MQ=255)
                                                                                            ccggggcg.g.GGGGGCGCTC.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgtg         <  2:1426043/100‑1 (MQ=255)
                                                                                             cggggca.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGc        <  1:1004314/101‑1 (MQ=255)
                                                                                                   a.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCATGATTAC.TGAAAAAGGCGCTCCGCTAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  2:788158/1‑101 (MQ=255)
                                                                                                   a.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  2:109755/1‑101 (MQ=255)
                                                                                                   a.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  1:2769593/1‑101 (MQ=255)
                                                                                                   a.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  1:1935437/1‑101 (MQ=255)
                                                                                                   a.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  2:4661621/1‑101 (MQ=255)
                                                                                                   a.g.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  2:3755716/1‑101 (MQ=255)
                                                                                                    |                                                                                                         
CGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.G.GGGGGCGCTG.CCCCCCGCCTTTGG..CCCCCCGGGATATTTTGGGCAAGATTAC.TGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCT  >  NC_014034/1000880‑1001079

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG