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

Only 100 of 200 total aligned reads displayed.

GCGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCT  >  NC_014034/1000879‑1001079
                                                                                                     |                                                                                                        
gCGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGgctgacacggggcagctgacacggggca                                                                                                           <  2:2878624/101‑1 (MQ=255)
 cGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGTCGCGCGAGATCCTGGACCTGATCCGCGCGGCGCGGGGctgacacggggcagctgacacggggca.g                                                                                                         >  2:2817784/1‑101 (MQ=255)
 cGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGctgacacggggcagctgacacggggca.g                                                                                                         <  2:301155/101‑1 (MQ=255)
  gATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.gg                                                                                                        <  1:378561/101‑1 (MQ=255)
   aTCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.ggg                                                                                                       >  2:773835/1‑101 (MQ=255)
     cTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.ggggg                                                                                                     >  1:1226478/1‑101 (MQ=255)
     cTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.ggggg                                                                                                     >  2:2124373/1‑101 (MQ=255)
      ttCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.gggggg                                                                                                    >  2:1424750/1‑101 (MQ=255)
       tCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGgc                                                                                                   <  1:2514195/101‑1 (MQ=255)
        cgggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGgcg                                                                                                  <  2:1868351/101‑1 (MQ=255)
         gggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.ggggggggg                                                                                                 >  1:3567926/1‑98 (MQ=255)
          ggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGgcgcg                                                                                                >  1:1972051/1‑100 (MQ=255)
          ggccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGACCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.t                                                                                               <  1:150467/101‑1 (MQ=255)
           gccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.g                                                                                             <  2:1783934/101‑1 (MQ=255)
           gccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.g                                                                                             <  2:2226326/101‑1 (MQ=255)
            ccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Gc                                                                                            <  1:1809545/101‑1 (MQ=255)
            ccgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Gc                                                                                            <  2:4206698/101‑1 (MQ=255)
             cgggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.gggggggggcggg                                                                                             >  2:23994/1‑94 (MQ=255)
              gggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.gggggggggcgccc                                                                                            >  2:3632610/1‑93 (MQ=255)
              gggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Gccc                                                                                          <  2:614588/101‑1 (MQ=255)
              gggcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Gccc                                                                                          >  2:4121632/1‑101 (MQ=255)
                gcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGGGG.T.Gccccc                                                                                        >  1:609154/1‑101 (MQ=255)
                gcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Tccccc                                                                                        >  2:1669747/1‑101 (MQ=255)
                gcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Gccccc                                                                                        <  2:1740877/101‑1 (MQ=255)
                gcACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Gccccc                                                                                        <  1:1003908/101‑1 (MQ=255)
                 cACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGGGC.G.Ggccccc                                                                                       >  2:814355/1‑101 (MQ=255)
                 cACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Gcccccc                                                                                       <  1:266730/101‑1 (MQ=255)
                 cACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.Gcccccc                                                                                       <  2:3241812/101‑1 (MQ=255)
                       cATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.ttt                                                                                >  2:3398449/1‑101 (MQ=255)
                        aTTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.CTTg                                                                               >  1:2046090/1‑101 (MQ=255)
                         ttCCGGCGGCGGCGCGGGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTgg                                                                              <  2:344923/101‑1 (MQ=255)
                             ggcggcggcGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGcccc                                                                          <  2:3926878/101‑1 (MQ=255)
                               gggcggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGcccccc                                                                        <  2:1719025/100‑1 (MQ=255)
                               cggcggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGcccccc                                                                        >  2:2022297/1‑101 (MQ=255)
                                ggcggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GACCCCCGCC.TTTGGCCCCCCg                                                                       >  2:3347380/1‑101 (MQ=255)
                                  gggcgCGCGAGATCCTGGATCTGATCCGCGCGGCGGGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCggg                                                                     <  1:3516127/100‑1 (MQ=255)
                                  cggcgCGCGAGATCCGGGATCTGATCCGGGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCggg                                                                     <  1:3379862/101‑1 (MQ=255)
                                     cccgcgAGACCCTGGATCTGCTCCGCGCGGCGCGGTGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGata                                                                  <  2:2387004/99‑1 (MQ=255)
                                         cgAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATAtttt                                                              >  2:2467133/1‑101 (MQ=255)
                                          gagaTCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTg                                                             <  2:261328/101‑1 (MQ=255)
                                           agaTCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.CTTGGCCCCCCGGGATATTTTgg                                                            >  1:1320148/1‑101 (MQ=255)
                                            gaTCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTggg                                                           >  2:2641117/1‑101 (MQ=255)
                                            gaTCCTGGATCCGATCCGCGCGGGGGGGGGCTGACCCGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTggg                                                           <  2:3153975/101‑1 (MQ=255)
                                             acccTGGATGGGACCCGGGGGGGGGGGGGCTGACACGGGGCAGCTGACACGGGGCG.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGc                                                          <  1:2932068/99‑1 (MQ=255)
                                             aTCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGGGC.TGGCCCCCCCCC.CTGGGCCCCCCgggggtctttggg                                                           >  1:2320625/1‑91 (MQ=255)
                                             aTCCGGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGc                                                          <  1:2617217/101‑1 (MQ=255)
                                              tCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCa                                                         <  2:2666822/101‑1 (MQ=255)
                                               ccTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.C.GCGCCCCGGC.CTTGGGGCCCCGGGATATTTTGGGAaa                                                        >  2:3424083/1‑101 (MQ=255)
                                                cTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.G.GGCCCCCGCC.CTTGGCCCCCCCGGGTATTTTGGGCAAg                                                       >  2:2543446/1‑101 (MQ=255)
                                                  ggACCTGATCCGCGCGGCGCGGGGCTGACAGGGGGCAGCAGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.At                                                    <  2:2142307/101‑1 (MQ=255)
                                                   gATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.Att                                                   >  1:718295/1‑101 (MQ=255)
                                                   gATCTGACCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.Att                                                   <  1:2124518/101‑1 (MQ=255)
                                                    aTCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTa                                                  <  2:1938309/101‑1 (MQ=255)
                                                     tCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCAGGGGGGGGGC.G.CCCCCCCCCC.TTTGGCCCCCCGGGGCATTTTGGGCAGG.Ataa                                                  >  1:4228600/1‑99 (MQ=255)
                                                     tCTGACCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTAc                                                 <  1:1856511/101‑1 (MQ=255)
                                                       tGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGGGCGG.GCCCCCCCCC.TTGGGCCCCCCGGGTGTTTTTGGGCaggtttctt                                                 >  1:1663882/1‑93 (MQ=255)
                                                       tGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTg                                               >  1:1195057/1‑101 (MQ=255)
                                                       tGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.CTTGGCCCCCCCGGGTATTTTGGGCAAG.ATTACTg                                               >  1:2281003/1‑101 (MQ=255)
                                                         aTCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGGGC.T.GCCCCCCCCC.TTTGGACCCCCCGGATATTTTTGGGAAG.AGTACTGaa                                             >  1:2829526/1‑101 (MQ=255)
                                                          tCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGaaa                                            <  2:4063509/101‑1 (MQ=255)
                                                           ccGCGGGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.CCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGaaaa                                           <  1:2543637/101‑1 (MQ=255)
                                                           ccGCGCGGCGGGGGGCTGACCCGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.CCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGaaaa                                           <  1:4121893/101‑1 (MQ=255)
                                                              cgcggcgcggGGCTGACACGGGGCAGCTGACACGGGGCAGGGGGGGGGG.G.GCCCCCCCCC.TTTGGCCCCCCGGG.TATTTTTGGGAAGAATTATTAAAAAAg                                         >  2:1450196/1‑101 (MQ=255)
                                                               ggtgcgcggGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGgc                                       <  2:3567673/98‑1 (MQ=255)
                                                               gcggcggggGGCTGACGCGGGGCAGCTGACACGGGGCG.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGgc                                       <  1:859148/101‑1 (MQ=255)
                                                                  gcggggGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCt                                    <  2:1120913/101‑1 (MQ=255)
                                                                   cggggggCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTc                                   <  2:1002625/98‑1 (MQ=255)
                                                                   cgcggGGCTGACACGGGGCAGCTGACACGGGCCA.GGGGGGGGC.G.GCCCCCCGCC.CCTTGCCCCCCCCGGTATTTTTGGCAAG.ATTAATGAAAAAAgcgccc                                   >  1:3752647/1‑99 (MQ=255)
                                                                    ggggggCTGCCCGGGGGCAGCTGACACGGGGGG.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGAAAG.ATTACTGAAAAAGGCGCTcc                                  <  1:797908/99‑1 (MQ=255)
                                                                    gcggGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTcc                                  <  1:2436979/101‑1 (MQ=255)
                                                                     cggGGCTGACACGGGGCAGCTGACACGGGGCAGGGGGGGGGG.G.GCCCCCCCCC.CTTGCCCCCCCGGGGATTTTTGGGGAGG.ATTTTTAaaaaaagggcccc                                  >  1:4098321/1‑93 (MQ=255)
                                                                      ggggCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGGGC.T.GCCCCCCCCC.CTTGGCCCCCCCGGGGATTTTGGGGAAG.GTTACTGAAAAAgagggcccgc                                >  1:240380/1‑92 (MQ=255)
                                                                      ggggCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGG.T.GCCCCCCGCCTTTTGGCCCCCCGGG.GATTTTGGGCAAG.GTTACCGAAAAAGGCGCTCcgc                                >  2:3836007/1‑101 (MQ=255)
                                                                       gggCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCcgcg                               >  2:2072310/1‑101 (MQ=255)
                                                                        ggCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGG.C.GCCCCCCGCC.TTTCGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGa                              >  2:2721445/1‑101 (MQ=255)
                                                                        ggCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGa                              >  2:376196/1‑101 (MQ=255)
                                                                          ctgacacggggcagctgacacggggca.gGGGGGGGC.T.GCGCCCCGCC.GTTGGCCCCCCGGGATATTTTGGGCCAG.ATTACTGAAAAAAGCGCTCCGCGTgg                            >  1:2064735/1‑101 (MQ=255)
                                                                           tgacacggggcagctgacacggggca.gGGGGGGGC.T.GCCCCCCGCC.CTTTGCCCCCCGGGGTATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCCCGAggg                           >  1:1128447/1‑101 (MQ=255)
                                                                           tgacacggggcagctgacacggggca.gGGGGGGGC.C.GCCCCCCCCC.CTTGGCCCCCCCGGGGATTTTGGGGGAG.ATTATTGAAAAAAGGGGTCCGCGCggg                           >  1:133799/1‑101 (MQ=255)
                                                                            gacacggggcagctgacacggggca.gGGGGGGGG.T.GCCCCCCCCC.CTTTGCCCCCCCGGGTATTTTTGGCAAG.GTTTCTGAAAAAAGGGGCCCGCGGgggg                          >  1:3100665/1‑101 (MQ=255)
                                                                             gcagggggcagttgacccgggggg.gGGGGGCCC.T.CCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGgc                         <  1:3217796/97‑1 (MQ=255)
                                                                             acccggggcagctgacgcggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGgc                         <  2:3893531/98‑1 (MQ=255)
                                                                             acacggggcagctgacacggggca.gGGGGGGGC.T.GCCCCCCGCC.TTTTGCCCCCCGGGGTATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAggggg                         >  1:2110753/1‑100 (MQ=255)
                                                                             acacggggcagctgacacggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGgc                         <  2:1320048/101‑1 (MQ=255)
                                                                             acacggggcagctgacacggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGgc                         <  2:11261/101‑1 (MQ=255)
                                                                               gggggggagctgaccgggggcg.gGGGGGCGC.C.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGgcgc                       <  1:448628/99‑1 (MQ=255)
                                                                                 ggggcagctgacacggggca.gGGGGGGGG.T.GCCCCCCCCC.CTTTGCCCCCCCGGGGACTTTTGGCGAG.ATTAATTAAAAAAGGGGCCCCCGCGGGgcggcc                     >  2:3691387/1‑98 (MQ=255)
                                                                                 ggggcagctgacacggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCt                     >  1:3352787/1‑101 (MQ=255)
                                                                                  gggcagctgacacggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCtt                    <  2:888946/101‑1 (MQ=255)
                                                                                   ggcagctgtgacggggcg.gGGGGCCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCttt                   <  1:96669/101‑1 (MQ=255)
                                                                                     cagctgacacggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGa                 >  2:3287545/1‑101 (MQ=255)
                                                                                        ctgacacggggcg.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCt              <  2:2320457/101‑1 (MQ=255)
                                                                                        ctgacacggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCt              >  2:521153/1‑101 (MQ=255)
                                                                                         tgacacggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGAACtt             >  2:4027765/1‑101 (MQ=255)
                                                                                          gacagggggcg.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCttt            <  1:388365/101‑1 (MQ=255)
                                                                                          gacacggggcg.gGGGGGCGC.T.CCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCttt            <  2:1831937/101‑1 (MQ=255)
                                                                                               ggggca.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTtgtggg       >  2:1201314/1‑99 (MQ=255)
                                                                                                  gca.gGGGGGCGC.T.GCCCCTCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCgggg    >  1:1679459/1‑101 (MQ=255)
                                                                                                    a.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  1:2132902/1‑101 (MQ=255)
                                                                                                    a.gGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCt  >  1:1723494/1‑101 (MQ=255)
                                                                                                     |                                                                                                        
GCGATCTTCGGGCCGGGCACCAACATTCCGGCGGCGGCGCGCGAGATCCTGGATCTGATCCGCGCGGCGCGGGGCTGACACGGGGCAGCTGACACGGGGCA.GGGGGGCGC.T.GCCCCCCGCC.TTTGGCCCCCCGGGATATTTTGGGCAAG.ATTACTGAAAAAGGCGCTCCGCGAGGGGCGCCTTTGACCTTTGTGCGGGGCT  >  NC_014034/1000879‑1001079

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG