Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,358 Δ1 bp 2.2% intergenic (‑154/‑151) 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,3580G→. 2.2% 6.8 ‑1.0 ‑0.0228intergenic (‑154/‑151)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (123/100):  total (127/101)
Fisher's exact test strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.00e-01

Only 100 of 252 total aligned reads displayed.

CGGGGCAAGAAGGGGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACC  >  NC_014034/2249260‑2249457
                                                                                                         |                                                                                                          
  gggCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggg                                                                                                         >  2:4510651/1‑101 (MQ=255)
    gCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                                                       <  1:2842509/101‑1 (MQ=255)
     cAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGt                                                                                                      <  2:4790983/101‑1 (MQ=255)
      aagaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTg                                                                                                     >  1:308160/1‑101 (MQ=255)
      aagaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTg                                                                                                     >  2:365312/1‑101 (MQ=255)
      aagaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTg                                                                                                     >  1:2483520/1‑101 (MQ=255)
      aagaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTg                                                                                                     >  2:4681391/1‑101 (MQ=255)
          agGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.Tg                                                                                                >  1:25854/1‑101 (MQ=255)
           gggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.Tgg                                                                                               <  2:3118500/101‑1 (MQ=255)
            ggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgc                                                                                              >  2:4840296/1‑101 (MQ=255)
            ggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgc                                                                                              <  2:2188273/101‑1 (MQ=255)
            ggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgc                                                                                              <  2:2296008/101‑1 (MQ=255)
             gggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcg                                                                                             <  1:255547/101‑1 (MQ=255)
             gggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcg                                                                                             <  2:4028024/101‑1 (MQ=255)
              ggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGgcgc                                                                                            >  1:1234036/1‑101 (MQ=255)
                ctctGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..g                                                                                        <  1:4295982/101‑1 (MQ=255)
                ctctGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..g                                                                                        <  1:1751733/101‑1 (MQ=255)
                  ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCNNNNNTATCCCGGGGGTTGAA.TCGCGCA..Gcc                                                                                      >  1:3250151/5‑101 (MQ=25)
                  ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                                                      <  2:4878389/97‑1 (MQ=255)
                  ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..Gcc                                                                                      >  1:3524313/5‑101 (MQ=255)
                  ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..Gcc                                                                                      >  1:3458951/5‑101 (MQ=255)
                    gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                    >  1:3177311/3‑101 (MQ=255)
                    gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                    >  1:3994163/3‑101 (MQ=255)
                    gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                    <  2:4133381/99‑1 (MQ=255)
                    gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                    <  1:1942502/99‑1 (MQ=255)
                         cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                               >  1:2776911/1‑101 (MQ=255)
                         cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAACATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                               >  2:579604/1‑101 (MQ=255)
                          gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGAGAA.TGGCGCA..GCCAGAggggg                                                                              >  2:1813618/1‑101 (MQ=255)
                           tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                             <  2:4118106/101‑1 (MQ=255)
                           atctgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                             <  1:2593148/97‑1 (MQ=255)
                            cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                            <  1:4980737/101‑1 (MQ=255)
                            cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                            >  1:2415653/1‑101 (MQ=255)
                            cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                            <  1:2111314/101‑1 (MQ=255)
                            cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGc                                                                            >  2:2443535/1‑101 (MQ=255)
                             ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                           <  2:608220/101‑1 (MQ=255)
                              cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                          >  1:396248/1‑101 (MQ=255)
                              cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                          <  1:3904232/101‑1 (MQ=255)
                              cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                              >  2:1685360/8‑101 (MQ=21)
                                caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                            >  2:1143867/6‑101 (MQ=18)
                                caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                            <  2:902873/96‑1 (MQ=18)
                                 gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                                                     >  2:2661997/1‑101 (MQ=255)
                                   gcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c                                                                        <  1:2286656/99‑1 (MQ=18)
                                   cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc                                                                   >  1:1052733/1‑101 (MQ=255)
                                    gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.c                                                                 >  1:4451631/1‑101 (MQ=255)
                                    gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.c                                                                 >  1:929201/1‑101 (MQ=255)
                                     ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.cg.c                                                                     <  2:3846146/101‑1 (MQ=17)
                                       ccT.CCCCAAGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.ccct                                                              <  2:3823789/100‑3 (MQ=255)
                                       aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.ccct                                                              >  1:485248/1‑99 (MQ=255)
                                        cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGAGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.ccctt                                                             >  2:1626094/1‑98 (MQ=255)
                                         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.cccttc                                                            <  1:2181228/101‑5 (MQ=255)
                                         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.cccttc                                                            <  2:1306701/101‑5 (MQ=255)
                                         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.cccttc                                                            >  1:4658426/1‑97 (MQ=255)
                                            ccccGGGATATTTGCGGCAAGGTGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCC.CC.tg                                                             >  1:2863925/1‑100 (MQ=2)
                                            ccccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCC.CC.tg                                                             >  1:3435718/1‑100 (MQ=0)
                                            ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTcgc                                                          <  2:4298079/101‑1 (MQ=255)
                                             cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCC.CC.tgc                                                            <  2:1204802/101‑3 (MQ=0)
                                                gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCgcgggc                                                      <  1:654632/101‑1 (MQ=255)
                                                  gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCC..GCCAGAGGGGGGGAG.CGCCCC.CC.TTCgccggcgc                                                    >  2:1464777/1‑100 (MQ=25)
                                                  gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCgcgggcgg                                                    >  1:4569746/1‑101 (MQ=255)
                                                   atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGAAG.CGCCCC.CC.TTCCGGggcggc                                                   >  1:1142867/1‑101 (MQ=25)
                                                   atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGggcggc                                                   >  1:2251467/1‑101 (MQ=255)
                                                   atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGggcggc                                                   >  1:74012/1‑101 (MQ=255)
                                                    tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGcc                                                  <  2:2757136/101‑1 (MQ=255)
                                                      tttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCcgc                                                >  1:4305555/1‑101 (MQ=255)
                                                       ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCcgcg                                               >  2:1926736/1‑101 (MQ=255)
                                                       ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCcgcg                                               >  2:508930/1‑101 (MQ=255)
                                                       ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCcgcg                                               >  1:4922108/1‑101 (MQ=255)
                                                       ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCcgcg                                               >  2:535086/1‑101 (MQ=255)
                                                        tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGa                                              <  1:3864265/101‑1 (MQ=255)
                                                        tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGGGGGCGGCCGCGa                                              >  1:477324/1‑101 (MQ=37)
                                                             cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCACC.TTCGCGGGCGGCCGCGAcgcgc                                         >  2:1280495/1‑101 (MQ=255)
                                                                gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCC.CCATTCGCGGGCGGCCGCGACGCGcgcc                                      >  2:225834/1‑101 (MQ=255)
                                                                   gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                   >  2:2560414/1‑101 (MQ=255)
                                                                    aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                  >  2:1052694/1‑101 (MQ=255)
                                                                    aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                  >  1:2076968/1‑101 (MQ=255)
                                                                    aa...GGGC..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                  >  2:3693541/1‑101 (MQ=255)
                                                                     a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                 <  2:2403842/101‑1 (MQ=255)
                                                                     a...GGGG..GGCTTTCATCTTGCCAAAAATGTCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                 <  2:396183/101‑1 (MQ=37)
                                                                         gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                <  2:349745/101‑1 (MQ=255)
                                                                         gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                <  2:2897380/101‑1 (MQ=255)
                                                                         gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                <  1:4133575/101‑1 (MQ=255)
                                                                                   ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                        <  2:4382807/101‑1 (MQ=255)
                                                                                     cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                      >  1:2946206/1‑101 (MQ=255)
                                                                                       tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                    >  2:2699450/1‑101 (MQ=255)
                                                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                   >  1:4349614/1‑101 (MQ=255)
                                                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                   <  2:3394470/101‑1 (MQ=255)
                                                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                   >  1:4263215/1‑101 (MQ=255)
                                                                                         ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                  >  2:1786414/1‑101 (MQ=255)
                                                                                         tggCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                  <  2:329169/99‑1 (MQ=255)
                                                                                             ggAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAT.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa              <  1:2340286/100‑1 (MQ=255)
                                                                                             ggAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa              <  2:1013479/100‑1 (MQ=255)
                                                                                             cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa              >  1:2150742/1‑101 (MQ=255)
                                                                                               aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt            >  2:5019650/1‑101 (MQ=255)
                                                                                               aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt            >  1:3784659/1‑101 (MQ=255)
                                                                                                 aaTTTCCGGGGGGGTGAA.TGGCGCA..GCCAAAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGCGCCCCTTTGTTATTCt          >  2:165061/1‑101 (MQ=37)
                                                                                                   tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt        >  2:3526695/1‑101 (MQ=255)
                                                                                                    atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc       >  2:2879455/1‑101 (MQ=255)
                                                                                                    atCCCGGGGGATGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc       >  2:432316/1‑101 (MQ=255)
                                                                                                    acccgggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc       <  1:1874073/96‑1 (MQ=255)
                                                                                                       ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCCC.CCC.CC.TTCGCGGGCGGCCGCGACCCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa    <  2:817038/101‑1 (MQ=255)
                                                                                                         |                                                                                                          
CGGGGCAAGAAGGGGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACC  >  NC_014034/2249260‑2249457

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