Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,348 A→G 5.9% intergenic (‑144/‑161) 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,3480A→G 5.9% 19.4 ‑1.5 ‑0.6185intergenic (‑144/‑161)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (8/3):  ref base (90/84):  total (98/87)
Fisher's exact test strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.62e-01

Only 100 of 206 total aligned reads displayed.

GCTTGGGCACCCGGGGCAAGAAGGGGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.C.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATT  >  NC_014034/2249249‑2249447
                                                                                                          |                                                                                                          
gCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAC...GGGG..GGCTTTCATCTTGCCGa                                                                                                            <  1:2398280/101‑1 (MQ=255)
gCTTGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGa                                                                                                            >  2:329499/1‑101 (MQ=255)
 tttGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                                           <  1:3389755/100‑1 (MQ=255)
  ttGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaa                                                                                                          <  2:641662/101‑1 (MQ=255)
  ttGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaa                                                                                                          <  2:472504/101‑1 (MQ=255)
  ttGGGCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaa                                                                                                          >  2:3389510/1‑101 (MQ=255)
    gggCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                                        <  1:1160251/101‑1 (MQ=255)
    gggCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                                        >  1:1706438/1‑101 (MQ=255)
      gCACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAatat                                                                                                      <  2:458970/101‑1 (MQ=255)
       cACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATc                                                                                                     <  1:1510934/101‑1 (MQ=255)
       cACCCGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATc                                                                                                     <  2:3329254/101‑1 (MQ=255)
         cccGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATccc                                                                                                   <  1:87582/101‑1 (MQ=255)
         cccGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATccc                                                                                                   <  2:1371457/101‑1 (MQ=255)
          ccGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.g                                                                                                 <  1:472521/101‑1 (MQ=255)
          ccGGGGCAAGAAGGGGGCCCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.g                                                                                                 <  2:1997830/101‑1 (MQ=255)
           cGGGGCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.gg                                                                                                >  2:2229856/1‑101 (MQ=255)
            ggggCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.ggg                                                                                               >  1:1677672/1‑101 (MQ=255)
               gCAAGAAGGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.gggggg                                                                                            <  2:1527196/101‑1 (MQ=255)
                 aagaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTg                                                                                          <  2:3511322/101‑1 (MQ=255)
                   gaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGaa                                                                                        <  2:911974/101‑1 (MQ=255)
                   gaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGaa                                                                                        >  2:3752254/1‑101 (MQ=255)
                     ggggggCTCTGCCCCCGCCCCGTGCCGGGCCT.CTGCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.Tg                                                                                     <  1:3633811/100‑1 (MQ=37)
                     agGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.Tg                                                                                     <  1:2616160/101‑1 (MQ=255)
                      gggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.Tgg                                                                                    <  1:119679/101‑1 (MQ=255)
                        gggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGgcg                                                                                  >  1:2206728/1‑101 (MQ=255)
                        gggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGgcg                                                                                  <  2:1312539/101‑1 (MQ=255)
                        gggCTCTGCCCCCGTCCCGTGACGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGgcg                                                                                  >  2:1753689/1‑101 (MQ=255)
                           ctctGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..g                                                                             <  1:106933/101‑1 (MQ=255)
                               gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCag                                                                         <  2:1149245/99‑1 (MQ=255)
                               gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..GCCag                                                                         >  1:2243123/3‑101 (MQ=255)
                                 ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCagag                                                                       >  2:2916409/1‑101 (MQ=255)
                                 ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCagag                                                                       >  1:2703894/1‑101 (MQ=255)
                                 ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..CCCagag                                                                       >  2:1681540/1‑101 (MQ=255)
                                  cccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                      >  1:4129767/1‑101 (MQ=255)
                                    cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                    >  1:833390/1‑101 (MQ=255)
                                    cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                    >  2:3683041/1‑101 (MQ=255)
                                        ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                <  2:436042/101‑1 (MQ=255)
                                        ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                <  1:560713/101‑1 (MQ=255)
                                         cgtTCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                               <  1:1813397/101‑1 (MQ=255)
                                         cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                               >  2:4102191/1‑101 (MQ=255)
                                         cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                   >  2:2990221/8‑101 (MQ=18)
                                          gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                             <  2:1711979/101‑1 (MQ=255)
                                          gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                             <  1:1001720/101‑1 (MQ=255)
                                          gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                  >  2:4122634/7‑101 (MQ=21)
                                           tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                            <  1:2892862/101‑1 (MQ=255)
                                           caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAggggggg                                                                 >  2:882955/6‑100 (MQ=18)
                                            gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                                          >  1:3537406/1‑101 (MQ=255)
                                              cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.cc                                                       >  1:1918090/1‑101 (MQ=255)
                                              cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.cc                                                       >  2:2919430/1‑101 (MQ=255)
                                              cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.cc                                                       >  1:2850463/1‑101 (MQ=255)
                                               gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.ccc                                                      >  1:1616232/1‑101 (MQ=255)
                                                gggctccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.cg.c                                                          <  1:3059470/95‑1 (MQ=17)
                                                     ccccccGGGTTTTTTGGGGGAAGAGGAA...GGGG..GGTTTTCATTTTGCCAAAAAAATCCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.ccccctt                                                  <  1:1422384/100‑4 (MQ=21)
                                                       ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTcgc                                                <  2:1633285/101‑1 (MQ=255)
                                                       ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTcgc                                                <  1:811842/101‑1 (MQ=255)
                                                       ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTcgc                                                >  1:2354331/1‑101 (MQ=255)
                                                        cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTcgcg                                               <  2:3190971/101‑1 (MQ=255)
                                                          cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGAGAAT.TGGCCAAAGGCCAAGAGGGGGCAG.CG.C.CCCCtgacc                                                >  2:1879142/1‑97 (MQ=2)
                                                          cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC..AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.C.CCCCtgccc                                                >  2:1663196/1‑97 (MQ=0)
                                                          cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCgcggg                                             <  2:3162182/101‑1 (MQ=255)
                                                          cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCgcggg                                             >  1:433227/1‑101 (MQ=255)
                                                          cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCgcggg                                             >  1:636164/1‑101 (MQ=255)
                                                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCgcgggc                                            >  1:1763632/1‑101 (MQ=255)
                                                            ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCgcgtgcg                                           >  1:1361193/1‑101 (MQ=255)
                                                            ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCgcgggcg                                           <  2:2533023/101‑1 (MQ=255)
                                                            gaatatTGGAAGAAAGATGAG...GGGG..TTTTTTCCTCGTGCAAAAAATCTCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCCC.C.CCCCTTCgcgggcg                                           <  1:762483/99‑1 (MQ=255)
                                                                atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCcg                                       >  1:3698587/1‑101 (MQ=255)
                                                                atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTCCGCGGGCGGCcg                                       >  1:991321/1‑101 (MQ=255)
                                                                 tttGAGGCAAGATGAA...GGGG..GCCTTTCATCTTGCCGAAAATACCCC.GGGGGGTGAA.TGGCGCA..GCCAAAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCcgc                                      <  1:1258836/101‑1 (MQ=37)
                                                                  ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGGGGGGGGCAG.CGCC.CCCCTTCGCGGCCGGCcgcg                                     >  2:1268224/1‑101 (MQ=37)
                                                                    gAGGCAAGATGAA...GGGG..GTCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGAc                                   <  2:1972989/101‑1 (MQ=255)
                                                                     aGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGAcg                                  >  2:3897094/1‑101 (MQ=255)
                                                                      ggCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGAcgc                                 >  2:370932/1‑101 (MQ=255)
                                                                       gCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGAcgcg                                <  2:1112260/101‑1 (MQ=255)
                                                                        cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAGAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGAcgcgc                               >  1:3765924/1‑101 (MQ=255)
                                                                         aaGATGAG...GGGG..GGCTTTCATCTTGCCGAAAATAGCGC.GGGGGGTGAA.TGGCGCA..GCGAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGAcgcgcg                              <  2:1066977/101‑1 (MQ=37)
                                                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGAcgcgcg                              >  2:1837501/1‑101 (MQ=255)
                                                                         aaGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGAcgcgcg                              <  1:2440432/101‑1 (MQ=255)
                                                                          aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACgcgcgc                             >  1:4232831/1‑101 (MQ=255)
                                                                          aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACgcgcgc                             >  2:1588134/1‑101 (MQ=255)
                                                                            aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGcgccg                           >  2:3618452/1‑101 (MQ=255)
                                                                             tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGcgccgc                          >  2:3799270/1‑101 (MQ=255)
                                                                              gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                         <  1:1461113/101‑1 (MQ=255)
                                                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                       <  2:2262820/101‑1 (MQ=255)
                                                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                       <  1:1899413/101‑1 (MQ=255)
                                                                                   ggggg..ggTTTTCTTCTTGCGAAAAATTCCGG.GGGGGGTAAA.TGGCGCG..GCCGGGGGGGGGCGGCCC.C.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                       <  1:3132242/100‑1 (MQ=255)
                                                                                    gggg..ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGAGAA..GCCAGAGGGGGGCAG.CG.CACGCCTTCGCGGGCGGGCGGGAAGCGcgccgccaaa                      >  2:4008941/1‑97 (MQ=25)
                                                                                     ggg..ggCTTTCATCTTGCCGAAAATATCCG.GGGGGGTGAA.TGGCTCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                     <  1:382674/101‑1 (MQ=37)
                                                                                     ggg..ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                     <  2:4223764/101‑1 (MQ=255)
                                                                                          ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                  >  2:2568943/1‑101 (MQ=255)
                                                                                           gCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                 <  1:3207496/101‑1 (MQ=255)
                                                                                           gCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                 <  1:3527858/101‑1 (MQ=255)
                                                                                               tCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTTGTGcc             >  2:2107527/1‑101 (MQ=255)
                                                                                                cATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc            <  1:2527817/101‑1 (MQ=255)
                                                                                                  tCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt          <  2:2302495/101‑1 (MQ=255)
                                                                                                   cTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGCCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt         >  1:2265057/1‑101 (MQ=255)
                                                                                                     tGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg       >  2:2885210/1‑101 (MQ=255)
                                                                                                      gCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt      >  2:3357811/1‑101 (MQ=255)
                                                                                                       ccGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc     >  1:3153813/1‑101 (MQ=255)
                                                                                                        cGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa    <  2:1009263/101‑1 (MQ=255)
                                                                                                          |                                                                                                          
GCTTGGGCACCCGGGGCAAGAAGGGGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.C.CCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATT  >  NC_014034/2249249‑2249447

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG