Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,374 +G 2.4% intergenic (‑170/‑135) 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,3742.→G 2.4% 23.1 ‑1.5 ‑0.8466intergenic (‑170/‑135)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (8/3):  ref base (233/222):  total (241/225)
Fisher's exact test strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.72e-01

Only 100 of 503 total aligned reads displayed.

GCTCTGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.CCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGC  >  NC_014034/2249275‑2249473
                                                                                                           |                                                                                                        
  ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                                                                      >  1:7867319/5‑101 (MQ=255)
    gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                    <  1:8949544/99‑1 (MQ=255)
    gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                    >  1:5875644/3‑101 (MQ=255)
    gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                    >  2:6671115/3‑101 (MQ=255)
    gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                    >  1:2691587/3‑101 (MQ=255)
     ccccgGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCG.AAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                  <  1:3658952/100‑1 (MQ=255)
      ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                  <  2:8198707/101‑1 (MQ=255)
          gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                              >  2:2094068/1‑101 (MQ=255)
          gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                              >  1:3172563/1‑101 (MQ=255)
          gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                              >  2:9139202/1‑101 (MQ=255)
           tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGCTGAA...GGGG..GGCTTTCATCTTGCCGCAAATATCCCGGGGGGGGAA.TGGCGCC..GCCCGAgggggg                                                                                             >  1:2967431/1‑101 (MQ=37)
           tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                             >  2:4126638/1‑101 (MQ=255)
           tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                             <  1:4053185/101‑1 (MQ=255)
            cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.c                                                                                           >  1:1388815/1‑101 (MQ=255)
             ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.Ca                                                                                          <  2:1900228/101‑1 (MQ=255)
              cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg                                                                                         >  2:8447664/1‑101 (MQ=255)
              cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg                                                                                         >  1:4716201/1‑101 (MQ=255)
               gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg.c                                                                                       <  2:5204198/101‑1 (MQ=255)
               gtGACGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg.c                                                                                       <  1:2906103/101‑1 (MQ=255)
               gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                                             >  1:7916613/7‑101 (MQ=21)
                caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGG.c                                                                                           >  1:4364353/6‑101 (MQ=18)
                caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATACCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGG.c                                                                                           <  1:6454527/96‑1 (MQ=18)
                 gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg.cg.c                                                                                    >  2:5365/1‑101 (MQ=255)
                   ctggACT.CCCCCGGAATATTTGAGGCGAGATGAA...GGGG..GGCGTTCATCTTGCCGAATATATCCCGGGTGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAT.CG.ccc                                                                                  <  2:1110049/99‑1 (MQ=255)
                    gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.c                                                                                <  1:226673/101‑1 (MQ=255)
                    gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.c                                                                                <  2:7563904/101‑1 (MQ=255)
                     ggACT.CCCCCGGGATATTTGAGGCAGGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.cc                                                                               <  2:1779782/101‑1 (MQ=255)
                      gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.ccc                                                                              >  2:3306408/1‑100 (MQ=255)
                        cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.ccctt                                                                            <  2:1484146/101‑4 (MQ=255)
                        cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.ccctt                                                                            <  2:2795914/101‑4 (MQ=255)
                        cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.ccctt                                                                            <  1:2422173/101‑4 (MQ=255)
                         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.cccttc                                                                           <  2:1281900/101‑5 (MQ=255)
                         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc.cccttc                                                                           <  1:771031/101‑5 (MQ=255)
                           cccccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGG.CAG.CG.CCC.CCt                                                                              <  1:1892457/101‑1 (MQ=0)
                               cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGG.CAG.CG.CCC.CCtgccc                                                                          >  1:7636372/1‑97 (MQ=0)
                                  gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCgcgggcgg                                                                    >  1:5603106/1‑101 (MQ=255)
                                    tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGcc                                                                  >  2:3705213/1‑101 (MQ=255)
                                    tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGcc                                                                  >  2:5563506/1‑101 (MQ=255)
                                      tttGGGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCcgc                                                                <  2:5820082/101‑1 (MQ=255)
                                      tttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCcgc                                                                <  2:1360245/101‑1 (MQ=255)
                                         gAGGCAAGATGAA...GGGG..GTCTTTCATCTTGCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGAc                                                             <  2:9023061/101‑1 (MQ=37)
                                         gAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGAc                                                             <  2:2203705/101‑1 (MQ=255)
                                         gAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGAc                                                             >  1:8079337/1‑101 (MQ=255)
                                          aGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGAcg                                                            >  2:6290077/1‑101 (MQ=255)
                                             cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGAcgcgc                                                         >  2:2054705/1‑101 (MQ=255)
                                             cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGAcgcgc                                                         <  2:3021647/101‑1 (MQ=255)
                                                   gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   <  2:8915125/101‑1 (MQ=255)
                                                   gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   >  1:3299225/1‑101 (MQ=255)
                                                    ga...GGGG..GGCTTTCATCGTGCGGAAAATATCCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                  <  1:340630/100‑1 (MQ=37)
                                                     a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                 <  2:4794888/101‑1 (MQ=255)
                                                     a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                 <  2:5049705/101‑1 (MQ=255)
                                                        gaggggagtcTTTCATCGTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   <  2:3529689/92‑1 (MQ=25)
                                                          ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                               <  2:6919664/101‑1 (MQ=255)
                                                           gg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                              <  1:6664534/101‑1 (MQ=255)
                                                            g..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAggc                                             >  1:5012030/1‑100 (MQ=255)
                                                               ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCCCGACGCCCGCCGCAACAGGTc                                            >  1:104517/1‑101 (MQ=25)
                                                                gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                                           >  2:4727696/1‑101 (MQ=255)
                                                                 cTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                          <  2:2398581/101‑1 (MQ=255)
                                                                      aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                     <  1:5708370/101‑1 (MQ=255)
                                                                       tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                    >  1:7422879/1‑101 (MQ=255)
                                                                       tCTGGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                    <  1:5350800/101‑1 (MQ=255)
                                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                   >  2:3901468/1‑101 (MQ=255)
                                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                   >  2:6978659/1‑101 (MQ=255)
                                                                        cTTGCCAAAAATGTCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CCG.GGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                   <  2:729945/101‑1 (MQ=255)
                                                                        cTTACCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGTGGCCGCGACGCGAGCCGCAACAAGGCGTGCCCCtt                                   >  1:8804130/1‑101 (MQ=25)
                                                                          tGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                 <  2:7402025/101‑1 (MQ=255)
                                                                             cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                              >  2:2530373/1‑101 (MQ=255)
                                                                              gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                             <  1:9195375/101‑1 (MQ=255)
                                                                              gAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                             >  2:6147952/1‑101 (MQ=255)
                                                                               aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                            >  2:6523545/1‑101 (MQ=255)
                                                                               aaaaTATCCCGGGGGGTGAA.TAGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                            >  1:999284/1‑101 (MQ=255)
                                                                                aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                           <  2:7906586/101‑1 (MQ=255)
                                                                                 aaTATCCCGGGGGGTGAA.TGTCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                          <  1:7204665/101‑1 (MQ=255)
                                                                                 aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGAGGG.CG.CCC.CTTCCGCGGG.GGGCGCGACGCGCGCCGCAACAGGTCCGGCCCCTTTGTCAttttt                         >  2:4449441/1‑98 (MQ=25)
                                                                                 aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                          <  1:6897623/101‑1 (MQ=255)
                                                                                 aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                          <  1:168561/101‑1 (MQ=255)
                                                                                  atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.CCCACCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGGGCCCCTTTGTCATTGTg                         >  1:4225264/1‑101 (MQ=38)
                                                                                   tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                        <  2:2669662/101‑1 (MQ=255)
                                                                                   tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                        >  1:7581512/1‑101 (MQ=255)
                                                                                       ccGGGGGGTGAA.TGGCGTA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                    >  1:5448480/1‑101 (MQ=255)
                                                                                       ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                    >  1:1193650/1‑101 (MQ=255)
                                                                                        cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                   >  2:1617854/1‑101 (MQ=255)
                                                                                        cGGGGGCTGAA.TGGGGAA..GCTAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                   <  2:5507599/101‑1 (MQ=255)
                                                                                         ggggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                  <  1:3820811/101‑1 (MQ=255)
                                                                                           ggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                <  1:6292409/101‑1 (MQ=255)
                                                                                            gggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt               >  2:5881813/1‑101 (MQ=255)
                                                                                            ggcgcgccgggCGCG..CCCCGAGGGGGG.CAGCCC.CCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                <  2:8074864/92‑1 (MQ=255)
                                                                                                gAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa           <  2:7710173/101‑1 (MQ=255)
                                                                                                 aa.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt          >  2:6033012/1‑101 (MQ=255)
                                                                                                  a.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa         >  2:7149675/1‑101 (MQ=255)
                                                                                                  a.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa         <  1:2049720/101‑1 (MQ=255)
                                                                                                  a.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa         <  2:8597677/101‑1 (MQ=255)
                                                                                                   agggCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa         <  1:6240183/99‑1 (MQ=255)
                                                                                                    tGGCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa        >  2:4665931/1‑101 (MQ=255)
                                                                                                     ggCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg       <  1:4831057/101‑1 (MQ=255)
                                                                                                     ggCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg       <  1:1532694/101‑1 (MQ=255)
                                                                                                     ggCGCA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg       <  1:3533787/101‑1 (MQ=255)
                                                                                                       cgcA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt     <  1:8438038/101‑1 (MQ=255)
                                                                                                         cA..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg   <  2:6621588/101‑1 (MQ=255)
                                                                                                          a..GCCAGAGGGGGG.CAG.CGCCCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGCCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc  >  1:8536939/1‑101 (MQ=255)
                                                                                                           |                                                                                                        
GCTCTGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.CCC.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGC  >  NC_014034/2249275‑2249473

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