Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,209,679 1 bp→GG 100% intergenic (+11/‑260) cobO → / → cobQ1 cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,209,6781.→G100.0% inf 428intergenic (+10/‑261)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (117/212):  ref base (65/34):  total (182/246)
Fisher's exact test strand distribution p-value = 1.39e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01
*NC_0140342,209,6790C→G100.0% 256.5 428intergenic (+11/‑260)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (121/213):  ref base (61/33):  total (182/246)
Fisher's exact test strand distribution p-value = 9.69e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.66e-01

Only 100 of 492 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780
                                                                                                     ||||                                                                                                      
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                                        >  2:2618127/1‑101 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                                        >  2:3354393/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g.ggg                                                                                                    <  1:1006365/101‑2 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g.gggg                                                                                                   >  2:7991246/1‑99 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g.gggg                                                                                                   >  1:7768451/1‑99 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..G.Ggggcgac                                                                                               >  2:4288657/1‑98 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGcccc                                                                                            >  1:7045533/1‑101 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGcccc                                                                                            <  2:8809605/101‑1 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGccccc                                                                                           <  1:5146184/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCg                                                                                          <  1:1887213/101‑1 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCgg                                                                                         >  2:6847975/1‑101 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCgg                                                                                         <  1:785631/101‑1 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGa                                                                                        >  1:4645725/1‑101 (MQ=255)
                     aCCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCGGAcccc                                                                                    <  1:6087440/101‑1 (MQ=255)
                     aCCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCGGAcccc                                                                                    <  1:725528/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGAccccc                                                                                   >  1:313507/1‑101 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCg                                                                                  <  1:5719515/101‑1 (MQ=38)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCgg                                                                                 <  2:2334774/101‑1 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCggg                                                                                >  2:6291696/1‑101 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCggg                                                                                <  2:6603731/101‑1 (MQ=38)
                          gaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGGCCCCCGGGa                                                                               >  1:1748117/1‑101 (MQ=255)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGat                                                                              <  2:5832811/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGat                                                                              <  1:2695979/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGGTCTGATCTCAAGGCCGG..G.GGAGCTGCCCCCGGACCCCCggggt                                                                              >  2:2972108/1‑99 (MQ=255)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGata.t                                                                           >  2:1172217/1‑101 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGata.t                                                                           >  1:7586821/1‑101 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.tt                                                                          <  2:9086628/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.tt                                                                          <  1:5543922/101‑1 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.ttt                                                                         <  1:426940/101‑1 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.ttt                                                                         <  2:5605599/101‑1 (MQ=38)
                                  aGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTtgtg                                                                      >  1:6583823/1‑101 (MQ=37)
                                  aGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGACCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGGGCTGCCCCCGGACCCCCGGGATA.TTtgtg                                                                      >  1:8210876/1‑101 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGc                                                                    <  2:7536472/101‑1 (MQ=37)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCa                                                                   <  1:6860544/101‑1 (MQ=37)
                                      gAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCaa                                                                  >  2:6452168/1‑101 (MQ=37)
                                      gAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCaa                                                                  <  1:4068739/101‑1 (MQ=37)
                                      gAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCAGACCCCCGGGATA.TTTGTGGCaa                                                                  >  2:6196791/1‑101 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAg                                                                 <  2:1081670/101‑1 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGAt                                                               >  2:4977379/1‑101 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGTAAGATg                                                              <  2:1737304/101‑1 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATg                                                              <  2:2742826/101‑1 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGa                                                             >  1:8928543/1‑101 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGa                                                             <  1:3444831/101‑1 (MQ=25)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGaa                                                            <  1:9122390/101‑1 (MQ=25)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGaa                                                            <  1:6981901/101‑1 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACa                                                         >  2:3473020/1‑101 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACa                                                         >  1:6200571/1‑101 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACa                                                         >  2:2013796/1‑101 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACag                                                        >  1:587347/1‑101 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACag                                                        <  2:3200921/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACaga                                                       <  2:3745435/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACaga                                                       <  1:542916/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACaga                                                       >  1:6379611/1‑101 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACagag                                                      <  1:488252/101‑1 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGAGGCAAGATGAAACagag                                                      <  2:5369898/101‑1 (MQ=21)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGc                                                     >  1:4492696/1‑101 (MQ=25)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCg                                                    >  2:605788/1‑101 (MQ=25)
                                                     ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGt                                                   >  1:6700094/1‑101 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGtt                                                  <  2:5474910/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGtt                                                  <  2:7713291/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTAGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGtt                                                  <  1:3275431/101‑1 (MQ=21)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTa                                                 <  2:1878997/101‑1 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAg                                                >  2:3930623/1‑101 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAg                                                >  2:7376895/1‑101 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAg                                                <  1:5583615/101‑1 (MQ=25)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAgg                                               >  2:6717101/1‑101 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAggg                                              >  2:4710545/1‑101 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAggg                                              <  2:2868496/101‑1 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAggg                                              >  1:8630615/1‑101 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCCGGATA.TTTGTGGCAAGATGAAACAGAGCGGTAgag                                              >  2:102703/1‑99 (MQ=21)
                                                           caTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGa                                             >  1:1396884/1‑101 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGGTCTCAAGGCC.G..GCGGCGCTGGCAGCGGCCCCACGGGTTA.TTTGAGGCAAGATGAAACAGAGCGAGAGGGAGTTg                                         >  2:1595449/1‑101 (MQ=16)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG.G.GGGGGTGGCCCCGGACCCCCCGGGTATTTTGTGGGAGAATAAAACA.AGCGGTAgggggttg                                         >  2:3047173/1‑96 (MQ=255)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                                      <  1:7496799/101‑1 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCa                                  <  1:2561527/101‑1 (MQ=25)
                                                                       gCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                                 <  2:8610457/101‑1 (MQ=25)
                                                                       gCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                                 >  1:882051/1‑101 (MQ=25)
                                                                         cGGGGGTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGa                               >  1:846048/1‑101 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc                             <  2:5243102/101‑1 (MQ=25)
                                                                            gggTCGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCg                            >  1:2770281/1‑101 (MQ=25)
                                                                             ggTCGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCgg                           <  2:8720393/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCa                        <  2:4080818/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCa                        >  1:5022696/1‑101 (MQ=25)
                                                                                 gAGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa                       <  1:1172235/101‑1 (MQ=25)
                                                                                 gAGTTCTGATCTC.AAGGC.CGGG.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa                       >  1:7362438/1‑101 (MQ=25)
                                                                                  aGTTCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAg                      <  2:7821864/101‑1 (MQ=25)
                                                                                    ttCTGATCTCAAGGCCGG..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt                    <  1:3629084/101‑1 (MQ=25)
                                                                                      cTGATCTCAAGGCC.GG.G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTgg                  >  2:8168371/1‑101 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG.G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGc                 >  2:7254559/1‑101 (MQ=25)
                                                                                              tcaaggccggg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAacac            <  2:7054572/91‑1 (MQ=255)
                                                                                                  ggccggg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTa        <  1:1820002/95‑1 (MQ=255)
                                                                                                   gccggg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTAc       <  1:7129219/96‑1 (MQ=255)
                                                                                                   gccggg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTAc       <  1:3765381/96‑1 (MQ=255)
                                                                                                    ccggg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACCCGGGACa      >  2:8383809/5‑101 (MQ=255)
                                                                                                    ccggg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACa      <  2:2517886/97‑1 (MQ=255)
                                                                                                    atctg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACa      <  1:2090709/97‑1 (MQ=255)
                                                                                                     cggg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACag     <  1:1601074/98‑1 (MQ=255)
                                                                                                      ggg.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  2:1356385/3‑101 (MQ=255)
                                                                                                        g.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaggga  <  1:3878970/101‑1 (MQ=255)
                                                                                                        g.ggCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaagat  >  2:2090647/1‑97 (MQ=255)
                                                                                                     ||||                                                                                                      
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..G.GGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780

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