Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,381 G→A 2.7% intergenic (‑177/‑128) 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,3810G→A 2.7% 6.6 ‑0.4 ‑0.1222intergenic (‑177/‑128)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (118/98):  total (122/100)
Fisher's exact test strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.66e-01

Only 100 of 270 total aligned reads displayed.

CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCG  >  NC_014034/2249281‑2249480
                                                                                                              |                                                                                                         
ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                           >  1:817140/1‑101 (MQ=255)
   cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                                        >  1:2776911/1‑101 (MQ=255)
     tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                                      <  2:4118106/101‑1 (MQ=255)
     tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                                      <  1:392326/101‑1 (MQ=255)
     atctgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                                      <  1:2593148/97‑1 (MQ=255)
      cccgtGCCGGGACT.CCCCGGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                     <  2:107116/101‑1 (MQ=255)
      cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                     >  1:2415653/1‑101 (MQ=255)
      cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                     >  2:2443535/1‑101 (MQ=255)
       ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                    >  2:3711690/1‑101 (MQ=255)
        cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   <  1:3312132/101‑1 (MQ=255)
        cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   <  1:3498506/101‑1 (MQ=255)
        cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   <  1:3904232/101‑1 (MQ=255)
        cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   <  1:675225/101‑1 (MQ=255)
        cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   <  2:4307761/101‑1 (MQ=255)
        cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                                       >  2:1685360/8‑101 (MQ=21)
          caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                                                     <  2:902873/96‑1 (MQ=18)
           gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAg.ag..c                                                                                             >  2:1737242/1‑101 (MQ=255)
             gcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c                                                                                                 <  1:2286656/99‑1 (MQ=18)
             cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc                                                                                           >  1:1052733/1‑101 (MQ=255)
              gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.c                                                                                         >  1:4451631/1‑101 (MQ=255)
              gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.c                                                                                         >  1:929201/1‑101 (MQ=255)
                gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAG.CG..cc                                                                                            <  1:4232038/96‑1 (MQ=17)
                gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.ccc                                                                                       <  1:4752675/101‑2 (MQ=255)
                gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.ccc                                                                                       <  2:270951/101‑2 (MQ=255)
                 ccT.CCCCAAGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.ccct                                                                                      <  2:3823789/100‑3 (MQ=255)
                 aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.ccct                                                                                      >  1:485248/1‑99 (MQ=255)
                 aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.ccct                                                                                      >  1:2140057/1‑99 (MQ=255)
                   t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.cccttc                                                                                    <  1:2181228/101‑5 (MQ=255)
                   t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.cccttc                                                                                    <  2:1645830/101‑5 (MQ=255)
                     cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..ccc.cccttcg                                                                                   <  1:2330065/101‑6 (MQ=255)
                       cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTT.GCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG..CCC.CC..tgc                                                                                   <  2:1204802/101‑3 (MQ=0)
                          gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCgcgggc                                                                             <  1:654632/101‑1 (MQ=255)
                            gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCC..GCCAGAGGGGGGGAG.CGC.CCC.CC..TTCgccggcgc                                                                           >  2:1464777/1‑100 (MQ=25)
                             atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGggcggc                                                                          >  1:2251467/1‑101 (MQ=255)
                             atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGggcggc                                                                          >  1:74012/1‑101 (MQ=255)
                               atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCcg                                                                        >  2:1375132/1‑101 (MQ=255)
                                tttGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCcgc                                                                       >  1:4305555/1‑101 (MQ=255)
                                 ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCcgcg                                                                      >  2:535086/1‑101 (MQ=255)
                                 ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCcgcg                                                                      >  2:508930/1‑101 (MQ=255)
                                 ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCcgcg                                                                      >  1:4922108/1‑101 (MQ=255)
                                   gAGGAAAGATGAA...GGGG..GGCTTTCATTTT.GCCGAAAATTTCGGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGAc                                                                    <  2:1789784/101‑1 (MQ=37)
                                    aggtaggaggAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGAcg                                                                   <  2:1582857/92‑1 (MQ=25)
                                       cAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..CCCACC..TTCGCGGGCGGCCGCGAcgcgc                                                                >  2:1280495/1‑101 (MQ=255)
                                        aaGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGAcgcgcg                                                               >  2:3936860/1‑101 (MQ=255)
                                           gagga...GGGG..GGCTTTCATTTT.GGCAAAAATGCCCCGGGGGGGTAA.TGGGGGC..GCCCGGGGGGGGGCG.CC..CCC.CCCTTTCGCGGGCGGCCGCGACGCGcgcc                                                             <  2:4236344/97‑1 (MQ=255)
                                             gAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                          >  1:3437886/1‑101 (MQ=255)
                                             gAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                          <  2:4682148/101‑1 (MQ=255)
                                             gAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                          >  2:3637792/1‑101 (MQ=255)
                                             gAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                          >  2:2560414/1‑101 (MQ=255)
                                             gAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                          >  2:2673307/1‑101 (MQ=255)
                                             gAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                          >  1:1613799/1‑101 (MQ=255)
                                              ta...GGGG..GGGCTTCATATTCCCCGAAAATAGCGTGGGGGGTAAA.CGGCGCA..GCCAGGGGGGGGCGGCCC..CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                          <  1:3078535/100‑1 (MQ=255)
                                              aa...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGGCGCGCGCCGCaa                                                         >  2:1923001/1‑101 (MQ=255)
                                              aa...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCGGCaa                                                         >  1:2565231/1‑101 (MQ=255)
                                              aa...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                         >  1:2076968/1‑101 (MQ=255)
                                              aa...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                         >  1:1143316/1‑101 (MQ=255)
                                               a...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                        <  1:2855191/101‑1 (MQ=255)
                                                   gggg..ggCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                       <  2:349745/101‑1 (MQ=255)
                                                   gggg..ggCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                       <  1:365365/101‑1 (MQ=255)
                                                   gggg..ggCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                       <  1:4806444/101‑1 (MQ=255)
                                                    ggg..ggCTTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                      >  1:2062827/1‑101 (MQ=255)
                                                            tttCATCTT.GCCGAAAATATCCCGGGGGCTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                >  2:398831/1‑101 (MQ=255)
                                                             ttCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                               <  2:4382807/101‑1 (MQ=255)
                                                              tCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                              >  2:5009555/1‑101 (MQ=255)
                                                              tCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                              >  2:1030871/1‑101 (MQ=255)
                                                              tCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                              >  1:3201183/1‑101 (MQ=255)
                                                                  cTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                          <  2:3394470/101‑1 (MQ=255)
                                                                  cTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                          >  1:4349614/1‑101 (MQ=255)
                                                                   tt.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                         <  2:2997216/101‑1 (MQ=255)
                                                                    tggCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                         <  2:329169/99‑1 (MQ=255)
                                                                       ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                      <  1:3846955/101‑1 (MQ=255)
                                                                       ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                      <  2:2778809/101‑1 (MQ=255)
                                                                          aaaaGATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                   <  2:4846047/101‑1 (MQ=255)
                                                                           aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                  >  1:4295762/1‑101 (MQ=255)
                                                                           aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                  <  2:308107/101‑1 (MQ=255)
                                                                           aaaTATCCCGGGGGGTGAA.TGGCGCA..GACAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                  <  2:2399495/101‑1 (MQ=255)
                                                                              tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCC.CC..TTCGCGGGCGCCCGCACCCCCGGCCGCAAAAGGCCGGGCCCCTTGTTCTTTCtt                                >  2:588032/1‑100 (MQ=12)
                                                                               atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                              >  2:2879455/1‑101 (MQ=255)
                                                                                  ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCCC..CCC.CC..TTCGCGGGCGGCCGCGACCCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                           <  2:817038/101‑1 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTGGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                      >  1:1550806/1‑101 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                      >  2:4470104/1‑101 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                      >  1:4000052/1‑101 (MQ=255)
                                                                                           gAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                  >  2:2122268/1‑101 (MQ=255)
                                                                                            aa.TGGCGCA..GCCAGGGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                 <  1:4564322/101‑1 (MQ=255)
                                                                                            aa.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                 <  2:2274322/101‑1 (MQ=255)
                                                                                             a.TGGCGAA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGGCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa                <  2:2491210/101‑1 (MQ=255)
                                                                                               tGGCGCA..GCGGGGGGGGGGCGG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa               <  1:579683/101‑1 (MQ=255)
                                                                                               tGGCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa               <  2:2415479/101‑1 (MQ=255)
                                                                                                ggCGCA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg              <  1:4976029/101‑1 (MQ=255)
                                                                                                   gcA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg           <  1:2936855/101‑1 (MQ=255)
                                                                                                    cA..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg          >  2:3420182/1‑101 (MQ=255)
                                                                                                     cccgcaaggggggggCTG.CCC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg           <  1:1813751/92‑1 (MQ=255)
                                                                                                     a..GCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc         >  2:4630846/1‑101 (MQ=255)
                                                                                                      gggcggagggggggCCG.CTC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg          <  2:88663/93‑1 (MQ=255)
                                                                                                        gCCAGAGGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGcc        >  2:3168218/1‑101 (MQ=255)
                                                                                                          cAGAGGGGGCCGGCCC..CCC.CC..TCCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc      <  2:2003903/101‑1 (MQ=255)
                                                                                                           agagGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa     >  2:4053254/1‑101 (MQ=255)
                                                                                                            gagGGGGGCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt    >  1:4947990/1‑101 (MQ=255)
                                                                                                              ggggggCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCg  <  1:4470423/101‑1 (MQ=255)
                                                                                                              ggggggCAG.CGC.CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCg  >  1:2089694/1‑101 (MQ=255)
                                                                                                              |                                                                                                         
CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..CCC.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCG  >  NC_014034/2249281‑2249480

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG