Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,209,677 G→C 100% intergenic (+9/‑262) 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,6770G→C100.0% 160.7 202intergenic (+9/‑262)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (37/110):  ref base (31/24):  total (68/134)
Fisher's exact test strand distribution p-value = 5.04e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765
                                                                                                  |                                                                                              
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                          >  2:2501847/1‑101 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                          >  1:4405163/1‑101 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                          <  1:280127/101‑1 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                          <  1:2680755/101‑1 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggg                                                                                        <  1:1614595/101‑1 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                       <  2:759779/101‑2 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                       >  1:1263499/1‑100 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                       >  2:3163275/1‑100 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggggg                                                                                     >  2:706320/1‑98 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGCCGAGTTCTGATCTC.AAGGC.C..gtgggg                                                                                     >  1:2016939/1‑96 (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:2509657/101‑4 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgct                                                                                 >  2:2147763/1‑98 (MQ=255)
         gcCGATCTGGTTTCGGAGACGGGTGAGGTGAAGCTCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgctg                                                                                <  1:4533785/101‑5 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgctg                                                                                >  2:2540087/1‑97 (MQ=255)
          ccGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggggggcggc                                                                                >  2:262658/1‑93 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggggggggcgcc                                                                               >  2:2281577/1‑92 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGg.ggcgctggc                                                                              >  2:2399650/1‑95 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.ccc                                                                              >  1:1313477/1‑101 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.cccc                                                                             <  2:1083610/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.CCCCCg                                                                           <  2:1210631/101‑1 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCgg                                                                          >  2:4489670/1‑101 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.CCCCCgg                                                                          <  2:1454236/101‑1 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.CCCCCGGa                                                                         <  1:265331/101‑1 (MQ=255)
                  gTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.CCCCCGGAc                                                                        <  2:4648895/101‑1 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.CCCCCGGAccc                                                                      <  2:4859428/101‑1 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.CCCCCGGAccc                                                                      <  2:1855499/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGTCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGAccccc                                                                    <  1:5462417/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGAccccc                                                                    <  1:3169292/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGAccccc                                                                    <  2:2305786/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGAccccc                                                                    <  1:369951/101‑1 (MQ=38)
                          gaTGGGTCAGGTGAAGCGCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGa                                                                <  2:3274262/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGat                                                               <  2:2718975/101‑1 (MQ=38)
                            agggTCAGGTGAAGGCCCCTTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGata                                                              <  2:489689/100‑1 (MQ=37)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGatat                                                             <  2:2796907/101‑1 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCGG.CCCCCGGGCCCCCCGGgtat                                                             >  1:122753/1‑101 (MQ=255)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATAttt                                                           <  1:5205265/101‑1 (MQ=38)
                                tCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTtg                                                          <  1:1719988/101‑1 (MQ=37)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTtgt                                                         <  2:1656755/101‑1 (MQ=37)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTtgt                                                         <  1:1069589/101‑1 (MQ=37)
                                  aGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.GGCTGCCCCCCGGACCCCCGGGaatttttg                                                         >  2:4547295/1‑94 (MQ=255)
                                  aGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTtgtg                                                        <  1:2124917/101‑1 (MQ=37)
                                   gatGAAGCACCCGTTCCGGGCTGTCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTgg                                                       <  2:5705357/99‑1 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCCAGa                                                  >  1:3693607/1‑101 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGa                                                  >  2:833570/1‑101 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGa                                                  <  1:362298/101‑1 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.GGCTG.CCCCCCGGCCCCCCGGGGATTTTTGGGCAGAt                                                 >  1:477968/1‑101 (MQ=255)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGAt                                                 >  2:3362457/1‑101 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGAt                                                 <  1:2400508/101‑1 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                                <  1:295320/101‑1 (MQ=25)
                                           ccccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                               <  2:1735120/100‑1 (MQ=25)
                                            cccGTTCCGTGCAGGCATCAAGGGCCTGCCGGGGGTCGAGTTCTGATCGCAAGGCC.G..GGGCCTCTG.CCCCCGGACCCCGGGGATATTTGTGGCAAGATGaa                                              <  2:3946524/101‑1 (MQ=21)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGaa                                              <  1:4558394/101‑1 (MQ=25)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGaa                                              <  2:399647/101‑1 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                            <  1:2208365/101‑1 (MQ=25)
                                               gtcccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACa                                           <  1:1473111/98‑1 (MQ=21)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                          <  2:1313389/101‑1 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                          <  1:1809756/101‑1 (MQ=25)
                                                tgccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                          <  2:406797/99‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                         <  1:2786316/101‑1 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                       >  2:5176689/1‑101 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                       <  2:1012128/101‑1 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                       >  1:693723/1‑101 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                       <  1:3143905/101‑1 (MQ=25)
                                                     ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGt                                     <  2:1547076/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                    >  2:585836/1‑101 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                    <  1:2281703/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                   <  2:2681015/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                   <  2:546848/101‑1 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.g                                 <  2:4449519/101‑1 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCCGGGAATTTGTGGCAAGATGAAACAGAGCGTTA.ggg                               >  1:5233832/1‑101 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCCG.CCCCCGGGCCCCCCGGATATTTGTGGCCAGATGAAACAGAGCGTTA.GGGAt                             >  1:3572145/1‑101 (MQ=37)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGAttt                           >  2:3729654/1‑101 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTTGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTg                          <  2:5198397/101‑1 (MQ=21)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.GGCGG.CCCCCGGACCCCCGGGAGATTTGTGGGGAGATGAAACAGAGCGCTA.GGGGTTTg                          >  1:3606128/1‑101 (MQ=21)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTg                          <  2:4337325/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTa                       <  1:2936203/101‑1 (MQ=25)
                                                                   gCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTaa                      <  1:5284617/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG.GGG.CGTTG.CCCCCGGACCCCCGGGATATTTGGGGCAGGATGAAACAGAGCGTTA.GGGATTTGGTAAc                     >  2:4761520/1‑101 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCa                   <  2:4541633/101‑1 (MQ=25)
                                                                       gCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCAt                  <  1:3992245/101‑1 (MQ=25)
                                                                         cGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGa                <  2:4607107/101‑1 (MQ=25)
                                                                         cGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGG.GGCCG.CCCCCCGGCCCCCCGGGTATTTGTGGGAAGATGAAAAAAAGGGGTAGGGGGTTTGGTACcctga                 >  2:4091135/1‑98 (MQ=255)
                                                                          gggggTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGAt               >  1:3503289/1‑101 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTC.AAGGC.CGGGGG.GGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATc              >  1:3354569/1‑101 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTC.AAGGC.CGGGGG.GGCTG.CCCCCCGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGGTTTGGTAACCATGATc              >  1:3780871/1‑101 (MQ=21)
                                                                           ggggTCGACTTCTGATCTCAAGGCCGG..GGG.CGCGG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATc              <  1:4996253/101‑1 (MQ=21)
                                                                            gggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCg             <  1:5339345/101‑1 (MQ=25)
                                                                            gggTCGAGTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCg             <  2:3066006/101‑1 (MQ=25)
                                                                             ggTCGAGTTCTGATCTC.AAGGC.CGGGGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATcag            >  2:3218442/1‑99 (MQ=25)
                                                                                  aGTTCTGATCTC.AAGGC.CGGGGG.GGCTG.CCCCCGGACCCCCCGGGAATTTTTGGGAAAATGAAAAAGAGCGTTA.GGGGTTTTGTAAACATGATTggggcaa       >  2:5151356/1‑97 (MQ=12)
                                                                                   gTTCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGc      <  1:3866575/101‑1 (MQ=25)
                                                                                    ttCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCt     <  1:833638/101‑1 (MQ=25)
                                                                                     tCTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTg    <  1:3451907/101‑1 (MQ=25)
                                                                                      cTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTgg   >  2:3588065/1‑101 (MQ=25)
                                                                                      cTGATCTCAAGGCCGG..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTgg   <  2:638526/101‑1 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG.GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGc  <  1:2540229/101‑1 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG.GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGc  <  2:4174532/101‑1 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG.GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGc  <  2:3693379/101‑1 (MQ=25)
                                                                                                  |                                                                                              
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGG.CGCTG.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765

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