Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,327 +C 2.0% intergenic (‑123/‑182) 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,3273.→C 2.0% 12.8 ‑0.3 ‑0.1357intergenic (‑123/‑182)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (3/4):  ref base (177/173):  total (180/177)
Fisher's exact test strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.21e-01

Only 100 of 396 total aligned reads displayed.

TTGAGATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTG.CCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAC  >  NC_014034/2249228‑2249426
                                                                                                       |                                                                                                            
ttGAGATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGaa                                                                                                               >  1:3201834/1‑101 (MQ=255)
 tGAGATGCACAGATACCGCAGCCTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...g                                                                                                           <  2:2056521/101‑1 (MQ=255)
  gagaTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gg                                                                                                          <  2:4117691/101‑1 (MQ=255)
    gaTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg                                                                                                        >  1:5995169/1‑101 (MQ=255)
      tGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg..gg                                                                                                    >  2:5094932/1‑101 (MQ=255)
        cacaGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCt                                                                                                  <  2:6931060/101‑1 (MQ=255)
        cacaGATACCGCAGCTGGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCt                                                                                                  <  2:6541101/101‑1 (MQ=255)
         acaGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCtt                                                                                                 <  2:3055129/101‑1 (MQ=255)
          caGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCttt                                                                                                <  1:5867733/101‑1 (MQ=255)
          caGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCttt                                                                                                <  1:672067/101‑1 (MQ=255)
            gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCa                                                                                              >  1:1524407/1‑101 (MQ=255)
             aTACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCAt                                                                                             <  1:4117337/101‑1 (MQ=255)
                ccGCAGCTTGGGCCCCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCtt                                                                                          <  2:3509808/101‑1 (MQ=255)
                ccGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCtt                                                                                          <  1:2198946/101‑1 (MQ=255)
                    aGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCg                                                                                      >  1:2310053/1‑101 (MQ=255)
                      cgtGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                    <  2:735009/99‑1 (MQ=255)
                      cTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                    >  2:2951317/1‑101 (MQ=255)
                      cTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                    <  1:1312018/101‑1 (MQ=255)
                        tGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                  >  2:1649286/1‑101 (MQ=255)
                        tGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                  >  2:3886469/1‑101 (MQ=255)
                        tGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTAATCTTGCCGagaa                                                                                  >  2:3300747/1‑101 (MQ=38)
                         gggCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                 <  1:1715073/101‑1 (MQ=255)
                         gggCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                 >  1:4807769/1‑101 (MQ=255)
                         gggCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                 >  2:5280081/1‑101 (MQ=255)
                         gcgCACCCGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGCATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                 <  1:3886632/99‑1 (MQ=255)
                              cccGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATccc                                                                            >  2:859900/1‑101 (MQ=255)
                               ccGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCg                                                                           >  1:300406/1‑101 (MQ=255)
                               ccGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCg                                                                           <  2:4546550/101‑1 (MQ=255)
                                cGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                          <  1:6911051/101‑1 (MQ=255)
                                cGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                          <  2:461277/101‑1 (MQ=255)
                                cGGGGCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                          <  2:561713/101‑1 (MQ=255)
                                   ggCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCggggg                                                                       >  1:6303693/1‑101 (MQ=255)
                                   ggCAAGAAGGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATAGTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCggggg                                                                       >  1:4204837/1‑101 (MQ=255)
                                    gCAAGAAGGGGGCTCTG.CCCCTGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                      <  2:1181186/101‑1 (MQ=255)
                                      aagaagGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTg                                                                    >  2:2169379/1‑101 (MQ=255)
                                       agaagGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGa                                                                   >  2:2681487/1‑101 (MQ=255)
                                        gaagGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGaa                                                                  <  1:3082853/101‑1 (MQ=255)
                                         gagGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAt                                                                 <  1:641727/100‑1 (MQ=255)
                                         aagGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAt                                                                 >  1:5555789/1‑101 (MQ=255)
                                         aagGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAt                                                                 <  2:1537081/101‑1 (MQ=255)
                                         aagGGGGCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAt                                                                 <  2:6189599/101‑1 (MQ=255)
                                            agggCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGgc                                                              <  1:5466088/100‑1 (MQ=255)
                                             gggCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGgcg                                                             <  1:1009599/101‑1 (MQ=255)
                                              ggCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGTGGTGAATGgcgc                                                            <  1:5131447/101‑1 (MQ=255)
                                              ggCTCTGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGgcgc                                                            >  2:6075470/1‑101 (MQ=255)
                                                ctctGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...g                                                       >  2:440226/1‑101 (MQ=255)
                                                ctctGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...g                                                       <  1:3419422/101‑1 (MQ=255)
                                                ctctGCCCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...g                                                       >  1:4492496/1‑101 (MQ=255)
                                                 tctgccccc.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAATGGCGCA...Gc                                                      >  1:4337646/6‑101 (MQ=255)
                                                  ctgccccc.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...Gcc                                                     <  1:84478/97‑1 (MQ=255)
                                                    gccccc.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCag                                                   >  2:5410155/3‑101 (MQ=255)
                                                    gccccc.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCag                                                   >  1:1891169/3‑101 (MQ=255)
                                                    gacccc.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCag                                                   >  2:2783494/3‑101 (MQ=255)
                                                      cccc.CTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCagag                                                 <  2:2304837/101‑1 (MQ=255)
                                                           gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAg.gggg                                            >  1:1163344/1‑101 (MQ=255)
                                                             cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGc                                          <  1:2169459/101‑1 (MQ=255)
                                                             cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGc                                          >  1:2832057/1‑101 (MQ=255)
                                                             cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGc                                          >  1:1399096/1‑101 (MQ=255)
                                                               cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCGTCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAg                                        >  1:3616988/1‑101 (MQ=255)
                                                               cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAg                                        <  2:6706884/101‑1 (MQ=255)
                                                               cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAAGGCGCA...GCCAGAG.GGGGGCAg                                        >  1:497490/1‑101 (MQ=255)
                                                                gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAgc                                       >  2:6251473/1‑101 (MQ=255)
                                                                 tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAgcg                                      >  1:3719900/1‑101 (MQ=255)
                                                                 caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAA...GCCAGAG.GGGGGc                                          >  1:1071622/6‑101 (MQ=18)
                                                                 caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAAGATCCCGGGGGGTGAATGGCGAA...GCCAGAG.GGGGGc                                          >  1:3373146/6‑101 (MQ=18)
                                                                  tccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAgcg.c                                    <  2:1031953/100‑1 (MQ=255)
                                                                  gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAgcg.c                                    >  2:4909331/1‑101 (MQ=255)
                                                                  gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAATGGAGAA...GCCAGAG.GGGGGCAgcg.c                                    >  1:5228563/1‑101 (MQ=255)
                                                                   ccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.cc                                   <  2:1744695/101‑1 (MQ=255)
                                                                     gggGCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.cccc                                 <  2:2636775/101‑1 (MQ=255)
                                                                     gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.cccc                                 <  2:5556588/101‑1 (MQ=255)
                                                                     gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.cccc                                 >  2:1144738/1‑101 (MQ=255)
                                                                      ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.ccccc                                <  2:1092191/101‑1 (MQ=255)
                                                                       gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAA...GCCAGAG.GGGGGCAGCG.cc                                   <  1:202820/96‑1 (MQ=17)
                                                                       gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.cccccc                               <  1:294354/101‑2 (MQ=255)
                                                                        aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.cccccct                              <  2:68211/101‑3 (MQ=255)
                                                                         cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATTGCGCA...GCCAGAG.GGGGGCAGCG.cccccctt                             <  2:1314294/101‑4 (MQ=255)
                                                                         cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.cccccctt                             <  2:6863426/101‑4 (MQ=255)
                                                                         cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.cccccctt                             <  2:1236503/101‑4 (MQ=255)
                                                                          t.CCCCCGGGATATCTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.ccccccttc                            <  2:4337424/101‑5 (MQ=255)
                                                                            cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.ccccccttcg                           >  2:3151444/1‑96 (MQ=255)
                                                                            cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.ccccccttcg                           >  1:4684026/1‑96 (MQ=255)
                                                                            cccccGGAATATTTGAGGCAAGATGAA...GGGG..GGCTTTCAGCTGGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.ccccccttcg                           <  1:6161009/101‑6 (MQ=25)
                                                                               ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCgcgg                         <  1:3788549/101‑1 (MQ=255)
                                                                               ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCgcgg                         >  2:5860343/1‑101 (MQ=255)
                                                                                cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCA.AGAGGGGGCAGCG.CCCCCTTccc                           <  1:6832091/101‑3 (MQ=0)
                                                                                 gggATATTTGAGCAAAGTTGAA...GGGG..GGCTTTAATCTTGCCGAAAAAATCCCGGGGGGGGAATGGGGCA...GCCAGAG.GGGGGCGGCGCCCCCCTTCgcgggc                       >  2:3741255/1‑101 (MQ=255)
                                                                                  ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCgcgggcg                      >  2:4029409/1‑101 (MQ=255)
                                                                                  ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCgcgggcg                      >  2:191590/1‑101 (MQ=255)
                                                                                    atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGggcggc                    >  1:6078253/1‑101 (MQ=255)
                                                                                     tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGcc                   <  2:4227588/101‑1 (MQ=255)
                                                                                      atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCcg                  <  1:2418644/101‑1 (MQ=255)
                                                                                      atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCcg                  >  1:533331/1‑101 (MQ=255)
                                                                                       tttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCTCGGGGGGGGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCcgc                 >  1:7068660/1‑101 (MQ=37)
                                                                                           aGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCCGCGAcg             <  2:547596/101‑1 (MQ=255)
                                                                                           aGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCCGCGAcg             >  2:5770225/1‑101 (MQ=255)
                                                                                                aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCCGCGACgcgcgc        >  1:664077/1‑101 (MQ=255)
                                                                                                  aTGAA...GGGG..TGCGTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg      <  2:1376315/101‑1 (MQ=25)
                                                                                                  aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg      >  2:776742/1‑101 (MQ=255)
                                                                                                      a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc  >  2:2892132/1‑101 (MQ=255)
                                                                                                       |                                                                                                            
TTGAGATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTG.CCCC.GTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCA...GCCAGAG.GGGGGCAGCG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAC  >  NC_014034/2249228‑2249426

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG