Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,330 G→A 4.2% intergenic (‑126/‑179) 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,3300G→A 4.2% 16.0 ‑0.7 ‑0.0236intergenic (‑126/‑179)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/6):  ref base (124/102):  total (128/108)
Fisher's exact test strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.70e-01

Only 100 of 273 total aligned reads displayed.

GATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTG.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAG  >  NC_014034/2249232‑2249428
                                                                                                         |                                                                                                          
 aTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg..g                                                                                                        <  1:4093201/101‑1 (MQ=255)
 aTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg..g                                                                                                        >  1:341858/1‑101 (MQ=255)
   gCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGc                                                                                                      <  2:4448752/101‑1 (MQ=255)
    cacaGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCt                                                                                                     <  1:1417376/101‑1 (MQ=255)
    cacaGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCt                                                                                                     <  2:1330604/101‑1 (MQ=255)
       aGATACCGCAGCTTGGGCACCCGGGGCAAGNAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTc                                                                                                  <  1:975538/101‑1 (MQ=255)
       aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTc                                                                                                  >  2:1270258/1‑101 (MQ=255)
       aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTc                                                                                                  <  1:3049482/101‑1 (MQ=255)
       aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTc                                                                                                  <  1:3110599/101‑1 (MQ=255)
        gATACCGCAGCTTGGGCNCCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCa                                                                                                 <  2:2347573/101‑1 (MQ=255)
        gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.TCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCa                                                                                                 <  2:2222766/101‑1 (MQ=255)
        gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCa                                                                                                 <  2:613945/101‑1 (MQ=255)
         aTACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCAt                                                                                                >  2:4133258/1‑101 (MQ=255)
           aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                              >  2:4356992/1‑101 (MQ=255)
            ccGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCtt                                                                                             <  1:1092739/101‑1 (MQ=255)
             cGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTg                                                                                            <  2:3337026/101‑1 (MQ=255)
             cGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTg                                                                                            >  1:4236661/1‑101 (MQ=255)
              gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGc                                                                                           >  1:3293908/1‑101 (MQ=255)
              gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGc                                                                                           <  1:3890979/101‑1 (MQ=255)
              gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGc                                                                                           >  1:3509049/1‑101 (MQ=255)
                aGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCg                                                                                         <  2:1443510/101‑1 (MQ=255)
                  cTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGCTATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                       >  1:1273009/1‑101 (MQ=255)
                  cTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                       >  1:4147154/1‑101 (MQ=255)
                   ttGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaa                                                                                      >  1:3827198/1‑101 (MQ=255)
                    tGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGGTTTCATCTTGCCGaaaa                                                                                     >  2:4819187/1‑101 (MQ=255)
                    tGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                     <  2:4579861/101‑1 (MQ=255)
                     gggCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                    >  1:3198249/1‑101 (MQ=255)
                     gggCACCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                    <  2:322504/101‑1 (MQ=255)
                         aCCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATcc                                                                                <  1:431164/101‑1 (MQ=255)
                         aCCCGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATcc                                                                                <  2:157007/101‑1 (MQ=255)
                          cccGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATccc                                                                               >  2:3515121/1‑101 (MQ=255)
                          cactcggCAAGAAGCGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATacc                                                                               >  2:3212936/6‑101 (MQ=37)
                            cGGGGCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                             >  1:1218484/1‑101 (MQ=255)
                              gggCAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggg                                                                           >  2:4510651/1‑101 (MQ=255)
                                 cAAGAAGGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGt                                                                        <  2:4790983/101‑1 (MQ=255)
                                  gagaagGGGGCTGCGC.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGTGTg                                                                       <  2:1110866/100‑1 (MQ=38)
                                  aagaagGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTg                                                                       >  2:365312/1‑101 (MQ=255)
                                     aagGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.t                                                                   <  1:1933932/101‑1 (MQ=255)
                                     aagGGGGCT.CTG.CCCACGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.t                                                                   >  2:3903940/1‑101 (MQ=255)
                                      agGGGGCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.Tg                                                                  >  1:25854/1‑101 (MQ=255)
                                        ggggCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgc                                                                <  2:2188273/101‑1 (MQ=255)
                                         gggCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcg                                                               <  2:4028024/101‑1 (MQ=255)
                                         gggCT.CTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcg                                                               <  1:255547/101‑1 (MQ=255)
                                            ct.ctGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..g                                                          <  1:4295982/101‑1 (MQ=255)
                                            ct.ctGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..g                                                          <  1:1751733/101‑1 (MQ=255)
                                               ctgcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCNNNNNTATCCCGGGGGTTGAA.TCGCGCA..Gcc                                                        >  1:3250151/5‑101 (MQ=25)
                                               ctgcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                        <  2:4878389/97‑1 (MQ=255)
                                                 gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                      >  1:3994163/3‑101 (MQ=255)
                                                 gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                      >  1:3177311/3‑101 (MQ=255)
                                                 gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGAAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATACCGGGGGGGGAA.GGGCGCA..GCCag                                                      >  1:2491384/3‑101 (MQ=38)
                                                  cccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                     <  1:2904333/100‑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..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                >  1:4881073/1‑101 (MQ=255)
                                                        gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGAGAA.TGGCGCA..GCCAGAggggg                                                >  2:1813618/1‑101 (MQ=255)
                                                         tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                               <  1:392326/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..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                              <  1:4980737/101‑1 (MQ=255)
                                                          cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGc                                              >  2:2443535/1‑101 (MQ=255)
                                                           ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                             >  2:3711690/1‑101 (MQ=255)
                                                           ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GACAGAGGGGGGCa                                             <  1:2728576/101‑1 (MQ=255)
                                                            cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                            <  1:675225/101‑1 (MQ=255)
                                                             gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggggagat                                           >  2:2592967/1‑97 (MQ=255)
                                                             gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAgc                                           >  1:4454983/1‑101 (MQ=255)
                                                               gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAgcg.c                                        >  2:2661997/1‑101 (MQ=255)
                                                                 gcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAgc                                           <  1:2286656/99‑1 (MQ=18)
                                                                 gcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAgc                                           <  1:2459501/99‑1 (MQ=18)
                                                                 cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccc                                      >  1:1052733/1‑101 (MQ=255)
                                                                   ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGTCAGCG.ccccc                                    <  1:4368513/101‑1 (MQ=255)
                                                                    gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.cccccc                                   <  2:270951/101‑2 (MQ=255)
                                                                    gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.cccccc                                   <  1:4752675/101‑2 (MQ=255)
                                                                     aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.cccccct                                  >  1:485248/1‑99 (MQ=255)
                                                                     aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.cccccct                                  >  1:2140057/1‑99 (MQ=255)
                                                                     aCT.CCCCAGGGTTATTTGAGGCAAGATGTA...GGGG..GGCTTTCATCTTGGCGAAAATTTGCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.cccccct                                  <  2:4101922/101‑3 (MQ=25)
                                                                       t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccccccttc                                <  2:1645830/101‑5 (MQ=255)
                                                                       t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccccccttc                                <  1:2181228/101‑5 (MQ=255)
                                                                       t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.ccccccttc                                >  1:4658426/1‑97 (MQ=255)
                                                                           cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAGCG.CCCCC..tgc                               <  2:1204802/101‑3 (MQ=0)
                                                                              gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCgcgggc                         <  1:654632/101‑1 (MQ=255)
                                                                               ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCgcggggg                        >  1:3003233/1‑99 (MQ=255)
                                                                                gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCgcgggcgg                       >  1:4569746/1‑101 (MQ=255)
                                                                                 atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTAAA.TGGCGCA..GACAGAGGGGGGCAGCGCCCCCC..TTCGCGggcggc                      >  2:901926/1‑101 (MQ=37)
                                                                                 atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGggcggc                      >  1:2251467/1‑101 (MQ=255)
                                                                                 atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGggcggc                      >  1:74012/1‑101 (MQ=255)
                                                                                     ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCcgcg                  >  2:535086/1‑101 (MQ=255)
                                                                                     ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCcgcg                  >  2:508930/1‑101 (MQ=255)
                                                                                              gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.CCCCCA.TTCGCGGGCGGCCGCGACGCGcgcc         >  2:225834/1‑101 (MQ=255)
                                                                                               gagga...GGGG..GGCTTTCATTTTGGCAAAAATGCCCCGGGGGGGTAA.TGGGGGC..GCCCGGGGGGGGGCGCC.CCCCCCTTTCGCGGGCGGCCGCGACGCGcgcc         <  2:4236344/97‑1 (MQ=255)
                                                                                                tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGcgccgc       >  2:1777897/1‑101 (MQ=255)
                                                                                                 gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCa      >  1:3437886/1‑101 (MQ=255)
                                                                                                 gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCa      <  2:4682148/101‑1 (MQ=255)
                                                                                                  aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGGCGCGCGCCGCaa     >  2:1923001/1‑101 (MQ=255)
                                                                                                  aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa     >  2:1052694/1‑101 (MQ=255)
                                                                                                  aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa     >  1:2076968/1‑101 (MQ=255)
                                                                                                   a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc    <  1:2855191/101‑1 (MQ=255)
                                                                                                   a...GGGG..GGCTTTCATCTTGCCAAAAATGTCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc    <  2:396183/101‑1 (MQ=37)
                                                                                                       gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa   <  1:365365/101‑1 (MQ=255)
                                                                                                       gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa   <  2:2897380/101‑1 (MQ=255)
                                                                                                       gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa   <  2:349745/101‑1 (MQ=255)
                                                                                                         |                                                                                                          
GATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTG.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCG.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAG  >  NC_014034/2249232‑2249428

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