Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,373 C→A 6.8% intergenic (‑169/‑136) 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,3730C→A 6.8% 31.2 ‑0.5 ‑1.9234intergenic (‑169/‑136)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (8/8):  ref base (123/95):  total (131/103)
Fisher's exact test strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.36e-02

Only 100 of 261 total aligned reads displayed.

GGCTCTG.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGG  >  NC_014034/2249274‑2249472
                                                                                                              |                                                                                                         
  ctctGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..g                                                                                                       <  1:4295982/101‑1 (MQ=255)
   tctgcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gc                                                                                                      <  2:2262889/96‑1 (MQ=255)
   tctgcccc.cGTCCCGTGCCGCGACT.CTCCCGGGATATATGAAGCAAGAAGAA...GGGG..GCC.TTACATCTT.GCCGAAAATATCCCGGGGCGTGAA.TGGCGCA..Gc                                                                                                      <  2:5020550/96‑1 (MQ=25)
    ctgcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..Gcc                                                                                                     >  1:3458951/5‑101 (MQ=255)
    ctG.CCCTCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                                                                     <  2:2212381/101‑1 (MQ=255)
      gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                   >  1:3994163/3‑101 (MQ=255)
      gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGAAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATACCGGGGGGGGAA.GGGCGCA..GCCag                                                                                                   >  1:2491384/3‑101 (MQ=38)
        ccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                 >  1:817140/1‑101 (MQ=255)
            cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                              >  1:2274477/1‑101 (MQ=255)
            cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                              >  1:231295/1‑101 (MQ=255)
             gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                             >  1:3213276/1‑101 (MQ=255)
             gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGAGAA.TGGCGCA..GCCAGAggggg                                                                                             >  2:1813618/1‑101 (MQ=255)
              tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                            <  1:392326/101‑1 (MQ=255)
               cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                           >  1:2415653/1‑101 (MQ=255)
                ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                          <  2:608220/101‑1 (MQ=255)
                ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GACAGAGGGGGGCa                                                                                          <  1:2728576/101‑1 (MQ=255)
                ccgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAAC.CGGA..TGCATTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggg                                                                                              >  1:4872835/9‑101 (MQ=21)
                 cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                         <  1:3904232/101‑1 (MQ=255)
                 cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                         <  2:4307761/101‑1 (MQ=255)
                 cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                             >  2:2381758/8‑101 (MQ=18)
                 cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                             >  2:1685360/8‑101 (MQ=21)
                   caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                                           >  2:1143867/6‑101 (MQ=18)
                   caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                                           <  2:902873/96‑1 (MQ=18)
                    gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                                                                    >  2:2661997/1‑101 (MQ=255)
                      cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc                                                                                  <  1:2591059/101‑1 (MQ=255)
                      cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc                                                                                  >  1:1052733/1‑101 (MQ=255)
                        ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCC.TTTCATCTT.GCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.cg.c                                                                                    <  2:3846146/101‑1 (MQ=17)
                         gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccc                                                                               <  1:2807368/101‑2 (MQ=255)
                         gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccc                                                                               <  1:4752675/101‑2 (MQ=255)
                          ccT.CCCCAAGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccct                                                                              <  2:3823789/100‑3 (MQ=255)
                           cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccctt                                                                             <  1:4041099/101‑4 (MQ=255)
                           cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGAGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccctt                                                                             >  2:1626094/1‑98 (MQ=255)
                            t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                                                            <  2:1645830/101‑5 (MQ=255)
                              cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttcg                                                                           <  1:2215211/101‑6 (MQ=255)
                               ccccGGGATATTTGCGGCAAGATGAA...GGGGCATGA.TTTCATCTT.GCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCC..tg                                                                            >  1:3435718/1‑100 (MQ=0)
                                cccGGGATATTTGCGGCAAGATGAA...GGGGCATGA.TTTCATCTT.GCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCC..tgc                                                                           <  2:1204802/101‑3 (MQ=0)
                                   gggATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCgcgggc                                                                     <  1:654632/101‑1 (MQ=255)
                                   gggATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCgcgggc                                                                     >  1:568043/1‑101 (MQ=255)
                                   gggATATTTGAGGAAAGATGAG...GGGG..GTT.TTTCATCTT.GCCGAAAATATCCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCgcgggc                                                                     <  1:3515391/101‑1 (MQ=21)
                                     gATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCC..GCCAGAGGGGGGGAG.CGCCCCCC..TTCgccggcgc                                                                   >  2:1464777/1‑100 (MQ=25)
                                     gATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCgcgggcgg                                                                   >  1:4569746/1‑101 (MQ=255)
                                      atatTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGggcggc                                                                  >  1:74012/1‑101 (MQ=255)
                                       tatTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGcc                                                                 <  2:2757136/101‑1 (MQ=255)
                                        atTTGCGGCAAGATGAA...GGGGCATGA.TTTCATCTT.GCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCC..tgccccgccgg                                                                   >  2:4122224/1‑91 (MQ=2)
                                        atTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCcg                                                                >  2:1375132/1‑101 (MQ=255)
                                         tttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCcgc                                                               <  2:4880746/101‑1 (MQ=255)
                                         tttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCcgc                                                               >  1:4305555/1‑101 (MQ=255)
                                          ttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCcgcg                                                              >  1:4922108/1‑101 (MQ=255)
                                          ttGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCcgcg                                                              >  2:1926736/1‑101 (MQ=255)
                                                    gagga...GGGG..GGC.TTTCATTTT.GGCAAAAATGCCCCGGGGGGGTAA.TGGGGGC..GCCCGGGGGGGGGCG.CC.CCCCCCTTTCGCGGGCGGCCGCGACGCGcgcc                                                     <  2:4236344/97‑1 (MQ=255)
                                                    aTGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGcgccg                                                    <  2:2697914/101‑1 (MQ=255)
                                                     tGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGcgccgc                                                   >  2:1777897/1‑101 (MQ=255)
                                                      gAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                  <  2:540454/101‑1 (MQ=255)
                                                      gAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                  <  2:4682148/101‑1 (MQ=255)
                                                      gAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                  >  2:515019/1‑101 (MQ=255)
                                                       ta...GGGG..GGG.CTTCATATTCCCCGAAAATAGCGTGGGGGGTAAA.CGGCGCA..GCCAGGGGGGGGCGGCCC.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                  <  1:3078535/100‑1 (MQ=255)
                                                       aa...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                 >  2:1052694/1‑101 (MQ=255)
                                                       aa...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                 >  1:1143316/1‑101 (MQ=255)
                                                       aa...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                 >  1:2076968/1‑101 (MQ=255)
                                                       aa...GGGC..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                 >  2:3693541/1‑101 (MQ=255)
                                                        a...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                >  2:4074801/1‑101 (MQ=255)
                                                        a...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                <  2:2403842/101‑1 (MQ=255)
                                                            gggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                               <  2:349745/101‑1 (MQ=255)
                                                            gggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                               >  2:3529936/1‑101 (MQ=255)
                                                            gggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                               <  1:2124760/101‑1 (MQ=255)
                                                            gggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                               <  1:365365/101‑1 (MQ=255)
                                                             ggg..ggC.TTTCATCTT.GCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                              >  1:2062827/1‑101 (MQ=255)
                                                                       ttCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                       <  2:4382807/101‑1 (MQ=255)
                                                                        tCATCTT.GCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                      <  2:1296295/101‑1 (MQ=255)
                                                                        tCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                      <  2:2483346/101‑1 (MQ=255)
                                                                           tCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                   <  1:2859550/101‑1 (MQ=255)
                                                                            cTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                  >  1:4263215/1‑101 (MQ=255)
                                                                            cTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                  >  1:4349614/1‑101 (MQ=255)
                                                                              tggCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                 <  2:329169/99‑1 (MQ=255)
                                                                                 ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                              <  1:3846955/101‑1 (MQ=255)
                                                                                 ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                              <  2:2778809/101‑1 (MQ=255)
                                                                                 ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                              <  1:4510971/101‑1 (MQ=255)
                                                                                  ggAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAT.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                             <  1:2340286/100‑1 (MQ=255)
                                                                                  cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                             >  1:2150742/1‑101 (MQ=255)
                                                                                   gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                            >  1:1794654/1‑101 (MQ=255)
                                                                                    aaaaGATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                           <  2:4846047/101‑1 (MQ=255)
                                                                                     aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTCTGTCATTc                          >  2:3328181/1‑101 (MQ=255)
                                                                                             gggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                  <  2:1009809/100‑1 (MQ=255)
                                                                                             cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                  >  2:2894041/1‑101 (MQ=255)
                                                                                              ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                 <  2:2439222/101‑1 (MQ=255)
                                                                                               gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                >  1:702558/1‑101 (MQ=255)
                                                                                               gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                <  2:3985285/101‑1 (MQ=255)
                                                                                                 gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              <  2:3054663/101‑1 (MQ=255)
                                                                                                 gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              >  2:2503832/1‑101 (MQ=255)
                                                                                                 gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              >  2:4470104/1‑101 (MQ=255)
                                                                                                  ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TGCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg             <  1:4829377/101‑1 (MQ=255)
                                                                                                     gAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGACCCTTTGTCATTCTGTCGCACCCATGTCa          >  2:824026/1‑101 (MQ=255)
                                                                                                       a.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa        >  1:1413744/1‑101 (MQ=255)
                                                                                                         tGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa       >  1:4413198/1‑101 (MQ=255)
                                                                                                          ggCGCA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg      <  1:4976029/101‑1 (MQ=255)
                                                                                                            cgcA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt    <  1:1560074/101‑1 (MQ=255)
                                                                                                             gcA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg   <  1:2936855/101‑1 (MQ=255)
                                                                                                              cA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg  >  2:3420182/1‑101 (MQ=255)
                                                                                                              cA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg  <  2:3100241/101‑1 (MQ=255)
                                                                                                              cA..GCCAGAGGGGGGCAG.CGCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg  >  1:3135105/1‑101 (MQ=255)
                                                                                                              |                                                                                                         
GGCTCTG.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGC.TTTCATCTT.GCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGG  >  NC_014034/2249274‑2249472

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