Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,332 G→T 2.5% intergenic (‑128/‑177) 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,3320G→T 2.5% 6.2 ‑1.5 ‑0.5236intergenic (‑128/‑177)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (123/107):  total (128/108)
Fisher's exact test strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.88e-01

Only 100 of 270 total aligned reads displayed.

ATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAG  >  NC_014034/2249233‑2249428
                                                                                                           |                                                                                                         
aTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg..g                                                                                                           >  1:341858/1‑101 (MQ=255)
  gCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGc                                                                                                         >  1:4565373/1‑101 (MQ=255)
   cacaGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.t                                                                                                       <  2:1330604/101‑1 (MQ=255)
   cacaGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.t                                                                                                       <  1:1417376/101‑1 (MQ=255)
      aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTc                                                                                                    >  2:1270258/1‑101 (MQ=255)
      aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTc                                                                                                    <  1:998889/101‑1 (MQ=255)
      aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTc                                                                                                    <  1:3110599/101‑1 (MQ=255)
       gATACCGCAGCTTGGGCNCCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCa                                                                                                   <  2:2347573/101‑1 (MQ=255)
       gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.TCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCa                                                                                                   <  2:2222766/101‑1 (MQ=255)
       gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCa                                                                                                   <  2:613945/101‑1 (MQ=255)
           ccGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCtt                                                                                               <  1:1092739/101‑1 (MQ=255)
            cGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.g                                                                                             <  2:3337026/101‑1 (MQ=255)
            cGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.g                                                                                             >  1:4236661/1‑101 (MQ=255)
             gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.Gc                                                                                            <  1:3890979/101‑1 (MQ=255)
             gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.Gc                                                                                            >  1:3509049/1‑101 (MQ=255)
             gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.Gc                                                                                            >  1:3293908/1‑101 (MQ=255)
               aGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCg                                                                                          <  2:1443510/101‑1 (MQ=255)
               aGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCg                                                                                          <  1:660105/101‑1 (MQ=255)
                 cTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGCTATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGaa                                                                                        >  1:1273009/1‑101 (MQ=255)
                 cTTGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGaa                                                                                        >  1:4147154/1‑101 (MQ=255)
                  ttGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGaaa                                                                                       >  1:3827198/1‑101 (MQ=255)
                  ttGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGaaa                                                                                       >  2:1220639/1‑101 (MQ=255)
                   tGGGCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGaaaa                                                                                      <  2:4579861/101‑1 (MQ=255)
                    gggCACCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAat                                                                                     <  2:322504/101‑1 (MQ=255)
                        aCCCGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATcc                                                                                 <  1:431164/101‑1 (MQ=255)
                         cccGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATccc                                                                                >  2:3515121/1‑101 (MQ=255)
                           cGGGGCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCgg                                                                              >  1:1218484/1‑101 (MQ=255)
                             gggCAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCgggg                                                                            >  2:4510651/1‑101 (MQ=255)
                                cAAGAAGGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGt                                                                         <  2:4790983/101‑1 (MQ=255)
                                 gagaagGGGGCTGCGC.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGTGTg                                                                        <  2:1110866/100‑1 (MQ=38)
                                 aagaagGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTg                                                                        >  2:4681391/1‑101 (MQ=255)
                                 aagaagGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTg                                                                        >  2:365312/1‑101 (MQ=255)
                                    aagGGGGCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.t                                                                    <  1:1933932/101‑1 (MQ=255)
                                       ggggCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGgc                                                                 <  2:2188273/101‑1 (MQ=255)
                                        gggCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGgcg                                                                <  1:255547/101‑1 (MQ=255)
                                        gggCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGgcg                                                                <  2:4028024/101‑1 (MQ=255)
                                         ggCT.CTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGgcgc                                                               >  1:2439396/1‑101 (MQ=255)
                                           ct.ctGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..g                                                           <  1:4295982/101‑1 (MQ=255)
                                           ct.ctGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..g                                                           <  1:1751733/101‑1 (MQ=255)
                                              ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCNNNNNTATCCCGGGGGTTGAA.TCGCGCA..Gcc                                                         >  1:3250151/5‑101 (MQ=25)
                                              ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                         <  2:4878389/97‑1 (MQ=255)
                                              ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..Gcc                                                         >  1:3458951/5‑101 (MQ=255)
                                                gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                       >  1:3177311/3‑101 (MQ=255)
                                                gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                       >  1:3994163/3‑101 (MQ=255)
                                                gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGAAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATACCGGGGGGGGAA.GGGCGCA..GCCag                                                       >  1:2491384/3‑101 (MQ=38)
                                                 cccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                      <  1:2904333/100‑1 (MQ=255)
                                                  ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                     >  1:817140/1‑101 (MQ=255)
                                                     cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                  >  1:2776911/1‑101 (MQ=255)
                                                     cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAACATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                  >  2:579604/1‑101 (MQ=255)
                                                      gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                 >  1:3213276/1‑101 (MQ=255)
                                                       tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                >  2:4829050/1‑101 (MQ=255)
                                                       tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                <  1:392326/101‑1 (MQ=255)
                                                       atctgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                <  1:2593148/97‑1 (MQ=255)
                                                        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                               >  1:2415653/1‑101 (MQ=255)
                                                        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                               <  1:4980737/101‑1 (MQ=255)
                                                        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                               <  1:2111314/101‑1 (MQ=255)
                                                        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGc                                               >  2:2443535/1‑101 (MQ=255)
                                                         ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                              >  2:3711690/1‑101 (MQ=255)
                                                         ccgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAAC.CGGA..TGCATTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggg                                                  >  1:4872835/9‑101 (MQ=21)
                                                          cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                             <  1:675225/101‑1 (MQ=255)
                                                          cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                             <  2:4307761/101‑1 (MQ=255)
                                                           gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggggagat                                            >  2:2592967/1‑97 (MQ=255)
                                                           gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                           >  1:2566735/1‑101 (MQ=255)
                                                           gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                           >  1:4454983/1‑101 (MQ=255)
                                                            tGCCTGGACT.CCCCCGGGATATTTGAGGCAAGATCAG...GGGG..GGT.TTTCATCTT.TCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                          <  1:2240441/101‑1 (MQ=37)
                                                             gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                        >  2:2661997/1‑101 (MQ=255)
                                                               gcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c                                           <  1:2286656/99‑1 (MQ=18)
                                                               cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc                                      >  1:1052733/1‑101 (MQ=255)
                                                                 ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.cg.c                                        <  2:3846146/101‑1 (MQ=17)
                                                                 ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGTCAG.CG.ccccc                                    <  1:4368513/101‑1 (MQ=255)
                                                                  gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccc                                   <  2:270951/101‑2 (MQ=255)
                                                                  gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccc                                   <  1:4752675/101‑2 (MQ=255)
                                                                   aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccct                                  >  1:2140057/1‑99 (MQ=255)
                                                                    cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccctt                                 <  1:4041099/101‑4 (MQ=255)
                                                                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                <  1:2181228/101‑5 (MQ=255)
                                                                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                <  2:1645830/101‑5 (MQ=255)
                                                                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                >  1:4658426/1‑97 (MQ=255)
                                                                           cGGGATATTTGCGGCAAGATGAA...GGGGCATGA.TTTCATCTT.GCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCC..tgccc                             >  2:3195667/1‑97 (MQ=2)
                                                                             ggATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCgcggggg                        >  1:3003233/1‑99 (MQ=255)
                                                                               atatTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTAAA.TGGCGCA..GACAGAGGGGGGCAG.CGCCCCCC..TTCGCGggcggc                      >  2:901926/1‑101 (MQ=37)
                                                                                 atTTGCGGCAAGATGAA...GGGGCATGA.TTTCATCTT.GCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCC..tgccccgccgg                       >  2:4122224/1‑91 (MQ=2)
                                                                                   ttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCcgcg                  >  2:535086/1‑101 (MQ=255)
                                                                                   ttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCcgcg                  >  2:508930/1‑101 (MQ=255)
                                                                                     gAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGAc                >  2:4379902/1‑101 (MQ=255)
                                                                                      aggtaggaggAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGAcg               <  2:1582857/92‑1 (MQ=25)
                                                                                       ggCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGAcgc              <  2:3458686/101‑1 (MQ=255)
                                                                                            gATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCCA.TTCGCGGGCGGCCGCGACGCGcgcc         >  2:225834/1‑101 (MQ=255)
                                                                                             gagga...GGGG..GGC.TTTCATTTT.GGCAAAAATGCCCCGGGGGGGTAA.TGGGGGC..GCCCGGGGGGGGGCG.CC.CCCCCCTTTCGCGGGCGGCCGCGACGCGcgcc         <  2:4236344/97‑1 (MQ=255)
                                                                                              tGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGcgccgc       >  2:1777897/1‑101 (MQ=255)
                                                                                                ta...GGGG..GGG.CTTCATATTCCCCGAAAATAGCGTGGGGGGTAAA.CGGCGCA..GCCAGGGGGGGGCGGCCC.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCa      <  1:3078535/100‑1 (MQ=255)
                                                                                                aa...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGGCGCGCGCCGCaa     >  2:1923001/1‑101 (MQ=255)
                                                                                                aa...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa     >  2:1052694/1‑101 (MQ=255)
                                                                                                 a...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc    <  1:2855191/101‑1 (MQ=255)
                                                                                                 a...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc    >  2:4074801/1‑101 (MQ=255)
                                                                                                 a...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc    >  2:1840746/1‑101 (MQ=255)
                                                                                                 a...GGGG..GGC.TTTCATCTT.GCCAAAAATGTCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc    <  2:396183/101‑1 (MQ=37)
                                                                                                     gggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa   <  1:365365/101‑1 (MQ=255)
                                                                                                     gggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa   <  2:2897380/101‑1 (MQ=255)
                                                                                                     gggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa   <  1:4806444/101‑1 (MQ=255)
                                                                                                      ggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg  >  1:2062827/1‑101 (MQ=255)
                                                                                                           |                                                                                                         
ATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCT.CTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAG  >  NC_014034/2249233‑2249428

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