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% 156.8 168intergenic (+10/‑261)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (32/91):  ref base (23/22):  total (55/113)
Fisher's exact test strand distribution p-value = 2.96e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.69e-01
*NC_0140342,209,6790C→G100.0% 115.9 168intergenic (+11/‑260)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (32/93):  ref base (23/20):  total (55/113)
Fisher's exact test strand distribution p-value = 1.24e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.81e-01

Only 100 of 195 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780
                                                                                                     ||||                                                                                                      
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                                        >  2:983916/1‑101 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                                        <  1:1109642/101‑1 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                                        <  1:3196282/101‑1 (MQ=255)
 tCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gg                                                                                                       <  2:2040983/101‑1 (MQ=255)
 tCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gg                                                                                                       <  1:2515611/101‑1 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.C..gg                                                                                                       >  1:179292/1‑101 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTCCTGATCTC.AAGGC.C..ggg                                                                                                      >  2:3803056/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                                     >  2:2825330/1‑100 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                                     <  2:3033731/101‑2 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggc                                                                                                   <  1:366962/101‑1 (MQ=255)
     ggcggcCGATCAAGTGACCGAGATGGGTCANTTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggc                                                                                                   <  1:662836/101‑1 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggggggct                                                                                                >  2:1859007/1‑95 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggggggggg                                                                                               >  2:2994868/1‑94 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgctgcc                                                                                             <  1:1590883/101‑7 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgctgcc                                                                                             <  2:2242895/101‑7 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGc.cc                                                                                             >  2:3423477/1‑101 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgcgcgcc                                                                                            >  2:39902/1‑94 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGc.cccc                                                                                           >  1:1405025/1‑101 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGc.cccc                                                                                           >  2:1801717/1‑101 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGc.cccc                                                                                           <  1:1058664/101‑1 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGc.cccc                                                                                           <  2:1705952/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCg                                                                                          >  2:2533287/1‑101 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCg                                                                                          >  1:2379256/1‑101 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCg                                                                                          <  2:2465617/101‑1 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCgg                                                                                         <  1:1350096/101‑1 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGAcc                                                                                      >  1:1297425/1‑101 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGAcc                                                                                      <  2:4208403/101‑1 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGAccc                                                                                     <  1:2792734/101‑1 (MQ=255)
                     aCCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGAcccc                                                                                    <  1:206580/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGAccccc                                                                                   <  1:3552908/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGAccccc                                                                                   <  2:3831136/101‑1 (MQ=38)
                        gagaTGGGTCAGGTGAAGCGCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCgg                                                                                 <  1:4086549/101‑1 (MQ=38)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGC.CCCCCGACCCCCgg                                                                                 >  1:3321387/1‑101 (MQ=255)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCgg                                                                                 >  2:3362928/1‑101 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCggg                                                                                >  2:606764/1‑101 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCggg                                                                                >  1:2673428/1‑101 (MQ=38)
                          gaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGa                                                                               <  2:2687348/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGat                                                                              <  1:2600575/101‑1 (MQ=38)
                            tGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGata                                                                             >  2:433168/1‑101 (MQ=38)
                            tGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGata                                                                             <  1:2995090/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATAtt                                                                           >  2:412672/1‑101 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATAtt                                                                           <  2:2409563/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATAtt                                                                           >  1:1922993/1‑101 (MQ=38)
                                tCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTtg                                                                         <  2:2700743/101‑1 (MQ=37)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTtgt                                                                        <  2:2859474/101‑1 (MQ=37)
                                   ggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTgg                                                                      <  2:1081468/101‑1 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGc                                                                     <  1:1874437/101‑1 (MQ=37)
                                     tGAAGCCCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCa                                                                    <  2:4104028/101‑1 (MQ=25)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCa                                                                    <  1:216122/101‑1 (MQ=37)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCa                                                                    <  1:2633094/101‑1 (MQ=37)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCa                                                                    <  1:3803313/101‑1 (MQ=37)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAg                                                                  <  2:74983/101‑1 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAg                                                                  <  2:259871/101‑1 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAg                                                                  <  2:571198/101‑1 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGa                                                                 <  1:852365/101‑1 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGa                                                                 <  2:373735/101‑1 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGAt                                                                <  1:3454090/101‑1 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCACGATGa                                                              >  2:531327/1‑101 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                                              >  2:1918427/1‑101 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                                           <  2:1277562/101‑1 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGACGG..GGGGGGTGCACCCCGGCCCCCCGGATAATTTGTGGCaaaagaaaaa                                                           >  1:4092268/1‑93 (MQ=255)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                        <  1:2886273/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                                        <  1:750155/101‑1 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                                      <  2:3556104/101‑1 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CGGGGGGGCGGC.CCCCGGCCCCCCGGGGAATTTTTGGCAAAAGGaaaagaag                                                       >  1:3785806/1‑96 (MQ=255)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGC.CCCCGGGCCCCCCGGATATTTGTGGGAAGATGAAACagaggg                                                     >  1:3469470/1‑99 (MQ=25)
                                                     ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.Gt                                                   >  2:1938589/1‑101 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTa                                                 <  2:114596/101‑1 (MQ=25)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAgg                                               <  1:1164002/101‑1 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGC.CCCCGGACCCCCCGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGAt                                            >  2:1001901/1‑101 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCCGACCCCCGGGATATTTGTGGCAAGATGAAAAAGAGC.CGTAggggt                                            >  2:1149381/1‑99 (MQ=18)
                                                                aGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTgg                                        <  2:883780/101‑1 (MQ=25)
                                                                 nnnnnnnnnnggggTCGNGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGt                                       <  2:617117/91‑1 (MQ=18)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTa                                      <  1:169210/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGGGCCGC.CCCCCGGCCCCCCGGGGATTTGTGGCCAGAAGAAACAGAGC.GGTAGGGGttttgga                                      >  1:413487/1‑97 (MQ=255)
                                                                   gCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTaa                                     <  2:1959642/101‑1 (MQ=25)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAAcc                                   >  2:3168301/1‑101 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTC.AAGGC.CGGGGGGGCTGC.CCCCGGGCCCCCCGGGTATTTGTGGGAAGATGAAACAGAGC.GTTAGGGGTTTGGGAACCATGAt                              >  1:474814/1‑101 (MQ=21)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATc                             <  2:830607/101‑1 (MQ=25)
                                                                            gggTCGAGTTCTGATCTC.AAGGC.CGGGGGGGGCGC.CCCCCGGCCCCCCGGGTATTTTTGGGAAGATGAAAAAGAGCGGTTGGGGATTTGG.AAACATGAttg                            >  2:1770582/1‑99 (MQ=255)
                                                                             ggTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCgg                           <  2:3728632/101‑1 (MQ=25)
                                                                              gTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCggg                          <  1:1039616/101‑1 (MQ=25)
                                                                               tCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGc                         <  2:3709666/101‑1 (MQ=25)
                                                                               tCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGc                         <  1:3043531/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCa                        <  1:606828/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCa                        >  2:4170324/1‑101 (MQ=25)
                                                                                 gAGTTCTGATCTC.AAGGC.CGGGGGGGCCGC.CCCCCGGCCCCCCGGGTATTTTTGGGAAGATGAAAAAGAGG.GGTAGGGGGTTTGGAAACATGATCGGGGaa                       >  1:3717245/1‑101 (MQ=255)
                                                                                  aGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAg                      >  2:2844692/1‑101 (MQ=25)
                                                                                    ttCTGATCTCGATGCC.GG.GGGCGCTGC.CCCAGGACACCAGGGCTATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATTGGGCAAGCt                    <  2:1669894/101‑1 (MQ=21)
                                                                                    ttCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCt                    <  2:1308789/101‑1 (MQ=25)
                                                                                              tcaaggccgggggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAacac            <  1:2756402/91‑1 (MQ=255)
                                                                                               caaggccgggggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACg           >  2:739926/10‑101 (MQ=255)
                                                                                                 aggccgggggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGt         <  1:1338851/94‑1 (MQ=255)
                                                                                                 aggccgggggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGt         <  1:433225/94‑1 (MQ=255)
                                                                                                 aggccgggggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGt         <  2:1555484/94‑1 (MQ=255)
                                                                                                    ccgggggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACa      <  1:887848/97‑1 (MQ=255)
                                                                                                      gggggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  1:2604407/3‑101 (MQ=255)
                                                                                                       ggggCGCTGC.CCCCGAACCCCAAGGATATTTGTGGCAAGATGAAGCAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACTCGGTACaggg   <  2:3928770/100‑1 (MQ=38)
                                                                                                        gggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaggga  >  1:2697328/1‑101 (MQ=255)
                                                                                                        gggCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaggga  <  2:3047919/101‑1 (MQ=255)
                                                                                                     ||||                                                                                                      
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGC.GTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780

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