Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,333 G→C 2.5% intergenic (‑129/‑176) 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,3330G→C 2.5% 5.2 0.0 ‑0.0253intergenic (‑129/‑176)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (3/4):  ref base (111/134):  total (114/139)
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 284 total aligned reads displayed.

TGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTG.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCG  >  NC_014034/2249234‑2249432
                                                                                                           |                                                                                                        
tGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg..gg                                                                                                          >  1:3335318/1‑101 (MQ=255)
tGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg..gg                                                                                                          <  2:1659645/101‑1 (MQ=255)
 gCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGc                                                                                                         >  2:5463407/1‑101 (MQ=255)
 gCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGc                                                                                                         >  1:3581464/1‑101 (MQ=255)
  cacaGATACCGCATCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCt                                                                                                        <  1:2854235/101‑1 (MQ=255)
  cacaGAAACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCt                                                                                                        <  1:959722/101‑1 (MQ=255)
    caGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCttt                                                                                                      >  1:5615551/1‑101 (MQ=255)
    caGATACCGCAGCTTGGCCCCCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCttt                                                                                                      <  2:2464731/101‑1 (MQ=38)
       aTACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCAt                                                                                                   >  2:49075/1‑101 (MQ=255)
       aTACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCAt                                                                                                   <  2:1781233/101‑1 (MQ=255)
         aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                                 <  1:5029491/101‑1 (MQ=255)
         aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                                 >  2:2059805/1‑101 (MQ=255)
         aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                                 <  2:4470016/101‑1 (MQ=255)
           cGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTg                                                                                               <  2:2963652/101‑1 (MQ=255)
           cGCAGCTGGGGCCCCGTGGGCAAGAAGGGGGATCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTg                                                                                               <  2:804596/101‑1 (MQ=37)
            gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGc                                                                                              >  1:4278022/1‑101 (MQ=255)
            gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGc                                                                                              <  1:5171188/101‑1 (MQ=255)
             cAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGcc                                                                                             >  1:1758832/1‑101 (MQ=255)
             cAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGcc                                                                                             >  2:4923532/1‑101 (MQ=255)
               gCTTGGGCACCCGGGGCAAGAAGGGGCCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGa                                                                                           <  2:3764183/101‑1 (MQ=255)
                cTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                          >  2:1073681/1‑101 (MQ=255)
                 ttGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaa                                                                                         <  2:507003/101‑1 (MQ=255)
                  tGGGCACGCGGGGCAAGAAATGGGCTCTGCCCC.CGTCCCGTGGCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGAG..GGCTTTCATCTTGCCGaaaa                                                                                        <  1:3279869/101‑1 (MQ=25)
                  tGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                        <  1:3061319/101‑1 (MQ=255)
                  tGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                        >  2:3555088/1‑101 (MQ=255)
                  tGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                        <  2:4824643/101‑1 (MQ=255)
                   gggCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                       <  2:3637712/101‑1 (MQ=255)
                     gCACCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAatat                                                                                     >  1:635718/1‑101 (MQ=255)
                       aCCCGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTCCCGAAAATATcc                                                                                   <  1:3126609/101‑1 (MQ=255)
                        cccGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATccc                                                                                  <  1:2900064/101‑1 (MQ=255)
                         ccGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCg                                                                                 <  1:5313274/101‑1 (MQ=255)
                          cGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                                >  2:3257542/1‑101 (MQ=255)
                          cGGGGCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                                <  1:4216584/101‑1 (MQ=255)
                             ggCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCggggg                                                                             <  2:4536711/101‑1 (MQ=255)
                              gCAAGAAGGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                            >  1:1531085/1‑101 (MQ=255)
                               gaagaagGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGt                                                                           <  1:1483221/100‑1 (MQ=255)
                                aagaagGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggggg                                                                          >  1:1957463/1‑99 (MQ=255)
                                aagaagGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggggg                                                                          >  1:2575173/1‑99 (MQ=255)
                                aagaagGGGGCTCTGCCCC.CGTCCCGGGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTg                                                                          <  1:2059922/101‑1 (MQ=255)
                                 agaagGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGa                                                                         <  1:2353482/101‑1 (MQ=255)
                                   aatggggCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.t                                                                      >  2:4143474/4‑101 (MQ=255)
                                   aagGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.t                                                                      >  2:1497481/1‑101 (MQ=255)
                                   aagGGGGCTCTGCCCC.CGTCCCGTGCCCGCACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.t                                                                      <  1:2256815/101‑1 (MQ=37)
                                    agGGGGCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.Tg                                                                     >  2:2174038/1‑101 (MQ=255)
                                      ggggCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgc                                                                   <  2:5601140/101‑1 (MQ=255)
                                      ggggCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgc                                                                   <  2:4654860/101‑1 (MQ=255)
                                        ggCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcgc                                                                 <  1:3463851/101‑1 (MQ=255)
                                        ggCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcgc                                                                 >  2:5143000/1‑101 (MQ=255)
                                          ctctGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..g                                                             >  2:1104055/1‑101 (MQ=255)
                                            ctG.CCCGCGTCCCGGGCCGGGACT.CCCGCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                           <  1:3555305/101‑1 (MQ=37)
                                             tgcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCa                                                          <  2:3915410/98‑1 (MQ=255)
                                              gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCGAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcccg                                                         >  2:3091455/3‑99 (MQ=255)
                                              gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                         >  1:3687868/3‑101 (MQ=255)
                                              gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                         >  2:3855116/3‑101 (MQ=255)
                                               cccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                        <  2:2317432/100‑1 (MQ=255)
                                                  cacatccagtACCGGGACT.CCCCCGGGATATTTGAGGCAAGCTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgg                                                      <  2:5704248/97‑1 (MQ=37)
                                                    cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                    >  2:475349/1‑101 (MQ=255)
                                                    cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                    >  2:251324/1‑101 (MQ=255)
                                                     gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                   <  1:251357/101‑1 (MQ=255)
                                                      tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                  >  1:4342072/1‑101 (MQ=255)
                                                        ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCCGAGGGGGGCa                                                >  2:2572536/1‑101 (MQ=255)
                                                        ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                <  1:4941439/101‑1 (MQ=255)
                                                        ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                >  2:2728021/1‑101 (MQ=255)
                                                        ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                >  1:3778898/1‑101 (MQ=255)
                                                         cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                               <  2:5103290/101‑1 (MQ=255)
                                                          gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                             <  1:2671209/101‑1 (MQ=255)
                                                          gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                             >  2:4357325/1‑101 (MQ=255)
                                                          gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                  <  2:5238796/95‑1 (MQ=18)
                                                             ccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc                                         <  2:171451/101‑1 (MQ=255)
                                                                ggccT.CCCCCGGGATATTGGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                      <  1:1024759/98‑1 (MQ=255)
                                                                ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                      <  2:5644335/101‑1 (MQ=255)
                                                                ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                      >  2:4023660/1‑101 (MQ=255)
                                                                ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                      >  2:4051320/1‑101 (MQ=255)
                                                                 gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccc                                     >  1:1672024/1‑100 (MQ=255)
                                                                        cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCCtgc                                   >  1:1041015/1‑99 (MQ=2)
                                                                        cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCCtgc                                   >  2:213729/1‑99 (MQ=2)
                                                                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggc                             >  1:4592779/1‑101 (MQ=255)
                                                                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggc                             <  2:4246316/101‑1 (MQ=255)
                                                                            ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggcg                            >  1:1103886/1‑101 (MQ=255)
                                                                             gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggcgg                           >  1:4712468/1‑101 (MQ=255)
                                                                               ttttggAGGAAAGAAAAG...GGGG..GGTTTTAATCTTGCCAAAAATATCCGGGGGGGTGAA.TGGCGCA..GCAAGGGGGGGGCAGCCC.CCCCCTTCGCGGGCGGcc                         <  2:1175788/96‑1 (MQ=255)
                                                                               tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGcc                         <  1:223948/101‑1 (MQ=255)
                                                                                   tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGa                     >  1:1217705/1‑101 (MQ=255)
                                                                                     aGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcg                   >  2:4950824/1‑101 (MQ=255)
                                                                                       gCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcg                 <  2:2212496/101‑1 (MQ=255)
                                                                                        cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgc                >  1:4014290/1‑101 (MQ=255)
                                                                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg               <  1:553878/101‑1 (MQ=255)
                                                                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg               >  1:85512/1‑101 (MQ=255)
                                                                                          aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACgcgcgc              >  1:5437517/1‑101 (MQ=255)
                                                                                           gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgcc             >  1:4795376/1‑101 (MQ=255)
                                                                                           gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgcc             <  2:3335113/101‑1 (MQ=255)
                                                                                               aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCaa         >  1:3178166/1‑101 (MQ=255)
                                                                                               aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCaa         >  1:1283618/1‑101 (MQ=255)
                                                                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc        <  2:1951910/101‑1 (MQ=255)
                                                                                                 gatggtg..ggCTTTCAACTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa          <  2:744295/95‑1 (MQ=25)
                                                                                                     ggg..ggTTTTCATCTTGCCGAAAAAACCCGGGGGGGTGAA.TGGCGCA..GCCAGGGGGGGCCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg      <  1:4228043/101‑1 (MQ=37)
                                                                                                      gg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg     <  2:284442/101‑1 (MQ=255)
                                                                                                       g..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGGGCA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCTCGGGCGGCAGGGACGCGGGACGCAACAgga    >  2:4375029/1‑100 (MQ=21)
                                                                                                          ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc   <  2:5398942/101‑1 (MQ=255)
                                                                                                           gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCACCAGGTCg  >  1:1329735/1‑101 (MQ=255)
                                                                                                           |                                                                                                        
TGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTG.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCG  >  NC_014034/2249234‑2249432

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG