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% 133.7 157intergenic (+10/‑261)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (26/92):  ref base (26/13):  total (52/105)
Fisher's exact test strand distribution p-value = 8.98e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.20e-01
*NC_0140342,209,6790C→G100.0% 93.1 157intergenic (+11/‑260)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (28/92):  ref base (24/13):  total (52/105)
Fisher's exact test strand distribution p-value = 6.23e-06
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.41e-01

Only 100 of 193 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780
                                                                                                     |||||                                                                                                    
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...ggg                                                                                                    >  1:4224579/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...gggg                                                                                                   <  1:4298825/101‑2 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...gggg                                                                                                   >  2:2745169/1‑100 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...ggggg                                                                                                  <  2:1811858/101‑3 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...ggggg                                                                                                  >  1:584919/1‑99 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...ggggg                                                                                                  >  1:229015/1‑99 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggc                                                                                                 <  1:3287080/101‑1 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggc                                                                                                 <  2:3918661/101‑1 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...gggggggcc                                                                                              >  1:1140512/1‑95 (MQ=255)
          ccGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.C...gggggggggg                                                                                             >  1:2429738/1‑94 (MQ=255)
          ccGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggcgctgc                                                                                            >  2:4464767/1‑96 (MQ=255)
          ccGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggcgctgc                                                                                            >  2:314818/1‑96 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggcgctgcc                                                                                           >  1:3817821/1‑95 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGG.GGGCGCTGcccc                                                                                           <  1:2719522/101‑1 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGG.GGGCGCTGccccc                                                                                          <  1:1503260/101‑1 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCgg                                                                                        >  2:3718537/1‑101 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCgg                                                                                        <  1:2682879/101‑1 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCggc                                                                                       >  1:2602060/1‑100 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGa                                                                                       <  2:4653768/101‑1 (MQ=255)
                  gTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGCAc                                                                                      <  1:3445745/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGAccccc                                                                                  <  2:3958112/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGAccccc                                                                                  <  1:2749778/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCg                                                                                 <  1:2729026/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCg                                                                                 >  2:718320/1‑101 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCg                                                                                 <  1:1003591/101‑1 (MQ=38)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCCGCCCCCCgg                                                                                >  2:1952760/1‑101 (MQ=255)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCggg                                                                               <  2:584858/101‑1 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCggg                                                                               <  1:1398446/101‑1 (MQ=38)
                          gaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGa                                                                              <  2:2536986/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGat                                                                             <  1:3403904/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGCCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCGGACCCCCGGGat                                                                             >  1:4191184/1‑101 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGatat                                                                           <  1:2573253/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATAtt                                                                          >  2:1083516/1‑101 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATAttt                                                                         <  2:917687/101‑1 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATAttt                                                                         <  1:3718812/101‑1 (MQ=38)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTtgt                                                                       >  1:4403621/1‑101 (MQ=37)
                                  aGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGGC.CGCGGGGGGCGGCCCCCGGCcccccccggtttttg                                                                        >  2:4026861/1‑91 (MQ=255)
                                   ggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCGGACCCCCGGGGTATTTGTgg                                                                     >  2:3479254/1‑101 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGc                                                                    >  2:1376382/1‑101 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCACCAGACCTCCGGGATATTTGTGGc                                                                    >  2:262681/1‑101 (MQ=255)
                                      gAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCGGACCCCCGGGGTATTTGTGGCaa                                                                  >  1:2830833/1‑101 (MQ=37)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGa                                                                <  2:1155648/101‑1 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAt                                                               <  2:2942606/101‑1 (MQ=25)
                                          ctcccGTTCCGGGCAGGCATCAAGGGCCAGTCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                                              <  2:3864780/99‑1 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCCGGGTATTTGTGGccaagtga                                                             >  2:328863/1‑94 (MQ=25)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCCGCCCCCGGACCCCCGGGATATTTGTGGCAAGACGaa                                                            >  1:630074/1‑101 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                                           <  1:3029147/101‑1 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                                          <  1:3345375/101‑1 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACa                                                         <  1:1705176/101‑1 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                                        <  1:2273547/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                       <  2:2417470/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                       <  2:2104184/101‑1 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGGTGCCCCCGGGCCCCCGGGGTATTTTTGGCAAGATGAAACagag                                                      >  1:403680/1‑101 (MQ=255)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                                     <  2:1581226/101‑1 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGCTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                                     <  1:3997815/101‑1 (MQ=25)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCGGACCCCCGGGCAAATTGGGGCCAGGTGAGACAGAGCg                                                    >  1:1222847/1‑101 (MQ=255)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCg                                                    <  1:4114484/101‑1 (MQ=25)
                                                     ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGGCCCCCGGGAAATTTGTGGCAAGATGAAACAGAGCGt                                                   >  1:3518028/1‑101 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                                  <  2:4115067/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                                  <  2:3830781/101‑1 (MQ=25)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAgg                                               >  1:2292590/1‑101 (MQ=25)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGGC.CG..GGGCGCTGCCCCCGGACCCCCGGGACATTTGTGGCAAGATTAAACATAGCGTTAAGGAttt                                          >  2:608704/1‑101 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                                       <  1:2200406/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGG.GGGGGCTGCCCCCCGACCCCCCGGGTATTTTTGGGAAGATGGAAAAGAGGGGTAGGGGttttgg                                       >  2:1860970/1‑98 (MQ=255)
                                                                   gCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGGGGCGCCCCCCGGCCCCCCGGGGAATTT.TGGCGAGAATAAAAAAAGCGTTACGGAGTTgggga                                     >  1:3754447/1‑98 (MQ=255)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCCGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTTGTAAcc                                   >  1:778155/1‑101 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCGCACGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAGCCa                                  <  2:1481456/101‑1 (MQ=21)
                                                                       gCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                                 <  1:2662954/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt                              <  1:1971111/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt                              <  1:4260316/101‑1 (MQ=25)
                                                                           ggtgTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc                             <  2:3095023/98‑1 (MQ=21)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc                             <  2:4277392/101‑1 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc                             <  2:2709089/101‑1 (MQ=25)
                                                                            gggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCg                            <  1:1905910/101‑1 (MQ=25)
                                                                            gggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCg                            <  1:916656/101‑1 (MQ=25)
                                                                               tCGAGTTCTGATCTCAACGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGc                         <  2:4190896/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCCGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGGAACCATGATCGGGCa                        >  1:3062235/1‑101 (MQ=21)
                                                                                cGAGTTCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCGGGGTATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCa                        >  1:3880997/1‑101 (MQ=25)
                                                                                 gAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa                       <  2:3817536/101‑1 (MQ=25)
                                                                                  aGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAg                      <  1:2851799/101‑1 (MQ=25)
                                                                                   gTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc                     <  1:1627159/101‑1 (MQ=25)
                                                                                   gTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc                     <  2:2068650/101‑1 (MQ=25)
                                                                                   gTTCTGATCTC.AAGGC.CGG.GGGGGCTGCCCCCGGACCCCCGGGGTATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc                     >  1:1420250/1‑101 (MQ=25)
                                                                                   gTTCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc                     >  2:1538136/1‑101 (MQ=25)
                                                                                    ttCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt                    <  1:99502/101‑1 (MQ=25)
                                                                                     tCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg                   >  2:2252888/1‑101 (MQ=25)
                                                                                     tCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg                   >  2:4639175/1‑101 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG..GGGCGCTGCCCCCGGGCCCCCGGGATATTTGTGGCAAGATTAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAACTgga                 >  2:3643447/1‑100 (MQ=21)
                                                                                                  aggccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGt         <  2:1843199/94‑1 (MQ=255)
                                                                                                  aggccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGt         <  2:2786853/94‑1 (MQ=255)
                                                                                                   ggccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTa        <  1:328895/95‑1 (MQ=255)
                                                                                                    gccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTAc       <  2:4216702/96‑1 (MQ=255)
                                                                                                    gccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTAc       <  1:1776012/96‑1 (MQ=255)
                                                                                                     ccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTaaa      >  1:1345856/5‑99 (MQ=255)
                                                                                                      cgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACag     <  2:1414528/98‑1 (MQ=255)
                                                                                                       gggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGGTTTGGTAAACATGATCGGGCAAGCTGGCAACACGGTACagg    >  2:879730/3‑101 (MQ=255)
                                                                                                       gggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  2:4543732/3‑101 (MQ=255)
                                                                                                       gggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  2:639726/3‑101 (MQ=255)
                                                                                                         gggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaggga  >  2:2516387/1‑101 (MQ=255)
                                                                                                         gggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaggga  >  1:4582050/1‑101 (MQ=255)
                                                                                                     |||||                                                                                                    
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG