Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,380 A→G 3.0% intergenic (‑176/‑129) parE ← / → RCAP_rcc02080 DNA topoisomerase IV subunit B/HAD superfamily hydrolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,249,3800A→G 3.0% 5.6 0.0 0.0236intergenic (‑176/‑129)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (5/3):  ref base (131/96):  total (137/99)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.CCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATC  >  NC_014034/2249281‑2249479
                                                                                                             |                                                                                                      
ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                       >  1:817140/1‑101 (MQ=255)
    gtcccgtGCCGGGACT.CCCCCGGGTTATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                                   >  2:3670845/1‑101 (MQ=255)
      cccgtGCCGGGACT.CCCCGGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                 <  2:107116/101‑1 (MQ=255)
      cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                 <  1:4980737/101‑1 (MQ=255)
       ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                <  2:608220/101‑1 (MQ=255)
       ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                >  2:3711690/1‑101 (MQ=255)
       ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                <  2:3319237/101‑1 (MQ=255)
       ccgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAAC.CGGA..TGCATTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggg                                                                                                    >  1:4872835/9‑101 (MQ=21)
        cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                               <  1:3498506/101‑1 (MQ=255)
        cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                               <  1:3904232/101‑1 (MQ=255)
        cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                               >  1:396248/1‑101 (MQ=255)
        cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                                   >  2:2381758/8‑101 (MQ=18)
          tGCCTGGACT.CCCCCGGGATATTTGAGGCAAGATCAG...GGGG..GGT.TTTCATCTTTCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAgcg                                                                                             <  1:2240441/101‑1 (MQ=37)
          caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                                                 <  2:902873/96‑1 (MQ=18)
           gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAgcg.c                                                                                           >  2:2661997/1‑101 (MQ=255)
           gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAgag.c                                                                                           >  2:1737242/1‑101 (MQ=255)
             gcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAgc                                                                                              <  1:2286656/99‑1 (MQ=18)
             cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc                                                                                         <  1:2591059/101‑1 (MQ=255)
             cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc                                                                                         >  1:1052733/1‑101 (MQ=255)
              gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc.c                                                                                       >  1:4451631/1‑101 (MQ=255)
              cggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAgcg                                                                                             <  1:276066/100‑1 (MQ=18)
                gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAGCG.cc                                                                                          <  1:4232038/96‑1 (MQ=17)
                 aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc.ccct                                                                                    >  1:2140057/1‑99 (MQ=255)
                 aCT.CCCCAGGGTTATTTGAGGCAAGATGTA...GGGG..GGC.TTTCATCTTGGCGAAAATTTGCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc.ccct                                                                                    <  2:4101922/101‑3 (MQ=25)
                   t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc.cccttc                                                                                  <  2:1645830/101‑5 (MQ=255)
                   t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc.cccttc                                                                                  <  1:2181228/101‑5 (MQ=255)
                     cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc.cccttcg                                                                                 <  1:2330065/101‑6 (MQ=255)
                      ccccGGGATATTTGCGGCAAGGTGAA...GGGGCATGA.TTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAGCG.CCC.CC.tg                                                                                   >  1:2863925/1‑100 (MQ=2)
                      ccccGGGATATTTGCGGCAAGATGAA...GGGGCATGA.TTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAGCG.CCC.CC.tg                                                                                   >  1:3435718/1‑100 (MQ=0)
                      ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTcgc                                                                                <  2:4298079/101‑1 (MQ=255)
                       cccGGGATATTTGCGGCAAGATGAA...GGGGCATGA.TTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAGCG.CCC.CC.tgc                                                                                  <  2:1204802/101‑3 (MQ=0)
                         cGGGATATTTGCGGCAAGATGAA...GGGGCATGA.TTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAGCG.CCC.CC.tgccc                                                                                >  2:3195667/1‑97 (MQ=2)
                          gggATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCgcgggc                                                                            <  1:654632/101‑1 (MQ=255)
                            gATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCgcgggcgg                                                                          >  1:4569746/1‑101 (MQ=255)
                             atatTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGggcggc                                                                         >  1:74012/1‑101 (MQ=255)
                              tatTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGcc                                                                        <  2:2757136/101‑1 (MQ=255)
                               atTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCcg                                                                       >  2:1375132/1‑101 (MQ=255)
                                tttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCcgc                                                                      <  2:4880746/101‑1 (MQ=255)
                                tttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCcgc                                                                      >  1:4305555/1‑101 (MQ=255)
                                 ttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCcgcg                                                                     >  1:4922108/1‑101 (MQ=255)
                                  tGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGa                                                                    <  1:3864265/101‑1 (MQ=255)
                                  tGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGGGGGCGGCCGCGa                                                                    >  1:477324/1‑101 (MQ=37)
                                   gAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGAc                                                                   >  2:4379902/1‑101 (MQ=255)
                                   gAGGAAAGATGAA...GGGG..GGC.TTTCATTTTGCCGAAAATTTCGGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGAc                                                                   <  2:1789784/101‑1 (MQ=37)
                                    aggtaggaggAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGAcg                                                                  <  2:1582857/92‑1 (MQ=25)
                                     ggCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGAcgc                                                                 <  2:3458686/101‑1 (MQ=255)
                                      gCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGAcgcg                                                                <  2:4658103/101‑1 (MQ=255)
                                       cAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGAcgcgc                                                               <  2:341805/101‑1 (MQ=255)
                                       cAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.CCCACC.TTCGCGGGCGGCCGCGAcgcgc                                                               >  2:1280495/1‑101 (MQ=255)
                                          gATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.CCC.CCATTCGCGGGCGGCCGCGACGCGcgcc                                                            >  2:225834/1‑101 (MQ=255)
                                           aTGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGcgccg                                                           >  1:2085237/1‑101 (MQ=255)
                                             gAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                         >  2:515019/1‑101 (MQ=255)
                                             gAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                         <  2:540454/101‑1 (MQ=255)
                                              aa...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCGGCaa                                                        >  1:2565231/1‑101 (MQ=255)
                                              aa...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                        >  2:249070/1‑101 (MQ=255)
                                              aa...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                        >  1:1143316/1‑101 (MQ=255)
                                               a...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                       >  2:1840746/1‑101 (MQ=255)
                                               a...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                       <  2:2403842/101‑1 (MQ=255)
                                               a...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                       <  1:2855191/101‑1 (MQ=255)
                                               a...GGGG..GGC.TTTCATCTTGCCAAAAATGTCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                       <  2:396183/101‑1 (MQ=37)
                                                   gggg..ggC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                      <  2:2897380/101‑1 (MQ=255)
                                                          gC.TTTCATCATGCCGAAAATATCCCGGGGGGTGAA.TGGCCCA..GCCAGGGGGGGGAAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGCCg                                                 >  2:1677257/1‑101 (MQ=37)
                                                           c.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                                <  2:3826913/101‑1 (MQ=255)
                                                             tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                               <  1:3109615/101‑1 (MQ=255)
                                                             tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                               >  1:2365250/1‑101 (MQ=255)
                                                             tttCATCTTGCCGAAAATATCCCGGGGGCTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                               >  2:398831/1‑101 (MQ=255)
                                                               tCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                             >  1:3201183/1‑101 (MQ=255)
                                                               tCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                             >  2:1030871/1‑101 (MQ=255)
                                                               tCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                             >  1:2075832/1‑101 (MQ=255)
                                                                cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                            >  1:2946206/1‑101 (MQ=255)
                                                                  tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                          >  2:2699450/1‑101 (MQ=255)
                                                                   cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                         >  1:4349614/1‑101 (MQ=255)
                                                                    ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                        <  2:4573046/101‑1 (MQ=255)
                                                                       ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                     <  1:4510971/101‑1 (MQ=255)
                                                                       ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                     <  2:4216837/101‑1 (MQ=255)
                                                                          aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                  >  2:5019650/1‑101 (MQ=255)
                                                                           aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                 >  1:4295762/1‑101 (MQ=255)
                                                                              tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                              >  2:3526695/1‑101 (MQ=255)
                                                                               acccgggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                             <  1:1874073/96‑1 (MQ=255)
                                                                                   gggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                         <  2:1009809/100‑1 (MQ=255)
                                                                                    ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                        <  2:2439222/101‑1 (MQ=255)
                                                                                     gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                       >  2:1458945/1‑101 (MQ=255)
                                                                                     gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                       >  1:702558/1‑101 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTGGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                     >  1:1550806/1‑101 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                     >  2:2503832/1‑101 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                     >  1:4000052/1‑101 (MQ=255)
                                                                                        ggTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TGCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                    <  1:4829377/101‑1 (MQ=255)
                                                                                           gAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                 >  2:2122268/1‑101 (MQ=255)
                                                                                           gAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                 >  1:3730874/1‑101 (MQ=255)
                                                                                           gAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGACCCTTTGTCATTCTGTCGCACCCATGTCa                 >  2:824026/1‑101 (MQ=255)
                                                                                              acgccgcA..GCCAGAGGGGGGCAGCGCCCC.CC.TCCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTCCCCCTTTGTCATTCTGTCGCACCCATGTCATa               <  2:4719794/97‑1 (MQ=255)
                                                                                               tGGCGCA..GCGGGGGGGGGGCGGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa              <  1:579683/101‑1 (MQ=255)
                                                                                               tGGCGCA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa              >  2:2808364/1‑101 (MQ=255)
                                                                                                    cA..GCCAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg         <  2:3100241/101‑1 (MQ=255)
                                                                                                         ccAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg      >  2:4386321/1‑101 (MQ=255)
                                                                                                         ccAGAGGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg      >  2:1941327/1‑101 (MQ=255)
                                                                                                           agagGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa    >  2:3004394/1‑101 (MQ=255)
                                                                                                            gagGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt   >  1:4947990/1‑101 (MQ=255)
                                                                                                            gagGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt   >  1:3453822/1‑101 (MQ=255)
                                                                                                            gagGGGGGCAGCGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt   >  1:1576532/1‑101 (MQ=255)
                                                                                                             |                                                                                                      
CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.CCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATC  >  NC_014034/2249281‑2249479

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG