Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 408,497 +GCAGA 100% intergenic (+115/+33) glnA1 → / ← dapA glutamine synthetase/dihydrodipicolinate synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034408,4971.→G100.0% 55.6 137intergenic (+115/+33)glnA1/dapAglutamine synthetase/dihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  new base (80/36):  ref base (5/13):  total (86/51)
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.00e+00
*NC_014034408,4972.→C100.0% 58.7 137intergenic (+115/+33)glnA1/dapAglutamine synthetase/dihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  new base (80/36):  ref base (5/14):  total (86/51)
Fisher's exact test strand distribution p-value = 6.15e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.60e-01
*NC_014034408,4973.→A100.0% 57.7 137intergenic (+115/+33)glnA1/dapAglutamine synthetase/dihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  new base (81/36):  ref base (5/14):  total (86/51)
Fisher's exact test strand distribution p-value = 5.73e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01
*NC_014034408,4974.→G100.0% 59.0 137intergenic (+115/+33)glnA1/dapAglutamine synthetase/dihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  new base (81/36):  ref base (5/14):  total (86/51)
Fisher's exact test strand distribution p-value = 5.73e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01
*NC_014034408,4975.→A100.0% 61.7 137intergenic (+115/+33)glnA1/dapAglutamine synthetase/dihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  new base (81/36):  ref base (5/15):  total (86/51)
Fisher's exact test strand distribution p-value = 2.82e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.64e-01

Only 100 of 145 total aligned reads displayed.

CTGCCCCCAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGCAA.AGACGGGG..G.....GC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAGCGC  >  NC_014034/408396‑408568
                                                                                                          |||||                                                                            
ctGCCCCCAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACgggg..g                                                                                   >  2:3364425/1‑101 (MQ=255)
  gCCCCCAGACCCCGGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gc                                                                            <  1:3096902/101‑1 (MQ=255)
   tcgcctcaCCCCCGGGATATTTGGGGCAAGATAAACCCCAT.CCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCA.AGAGgggg..gca                                                                                 >  1:3326953/8‑99 (MQ=25)
   cccccAGCCCCCCGGGATATTTGGGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGACTT.GCTGAAAGC.A.AGACGGGG..G.....Gca                                                                           <  1:4025070/101‑2 (MQ=25)
     cccAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcaga                                                                         <  2:1757956/101‑4 (MQ=255)
     cccAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcaga                                                                         <  1:1377659/101‑4 (MQ=255)
     cccAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcaga                                                                         <  2:3358852/101‑4 (MQ=255)
     cccAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcaga                                                                         <  1:2468328/101‑4 (MQ=255)
     cccAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcaga                                                                         <  1:3679316/101‑4 (MQ=255)
      ctaGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagag                                                                        <  1:3262211/99‑5 (MQ=255)
      ccAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagag                                                                        <  2:878839/101‑5 (MQ=255)
      ccAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagag                                                                        >  2:2967539/1‑97 (MQ=255)
        agcccccgGGGATTTTTGCGGAAAGAGAAAGCCGATGTCCTTAACTTTGCCGAAAATTTCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagagcc                                                                      <  1:3196961/98‑7 (MQ=255)
        aGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagagcc                                                                      >  2:4153993/1‑95 (MQ=255)
        aGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAG.AA.AGACGGGG..G.....Gcagagag                                                                      >  2:3261970/1‑95 (MQ=255)
           ccccgGGTATTTTGGCGGAAAAATGAACCCGTTTCCTTTCATCTTCCCAAAAATCCCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagagccccc                                                                   <  1:879360/101‑10 (MQ=255)
           cccccGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagagccccc                                                                   >  1:608292/1‑92 (MQ=25)
            ccccGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagagcccccc                                                                  <  1:2092163/101‑11 (MQ=25)
            ccccGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....Gcagagcccccc                                                                  >  2:1859944/1‑91 (MQ=25)
                          cGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaa                                                    >  2:1981973/1‑101 (MQ=21)
                          cGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaa                                                    >  1:586166/1‑101 (MQ=21)
                          cGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaa                                                    <  2:962094/101‑1 (MQ=21)
                           ggCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaaa                                                   >  2:1054739/1‑101 (MQ=21)
                           ggCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaaa                                                   >  1:2603345/1‑101 (MQ=21)
                            gaaaGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaaaa                                                  <  2:3597311/99‑1 (MQ=21)
                            gCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaaaa                                                  <  2:3675828/101‑1 (MQ=21)
                            gCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaaaa                                                  <  2:2454172/101‑1 (MQ=21)
                            gCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaaaa                                                  >  2:1244510/1‑101 (MQ=21)
                            gCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaaaa                                                  <  1:2277220/101‑1 (MQ=21)
                            gCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCaaaa                                                  >  2:3266524/1‑101 (MQ=21)
                             cAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAg                                                 <  2:3143514/101‑1 (MQ=21)
                             cAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAg                                                 >  2:2452975/1‑101 (MQ=21)
                              aaGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAgg                                                >  1:342787/1‑101 (MQ=21)
                              aaGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAgg                                                >  1:2762653/1‑101 (MQ=21)
                               aGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACGGGTCGGCAAAAGga                                               >  1:3656808/1‑101 (MQ=18)
                               aGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGCCGCCAAAAGga                                               >  2:460123/1‑101 (MQ=18)
                                gATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGGAGAGC.....CCCCCGCACAGGGCGCCAAAAGgac                                              >  2:1508479/1‑101 (MQ=21)
                                gATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGgac                                              <  2:206684/101‑1 (MQ=21)
                                gATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGgac                                              <  1:1949175/101‑1 (MQ=21)
                                gATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGgac                                              >  1:2759788/1‑101 (MQ=21)
                                 aTGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGgacc                                             >  1:17272/1‑101 (MQ=21)
                                   gAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGgaccga                                           >  2:2816777/1‑101 (MQ=21)
                                   gAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGgaccga                                           >  1:2783109/1‑101 (MQ=21)
                                     aGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGgaccgacc                                         >  1:1031664/1‑101 (MQ=21)
                                      gCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCt                                        >  2:2256487/1‑101 (MQ=21)
                                      gCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCt                                        >  1:1513656/1‑101 (MQ=21)
                                      gCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GACAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCt                                        >  2:2825734/1‑101 (MQ=21)
                                      gCCGATGTCCTTCATCTTGCCGCAAATACCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCt                                        <  2:1816251/101‑1 (MQ=21)
                                       ccGATGTCCTTCATTTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGCAA.GAAGGGGG..GCAGA.GC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTc                                       <  1:4182524/101‑1 (MQ=21)
                                        cGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCa                                      >  1:107444/1‑101 (MQ=21)
                                        cGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCa                                      <  2:608228/101‑1 (MQ=21)
                                        cGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCa                                      >  1:2055735/1‑101 (MQ=21)
                                         gATGTCCTTCATCTTGCGGCAAATTCCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAg                                     <  1:30138/101‑1 (MQ=255)
                                         gATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAg                                     <  2:2747460/101‑1 (MQ=21)
                                         gATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAg                                     <  2:3830015/101‑1 (MQ=21)
                                         gATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAg                                     >  1:1212333/1‑101 (MQ=21)
                                          aTGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGa                                    >  2:409701/1‑101 (MQ=21)
                                          aTGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGa                                    >  2:2332950/1‑101 (MQ=21)
                                          aTGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGa                                    <  1:2457365/101‑1 (MQ=21)
                                           tGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGAt                                   <  1:976197/101‑1 (MQ=21)
                                           tGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGAt                                   >  1:1249390/1‑101 (MQ=21)
                                           tGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGAt                                   >  2:2868143/1‑101 (MQ=21)
                                           ggcccTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGAt                                   <  2:2485896/98‑1 (MQ=18)
                                            gcccTTCATCTGGACGCAAATACCACGGGGGGAGGCTT.GCTGAAAGCAAGAGGGGGGGCAGA....GC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATc                                  <  2:1631936/99‑1 (MQ=255)
                                             tCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCa                                 <  1:432479/101‑1 (MQ=21)
                                             tCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCa                                 <  1:1169721/101‑1 (MQ=21)
                                              ccTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCTACCTCAGATCAg                                >  2:3778110/1‑101 (MQ=18)
                                              ccTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAg                                >  1:1998602/1‑101 (MQ=21)
                                              ccTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAg                                >  1:622243/1‑101 (MQ=21)
                                                ttCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGcc                              >  2:1434319/1‑101 (MQ=21)
                                                  cATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCg                            >  2:1424036/1‑101 (MQ=21)
                                                  cATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCg                            >  2:2769081/1‑101 (MQ=21)
                                                   aTCTTGCGGCAAATATCTCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCgg                           <  2:2608844/101‑1 (MQ=255)
                                                   aTCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAAccccg                           >  1:782389/1‑99 (MQ=21)
                                                    tCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGc                          >  1:3785937/1‑101 (MQ=21)
                                                    tCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGc                          <  2:1365626/101‑1 (MQ=21)
                                                       tGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..G.....GCAGAGGCCCCCCCACAGGTCGCCAAAAGGGCCGACCTCAGATCAGCCCGGCgtg                       >  1:14229/1‑101 (MQ=21)
                                                        gCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCgtgt                      <  1:223354/101‑1 (MQ=25)
                                                           gCAAATATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGc                   >  2:3577842/1‑101 (MQ=25)
                                                              aaTATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATc                >  1:701751/1‑101 (MQ=255)
                                                              aaTATCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATc                >  2:1930649/1‑101 (MQ=255)
                                                               atatCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATcg               >  2:1627531/1‑101 (MQ=255)
                                                               atatCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATcg               >  2:3954020/1‑101 (MQ=255)
                                                                 atCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATcgcg             >  1:2619662/1‑101 (MQ=255)
                                                                  tCCCGGGGGGAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCgcgg            <  1:472022/101‑1 (MQ=255)
                                                                      ggggggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTc        <  1:3803522/101‑1 (MQ=255)
                                                                      ggggggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTc        >  1:4003930/1‑101 (MQ=255)
                                                                      ggggggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGAACTCAGATCAGCCCGGCGTGTTGCATCGCGGCGta        >  2:3023796/1‑100 (MQ=255)
                                                                       gggggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGGGTCa       >  2:174247/1‑101 (MQ=255)
                                                                       gggggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCa       >  2:256970/1‑101 (MQ=255)
                                                                       gggggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCa       >  2:2983468/1‑101 (MQ=255)
                                                                       gggggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCa       >  2:193211/1‑101 (MQ=255)
                                                                        ggggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCaa      <  1:4154255/101‑1 (MQ=255)
                                                                         gggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAg     >  2:3722250/1‑101 (MQ=255)
                                                                         gggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAg     >  2:56671/1‑101 (MQ=255)
                                                                          ggAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAgc    >  1:380876/1‑101 (MQ=255)
                                                                           gAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAgcg   >  2:383005/1‑101 (MQ=255)
                                                                           gAGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAgcg   >  2:2020580/1‑101 (MQ=255)
                                                                            aGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAgcgc  >  1:1369964/1‑101 (MQ=255)
                                                                            aGGCTT.GCTGAAAGC.A.AGACGGGG..GGCAGAGC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAgcgc  >  1:893713/1‑101 (MQ=255)
                                                                                                          |||||                                                                            
CTGCCCCCAGACCCCCGGGATATTTGCGGCAAGATGAAGCCGATGTCCTTCATCTTGCCGCAAATATCCCGGGGGGAGGCTT.GCTGAAAGCAA.AGACGGGG..G.....GC.....CCCCCGCACAGGTCGCCAAAAGGACCGACCTCAGATCAGCCCGGCGTGTTGCATCGCGGCGTCAAGCGC  >  NC_014034/408396‑408568

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG