Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,330 G→A 3.2% intergenic (‑126/‑179) 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,3300G→A 3.2% 7.6 ‑0.4 ‑0.0187intergenic (‑126/‑179)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (93/87):  total (97/90)
Fisher's exact test strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.77e-01

Only 100 of 207 total aligned reads displayed.

AGATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAG  >  NC_014034/2249231‑2249428
                                                                                                        |                                                                                                            
agaTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...ggg                                                                                                              >  1:472561/1‑101 (MQ=255)
 gaTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg                                                                                                             >  2:2918460/1‑101 (MQ=255)
       caGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCttt                                                                                                     >  2:2944183/1‑101 (MQ=255)
        aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTc                                                                                                    <  2:2682658/101‑1 (MQ=255)
         gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCa                                                                                                   <  2:3723645/101‑1 (MQ=255)
         gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCa                                                                                                   >  1:136121/1‑101 (MQ=255)
           tACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATc                                                                                                 <  1:629049/101‑1 (MQ=255)
            aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                                >  2:484617/1‑101 (MQ=255)
            aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                                <  2:3361288/101‑1 (MQ=255)
            aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                                <  1:3145267/101‑1 (MQ=255)
                  gCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAC...GGGG..GGCTTTCATCTTGCCGa                                                                                          <  1:2398280/101‑1 (MQ=255)
                  gCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGa                                                                                          >  2:329499/1‑101 (MQ=255)
                   cTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                         <  1:2138596/101‑1 (MQ=255)
                    ttGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaa                                                                                        <  2:641662/101‑1 (MQ=255)
                    ttGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaa                                                                                        >  2:3389510/1‑101 (MQ=255)
                     tGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                       >  2:2259085/1‑101 (MQ=255)
                      gggCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                      >  1:1706438/1‑101 (MQ=255)
                      gggCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                      <  1:1160251/101‑1 (MQ=255)
                        gCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAatat                                                                                    <  2:1046824/101‑1 (MQ=255)
                         cACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATc                                                                                   <  1:4028292/101‑1 (MQ=255)
                         cACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATc                                                                                   <  2:2165523/101‑1 (MQ=255)
                         cACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATc                                                                                   >  2:3660498/1‑101 (MQ=255)
                         cACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATc                                                                                   <  2:3329254/101‑1 (MQ=255)
                         cACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATc                                                                                   <  2:2549880/101‑1 (MQ=255)
                           cccGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGTA...GGGG..GGCTTTCATCTTGCCGAAAATATccc                                                                                 <  1:3762461/101‑1 (MQ=255)
                           cccGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATccc                                                                                 <  2:1371457/101‑1 (MQ=255)
                            ccGGGGCAAGAAGGGGGCCCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.g                                                                               <  2:1997830/101‑1 (MQ=255)
                              ggggCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.ggg                                                                             >  1:1677672/1‑101 (MQ=255)
                                 gCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.gggggg                                                                          <  2:1527196/101‑1 (MQ=255)
                                  cAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.ggggggg                                                                         >  1:24659/1‑100 (MQ=255)
                                     gaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGaa                                                                      <  2:911974/101‑1 (MQ=255)
                                       ggggggCTCTGCCCCCGCCCCGTGCCGGGCCT.CTGCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.Tg                                                                   <  1:3633811/100‑1 (MQ=37)
                                       agGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.Tg                                                                   >  2:1251327/1‑101 (MQ=255)
                                        gggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.Tgg                                                                  <  1:119679/101‑1 (MQ=255)
                                        gggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCAC.GGGGGGTGAA.Tgg                                                                  >  1:450802/1‑101 (MQ=255)
                                         ggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGgc                                                                 >  2:368438/1‑101 (MQ=255)
                                          gggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGgc.g                                                               <  1:2264758/101‑1 (MQ=255)
                                          gggCTCTGCCCCCGTCCCGTGACGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGgc.g                                                               >  2:1753689/1‑101 (MQ=255)
                                             ctctGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..g                                                          <  1:106933/101‑1 (MQ=255)
                                              tctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..Gc                                                         >  1:4224026/6‑101 (MQ=255)
                                                tgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCa                                                       >  1:3225845/4‑101 (MQ=255)
                                                 gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCag                                                      <  2:1149245/99‑1 (MQ=255)
                                                 gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCag                                                      <  1:2604459/99‑1 (MQ=255)
                                                 gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGC.GCA..GCCag                                                      >  1:2243123/3‑101 (MQ=255)
                                                   ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCagag                                                    >  1:2703894/1‑101 (MQ=255)
                                                   ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..CCCagag                                                    >  2:1681540/1‑101 (MQ=255)
                                                    cccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAg.g                                                  <  1:368489/101‑1 (MQ=255)
                                                     ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAg.gg                                                 <  2:2269323/101‑1 (MQ=255)
                                                         ccgatctcGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGc                                             <  2:425147/94‑1 (MQ=255)
                                                         cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGc                                             >  2:3872364/1‑101 (MQ=255)
                                                          ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCa                                            <  1:560713/101‑1 (MQ=255)
                                                           cgtTCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAg                                           <  1:1813397/101‑1 (MQ=255)
                                                           cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAg                                           <  2:3415444/101‑1 (MQ=255)
                                                           cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAg                                           >  1:84139/1‑101 (MQ=255)
                                                           cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGC.GAA..GCCAGAg.gggg                                               >  2:2990221/8‑101 (MQ=18)
                                                            gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAg.c                                         <  2:1711979/101‑1 (MQ=255)
                                                            gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAg.c                                         >  2:4053323/1‑101 (MQ=255)
                                                            gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAg.c                                         <  1:1001720/101‑1 (MQ=255)
                                                             tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAg.cg                                        <  1:2892862/101‑1 (MQ=255)
                                                              gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAg.cg.c                                      >  1:3537406/1‑101 (MQ=255)
                                                                cGGGACT.CCTCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.ccc                                    <  2:459552/101‑1 (MQ=255)
                                                                cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAgatcgg                                      >  1:425204/1‑96 (MQ=255)
                                                                cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.ccc                                    >  1:2850463/1‑101 (MQ=255)
                                                                cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.ccc                                    >  1:995632/1‑101 (MQ=255)
                                                                  gggctccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGC.GAA..GCCAGAG.GGGGGCAg.cg.c                                      <  1:3059470/95‑1 (MQ=17)
                                                                  ggACT.CCCCCGGGATATTTGAGGCAAGATGAC...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.ccccc                                  >  1:3808942/1‑101 (MQ=255)
                                                                   gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.cccccc                                 <  2:3593020/101‑2 (MQ=255)
                                                                    aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.TGGGGCAG.CG.cccccct                                <  2:4129505/101‑3 (MQ=255)
                                                                       ccccccGGGTTTTTTGGGGGAAGAGGAA...GGGG..GGTTTTCATTTTGCCAAAAAAATCCCGGGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.cccccctt                               <  1:1422384/100‑4 (MQ=21)
                                                                       ccccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.ccccccttc                              <  2:2892488/100‑5 (MQ=255)
                                                                        cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.ccccccttcg                             <  2:3959378/101‑6 (MQ=255)
                                                                        cccccGGGATATTTGAGGCAAGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.ccccccttcg                             <  2:450745/101‑6 (MQ=255)
                                                                         ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTcgc                             <  1:1547977/101‑1 (MQ=255)
                                                                         ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTcgc                             <  2:1633285/101‑1 (MQ=255)
                                                                         ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTcgc                             <  1:811842/101‑1 (MQ=255)
                                                                         ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTcgc                             >  1:2354331/1‑101 (MQ=255)
                                                                          cccGGGATATTTGCAGCAAGGTGAAACCGGTT..GCCTTTCATCTGGCCCGAAATATCCC.GGGGGGTGAA.TGGC.GAA..GCCAGAG.GGGGGCAG.CG.ccctcctgc                              <  1:291216/101‑7 (MQ=14)
                                                                          cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTcgcg                            <  2:3190971/101‑1 (MQ=255)
                                                                          cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAAGATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTcgcg                            >  1:3579042/1‑101 (MQ=255)
                                                                           ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCgcgg                           >  1:82437/1‑101 (MQ=255)
                                                                            cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGAGAAT.TGGC.CAAAGGCCAAGA.GGGGGCAG.CG.CCCCCtgacc                             >  2:1879142/1‑97 (MQ=2)
                                                                             gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCgcgggc                         <  2:410638/101‑1 (MQ=255)
                                                                             gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCgcgggc                         >  1:1763632/1‑101 (MQ=255)
                                                                             gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCgcgggc                         >  2:3783320/1‑101 (MQ=255)
                                                                              ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCgcgggcg                        <  2:2533023/101‑1 (MQ=255)
                                                                                 tatTTGCGGAAAGAGGAA...GGGGCATGTTTTCATCTTGCCCAAAATATCCC..GGGGGTGAATTGGCAAAA..GGCAAAGAGGGGGCAG.CG.CCCCCtgccccgccg                        <  2:3139117/101‑10 (MQ=9)
                                                                                   tttGAGGCAAGATGAA...GGGG..GCCTTTCATCTTGCCGAAAATACCCC.GGGGGGTGAA.TGGC.GCA..GCCAAAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCcgc                   <  1:1258836/101‑1 (MQ=37)
                                                                                      gAGGCAAGATGAA...GGGG..GTCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAc                <  2:1972989/101‑1 (MQ=255)
                                                                                        ggCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgc              >  1:910587/1‑101 (MQ=255)
                                                                                         gCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcg             <  2:1112260/101‑1 (MQ=255)
                                                                                           aaGATGAG...GGGG..GGCTTTCATCTTGCCGAAAATAGCGC.GGGGGGTGAA.TGGC.GCA..GCGAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg           <  2:1066977/101‑1 (MQ=37)
                                                                                           aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg           >  2:1837501/1‑101 (MQ=255)
                                                                                            aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACgcgcgc          >  2:1588134/1‑101 (MQ=255)
                                                                                                gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa      <  1:1461113/101‑1 (MQ=255)
                                                                                                gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa      <  2:3521903/101‑1 (MQ=255)
                                                                                                gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa      >  1:735256/1‑101 (MQ=255)
                                                                                                  a...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc    <  1:1899413/101‑1 (MQ=255)
                                                                                                     ggggg..ggTTTTCTTCTTGCGAAAAATTCCGG.GGGGGGTAAA.TGGC.GCG..GCCGGGG.GGGGGCGGCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc    <  1:3132242/100‑1 (MQ=255)
                                                                                                      gggg..ggTTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa   <  2:899651/101‑1 (MQ=255)
                                                                                                      gggg..ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa   >  2:1269296/1‑101 (MQ=255)
                                                                                                        |                                                                                                            
AGATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAG  >  NC_014034/2249231‑2249428

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