Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,000,979 +G 2.6% 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 2.6% 4.6 ‑4.4 ‑0.0152intergenic (+23/+100)bhbA/bfrmethylmalonyl‑CoA mutase/bacterioferritin
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (47/101):  total (51/101)
Fisher's exact test strand distribution p-value = 1.17e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

GCGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCT  >  NC_014034/1000879‑1001079
                                                                                                      |                                                                                                    
 cGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGct.gacacggggcagctgacacggggca.g                                                                                                     <  2:859103/101‑1 (MQ=255)
 cGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGct.gacacggggcagctgacacggggca.g                                                                                                     <  1:3408958/101‑1 (MQ=255)
  gATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.gg                                                                                                    >  2:4122602/1‑101 (MQ=255)
     cTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.ggggg                                                                                                 <  1:2856376/101‑1 (MQ=255)
      ttCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.gggggg                                                                                                <  1:877141/101‑1 (MQ=255)
      ttCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.gggggg                                                                                                >  2:2713741/1‑101 (MQ=255)
      ttCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.gggggg                                                                                                >  2:3522068/1‑101 (MQ=255)
       tCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGgc                                                                                               <  2:728749/101‑1 (MQ=255)
        cgggccgggcACCAACATTCCGGCGGCGGCGCGCGGGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGgcg                                                                                              >  1:2688876/1‑101 (MQ=255)
        cgggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.gggggggg                                                                                              >  1:436490/1‑99 (MQ=255)
         gggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGgggc                                                                                             >  2:521410/1‑101 (MQ=255)
         gggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGgcgc                                                                                             <  2:1663763/101‑1 (MQ=255)
          ggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCTCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCt                                                                                            <  2:2953786/101‑1 (MQ=255)
           gccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTg                                                                                           <  2:3053734/101‑1 (MQ=255)
            ccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGc                                                                                          <  2:2255586/101‑1 (MQ=255)
             cgggcACCAACATTCCGGCTGCGGCGCGCGAGATCCTTGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGcc                                                                                         <  1:2418700/101‑1 (MQ=255)
             cgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGGT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGcc                                                                                         >  2:3338757/1‑101 (MQ=255)
             cgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.ggggggggggggc                                                                                         >  2:5118239/1‑94 (MQ=255)
             cgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGcc                                                                                         <  2:4651110/101‑1 (MQ=255)
               ggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.ggggggggggggccc                                                                                       >  2:1076523/1‑92 (MQ=255)
                 cACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGcccccc                                                                                     <  1:3798306/101‑1 (MQ=255)
                    cAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGTGTTGGCCTCCgtc                                                                                  >  1:3686882/1‑99 (MQ=255)
                    cAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGcc                                                                                  <  2:1118646/101‑1 (MQ=255)
                     aaCATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCt                                                                                 <  1:4957521/101‑1 (MQ=255)
                     aaCATCCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGGGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCt                                                                                 <  1:3938855/101‑1 (MQ=255)
                        atcccGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTg                                                                              <  2:3450496/98‑1 (MQ=255)
                        aTTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTg                                                                              <  2:2665385/101‑1 (MQ=255)
                         ttCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCAGGGGGGGGGCGGCCCCCCCCCttgg                                                                              >  2:4893178/1‑99 (MQ=255)
                           ccGGCGGCGGCGGGCGAGATCCTGGATCTGATCCGCGCGGCGGGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGcc                                                                           <  1:5572436/101‑1 (MQ=255)
                           ccGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGcc                                                                           <  2:1929949/101‑1 (MQ=255)
                             ggcggcggcGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGcccc                                                                         <  1:4892005/101‑1 (MQ=255)
                               cggcggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCCCCGTTGGcccccc                                                                       >  2:242765/1‑101 (MQ=255)
                                   ggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGa                                                                   >  2:129035/1‑101 (MQ=255)
                                    gcgcgcGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGat                                                                  <  2:4161114/101‑1 (MQ=255)
                                    gcgcgcGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGat                                                                  <  2:589382/101‑1 (MQ=255)
                                    gcgcgcGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCCTTGGCCCCCCGGGat                                                                  >  2:5377637/1‑101 (MQ=255)
                                     cgcgcgAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGata                                                                 <  1:1839510/101‑1 (MQ=255)
                                       cgggagaTCCTGGATCTGATCCGCGCGGCGGGGGGCT.GACCCGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATAtt                                                               <  1:1411089/98‑1 (MQ=255)
                                       cgcgAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATAtt                                                               <  1:814385/101‑1 (MQ=255)
                                        gcgAGATCCGGGATCTGACCCGCGGGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATAttt                                                              <  2:2091556/101‑1 (MQ=255)
                                          gagaTCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCAGGGGGGGGGGGGCCCCCCCCCTTTGGgccccccgctgttttt                                                             >  1:4872519/1‑91 (MQ=255)
                                          gagaTCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTg                                                            >  2:1703092/1‑101 (MQ=255)
                                          gagaTCCTGGATCTGACCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTg                                                            <  1:765085/101‑1 (MQ=255)
                                                  ggATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGGTGCCCCCCGCCTTTGGTCCCCCGGGAGATTTTGGGCAAGAt                                                    >  2:1811205/1‑101 (MQ=255)
                                                  ggATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGGCATTTTGGGCAAGAt                                                    >  1:1321300/1‑101 (MQ=255)
                                                  ggATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGAt                                                    <  2:3276381/101‑1 (MQ=255)
                                                   gATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGAtt                                                   <  2:5418573/101‑1 (MQ=255)
                                                    gtcggATCCGCGCGGGGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTa                                                  <  2:3942071/97‑1 (MQ=255)
                                                     tCTGATCCGCGGGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTAc                                                 <  2:5645101/101‑1 (MQ=255)
                                                     tCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTAc                                                 <  2:4465313/101‑1 (MQ=255)
                                                     tCTGACCCGCGCGGCGGGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTAc                                                 <  1:3054366/101‑1 (MQ=255)
                                                      cggaccggcggggggcgGGGCT.GACACGGGGAGGCTGCCCGGGGGCA.GGGGGGCACTCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACt                                                <  2:4037286/94‑1 (MQ=255)
                                                      cTGATCCGCGGGGCGCGGGGCT.GACACGGGGCAGCGGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACt                                                <  2:3691182/101‑1 (MQ=255)
                                                      cTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCGGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACt                                                <  2:5351991/101‑1 (MQ=255)
                                                       tGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTg                                               <  2:2177330/101‑1 (MQ=255)
                                                       tGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCCGCCCCCCGCCCTTTGCCCCCCCGGCGTTTTTGGGCAAGATTACta                                               >  1:1713654/1‑100 (MQ=255)
                                                        gATCCGCGCGGCGGGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGa                                              <  2:3576533/101‑1 (MQ=255)
                                                        gATCCGCGCGGCGCGGGGCT.GACAGGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGa                                              <  2:620839/101‑1 (MQ=255)
                                                        gATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGGGCTGCCCCCCCCCCTTGGCCCCCCCGGGGATTTTTGGCAAGATTACTGa                                              >  1:289304/1‑101 (MQ=255)
                                                         aTCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGaa                                             >  1:539882/1‑101 (MQ=255)
                                                           ccGCGCGGCGGGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGaaaa                                           <  1:5065264/101‑1 (MQ=255)
                                                           ccGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGGGCTGCCCCCCCCCCTTTGCCCCCCCCGGTATTTTTGGCCAGATTAATGaaaa                                           >  1:2213247/1‑101 (MQ=255)
                                                               gcggcgcggGGCT.GACAGGGGGCAGCTGACACGGGGCG.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGgc                                       <  2:3871026/101‑1 (MQ=255)
                                                                gggggcggGGGTGGACACGGGGCAGCGGACGCGGGGCA.GGGGGGCGCGGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGgc                                       <  1:1068672/97‑1 (MQ=255)
                                                                cggcggggGGCT.GACAGGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGgcg                                      <  2:1581108/101‑1 (MQ=255)
                                                                 ggcggggGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGgcgc                                     <  1:1140127/101‑1 (MQ=255)
                                                                   ggcggGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTc                                   <  2:3786150/100‑1 (MQ=255)
                                                                     cggGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCcg                                 >  1:596526/1‑101 (MQ=255)
                                                                      ggggCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCcgc                                <  2:4653534/101‑1 (MQ=255)
                                                                      ggggCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCcgc                                <  2:4276183/101‑1 (MQ=255)
                                                                      ggggCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCcgc                                <  1:3497154/101‑1 (MQ=255)
                                                                         gct.gacacggggcagctgacacggggca.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGGGATCCGCGAg                             >  1:1796963/1‑101 (MQ=255)
                                                                          gttcacgggggtcagcggacccggggca.gGGGGGCGCTTCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAg                             <  1:4034022/93‑1 (MQ=38)
                                                                            tgaccgggggcagctgacacggggca.gGGGGGCCCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAggg                           <  1:5427241/95‑1 (MQ=255)
                                                                             gacacggggcagctgacacggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCCCGAgggg                          >  2:702232/1‑101 (MQ=255)
                                                                              acggggggcagctgacggggggcg.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGgc                         <  1:4452487/97‑1 (MQ=255)
                                                                              acagggggcagctgaccgggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGgc                         <  2:2596910/101‑1 (MQ=255)
                                                                               cgcggggcagctgccgcggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGgcg                        <  1:4322320/99‑1 (MQ=255)
                                                                                  ggggcggctgacgggggggc.gGTGGGCCCTCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCt                     <  1:3251017/101‑1 (MQ=255)
                                                                                  ggggcagctgacccggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCt                     <  2:5147038/101‑1 (MQ=255)
                                                                                  ggggcagctgacacggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCt                     >  2:4727484/1‑101 (MQ=255)
                                                                                   gggcagctgacacggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCtt                    >  1:5724710/1‑101 (MQ=255)
                                                                                   gggcagctgacacggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCtt                    >  2:5447317/1‑101 (MQ=255)
                                                                                   gggcagctgacacggggca.gGGGGGCGCGCCCCCCCCCCTCTGGCCCCCCGGGAGATTTTGGGGAAGATTACTGAAAAAGGCGCTCCGCGCGGGGCGCCtt                    >  2:5157152/1‑101 (MQ=255)
                                                                                    ggcagctgacagggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCttt                   <  1:4651540/101‑1 (MQ=255)
                                                                                    ggcagcggacagggggca.gGGGGGCCCCCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCttt                   <  1:2089863/101‑1 (MQ=255)
                                                                                      gagccgaccgggggca.gGGGGGCCCTCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCCCGAGGGGCGCCTTTGa                 <  1:5054533/91‑1 (MQ=255)
                                                                                      cagctgacacggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGa                 <  2:279792/101‑1 (MQ=255)
                                                                                       agctgacacggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGAc                >  2:3773145/1‑101 (MQ=255)
                                                                                         ctgacacggggca.gGGGGGCGCTCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCt              <  1:4122851/101‑1 (MQ=255)
                                                                                          tgacacggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGGGCTCCGCGAGGGGCGCCTTTGGCCtt             >  2:5514837/1‑101 (MQ=255)
                                                                                             cgcggggcg.gGGGGGCGCCCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgt          <  1:1985355/99‑1 (MQ=255)
                                                                                              cgggggca.gGGGGGCGCGCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgtg         <  2:4394288/99‑1 (MQ=255)
                                                                                              acggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgtg         <  1:5294575/101‑1 (MQ=255)
                                                                                              acggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgtg         <  1:2828012/101‑1 (MQ=255)
                                                                                                ggggca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGGTATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgtggg       >  2:4210116/1‑99 (MQ=255)
                                                                                                ggggca.gGGGGGCGCTGCCCCCCCCCTTTGGCCCCCCCGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCCTTGACCTTTGTGCg       >  1:416704/1‑101 (MQ=255)
                                                                                                 gggcg.gGGGGGCGCTCCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCgg      <  2:4935151/101‑1 (MQ=255)
                                                                                                    ca.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGc   >  1:4233727/1‑101 (MQ=255)
                                                                                                     a.gGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  1:5495811/1‑101 (MQ=255)
                                                                                                      |                                                                                                    
GCGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCT.GACACGGGGCAGCTGACACGGGGCA.GGGGGGCGCTGCCCCCCGCCTTTGGCCCCCCGGGATATTTTGGGCAAGATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCT  >  NC_014034/1000879‑1001079

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