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% 199.5 198intergenic (+10/‑261)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (35/107):  ref base (29/27):  total (64/134)
Fisher's exact test strand distribution p-value = 3.66e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00
*NC_0140342,209,6790C→G100.0% 144.1 198intergenic (+11/‑260)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (36/110):  ref base (28/24):  total (64/134)
Fisher's exact test strand distribution p-value = 2.32e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

Only 100 of 234 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780
                                                                                                     ||||                                                                                                      
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                                        >  2:2501847/1‑101 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggg                                                                                                      <  1:1614595/101‑1 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggg                                                                                                      >  2:4169728/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGACGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                                     >  2:3092194/1‑100 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggggg                                                                                                    >  1:1747946/1‑99 (MQ=255)
     ggcggcCGATCTAGTAACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggc                                                                                                  <  1:2373463/101‑1 (MQ=255)
      gcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcg                                                                                                 <  1:226166/101‑2 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGCAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgct                                                                                               >  1:5162341/1‑98 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgct                                                                                               <  2:5786382/101‑4 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgct                                                                                               <  1:2039590/101‑4 (MQ=255)
         gctgagcTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTGCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgctg                                                                                              <  2:4588401/95‑5 (MQ=255)
         gcCGATCTGGTTTCGGAGACGGGTGAGGTGAAGCTCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgctg                                                                                              <  1:4533785/101‑5 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgctg                                                                                              >  2:2540087/1‑97 (MQ=255)
          ccGATCTTGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCACGGGC.G..gggggggcgc                                                                                               >  2:1194602/1‑94 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.C..ggggggggggc                                                                                              >  1:2939715/1‑93 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgctggc                                                                                            >  2:2399650/1‑95 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGccc                                                                                             >  1:1313477/1‑101 (MQ=255)
               ttAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCg                                                                                          <  2:1827304/100‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCg                                                                                          >  2:4643837/1‑101 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCgg                                                                                         >  2:4489670/1‑101 (MQ=255)
                tAGTGACCGAGATGGATCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCgg                                                                                         <  2:640700/101‑1 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTGCCCCCGGAccc                                                                                     <  2:1855499/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGAccccc                                                                                   <  1:3169292/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCg                                                                                  <  2:1478742/101‑1 (MQ=38)
                          gaTGGGTCAGGTGAAGCGCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGa                                                                               <  2:3274262/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGat                                                                              <  1:425235/101‑1 (MQ=38)
                            tGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGata                                                                             <  1:3698385/101‑1 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.GGGTGGCCCCCGACCCCCggggtgt                                                                            >  1:5622810/1‑97 (MQ=255)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGatat                                                                            <  1:2896640/101‑1 (MQ=38)
                              gggcAGGTGAAGCCCCCGTTCCGGGCGGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATAtt                                                                           <  1:5296803/98‑1 (MQ=25)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATAtt                                                                           <  2:1277438/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATAtt                                                                           >  1:343840/1‑101 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATAttt                                                                          <  1:5205265/101‑1 (MQ=38)
                                tCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTtg                                                                         <  1:1256643/101‑1 (MQ=37)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTtgt                                                                        <  1:2683434/101‑1 (MQ=37)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTtgt                                                                        <  1:1069589/101‑1 (MQ=37)
                                   ggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTgg                                                                      <  2:1339461/101‑1 (MQ=37)
                                   gatGAAGCACCCGTTCCGGGCTGTCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTgg                                                                      <  2:5705357/99‑1 (MQ=25)
                                      gAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCaa                                                                   >  1:2588544/1‑101 (MQ=37)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAg                                                                  <  2:629613/101‑1 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAg                                                                  <  1:3664521/101‑1 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGa                                                                 >  2:833570/1‑101 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.GGCTGCCCCCCGGCCCCCCGGGGATTTTTGGGCAGAt                                                                >  1:477968/1‑101 (MQ=255)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAt                                                                >  2:3362457/1‑101 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                                               >  2:2785660/1‑101 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                                               <  1:2533607/101‑1 (MQ=25)
                                            cccGTTCCGTGCAGGCATCAAGGGCCTGCCGGGGGTCGAGTTCTGATCGCAAGGCC.G..GGGCCTCTGCCCCCGGACCCCGGGGATATTTGTGGCAAGATGaa                                                             <  2:3946524/101‑1 (MQ=21)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaa                                                             <  1:4558394/101‑1 (MQ=25)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaa                                                             <  2:399647/101‑1 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                                            <  2:4869070/101‑1 (MQ=25)
                                               gtcccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACa                                                          <  1:1473111/98‑1 (MQ=21)
                                                tgccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                                         <  2:406797/99‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGGTGCCCCCGGGCCCCCGGGGGATTTGTGGCAAGATGaaaaaaa                                                        >  1:2485280/1‑97 (MQ=21)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                        <  2:4321425/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                        <  2:3724828/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                        <  2:4797315/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                        <  2:1970683/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                        <  1:2786316/101‑1 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACagag                                                       >  2:724195/1‑101 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                                      >  2:5176689/1‑101 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                                      >  1:693723/1‑101 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                                      <  2:1012128/101‑1 (MQ=25)
                                                     ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGt                                                    <  2:1547076/101‑1 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.g                                                <  2:4449519/101‑1 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCCGGGAATTTGTGGCAAGATGAAACAGAGCGTTA.ggg                                              >  1:5233832/1‑101 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCCGCCCCCGGGCCCCCCGGATATTTGTGGCCAGATGAAACAGAGCGTTA.GGGAt                                            >  1:3572145/1‑101 (MQ=37)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGAttt                                          <  1:2679385/101‑1 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTTGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTg                                         <  2:5198397/101‑1 (MQ=21)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.GGCGGCCCCCGGACCCCCGGGAGATTTGTGGGGAGATGAAACAGAGCGCTA.GGGGTTTg                                         >  1:3606128/1‑101 (MQ=21)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTg                                         >  2:2718136/1‑101 (MQ=25)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAAcc                                   <  2:2294565/101‑1 (MQ=25)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAAcc                                   <  2:2250574/101‑1 (MQ=25)
                                                                     cAGACGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAAcc                                   <  1:3092368/101‑1 (MQ=21)
                                                                      aGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCa                                  <  2:4541633/101‑1 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGAtagatcggaag                                  >  2:1697249/1‑91 (MQ=21)
                                                                       gCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCAt                                 <  2:4907954/101‑1 (MQ=25)
                                                                        cggggggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTTGTAACCATg                                <  2:2640341/99‑1 (MQ=21)
                                                                         cGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.GGCCGCCCCCCGGCCCCCCGGGTATTTGTGGGAAGATGAAAAAAAGGGGTAGGGGGTTTGGTACcctga                                >  2:4091135/1‑98 (MQ=255)
                                                                          gggggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGAt                              <  1:3729882/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCCGCCCCCGGAAGCCCACGGTATTTGTGGCAAGACTAAACAGAGAGCTA.GGGATTTGGTAACCATTAt                              >  1:3897847/1‑101 (MQ=255)
                                                                          gggggTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGAt                              >  1:3503289/1‑101 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGTTc                             <  1:3295726/101‑1 (MQ=21)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATc                             <  1:5451206/101‑1 (MQ=25)
                                                                           ggggTCGACTTCTGATCTCAAGGCCGG..GGG.CGCGGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATc                             <  1:4996253/101‑1 (MQ=21)
                                                                            gggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCg                            <  2:3066006/101‑1 (MQ=25)
                                                                             ggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCgg                           <  2:5622518/101‑1 (MQ=25)
                                                                             ggTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATcag                           >  2:3218442/1‑99 (MQ=25)
                                                                                 gAGTTCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCaa                       <  1:1875652/101‑1 (MQ=25)
                                                                                  aGTTCTGATCTC.AAGGC.CGGGGG.GGCTGCCCCCGGACCCCCCGGGAATTTTTGGGAAAATGAAAAAGAGCGTTA.GGGGTTTTGTAAACATGATTggggcaa                      >  2:5151356/1‑97 (MQ=12)
                                                                                    ttCTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCt                    <  1:833638/101‑1 (MQ=25)
                                                                                      cTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTgg                  <  1:4644113/101‑1 (MQ=25)
                                                                                      cTGATCTCAAGGCCGG..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTgg                  >  2:3588065/1‑101 (MQ=25)
                                                                                              tcaaggccggggg.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAacac            >  2:806836/11‑101 (MQ=255)
                                                                                                aaggccggggg.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACgg          <  1:1207631/93‑1 (MQ=255)
                                                                                                 aggccggggg.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGt         <  2:2908618/94‑1 (MQ=255)
                                                                                                      ggggg.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  2:4657653/3‑101 (MQ=255)
                                                                                                      ggggg.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  2:3688319/3‑101 (MQ=255)
                                                                                                      ggggg.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  1:4456376/3‑101 (MQ=255)
                                                                                                      ggggg.CGCTGCCCCCGGACCCCCGGGATATTTGAGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    <  1:2151431/99‑1 (MQ=255)
                                                                                                       gggg.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaggg   <  1:3826112/100‑1 (MQ=255)
                                                                                                     ||||                                                                                                      
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGG.CGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG