Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,380 A→G 2.7% intergenic (‑176/‑129) 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,3800A→G 2.7% 3.9 ‑0.4 ‑0.0209intergenic (‑176/‑129)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (3/4):  ref base (114/88):  total (117/92)
Fisher's exact test strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

Only 100 of 241 total aligned reads displayed.

GCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATC  >  NC_014034/2249280‑2249479
                                                                                                                |                                                                                                          
cccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                                                                            >  2:4300532/2‑101 (MQ=255)
 ccgcgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                           >  1:3510041/4‑101 (MQ=255)
 ccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                           <  1:3116734/101‑1 (MQ=255)
  cccgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                                                          >  1:2206632/1‑101 (MQ=255)
   ccgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCAGGGGGGTGAA.TGGCGCA..GCCAGAggg                                                                                                         <  2:1978608/101‑1 (MQ=255)
    cgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                                        >  1:4450149/1‑101 (MQ=255)
    cgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                                        >  2:45967/1‑101 (MQ=255)
    cgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAgggg                                                                                                        >  2:2556298/1‑101 (MQ=255)
     gtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                                       <  2:1738376/101‑1 (MQ=255)
       cccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                     <  2:1634417/101‑1 (MQ=255)
        ccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                    >  2:2215383/1‑101 (MQ=255)
         cgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   >  2:2785752/1‑101 (MQ=255)
         cgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   >  1:296729/1‑101 (MQ=255)
         cgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   >  1:1718570/1‑101 (MQ=255)
         cgcaagcggACTCCCCCCGGG.ATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                                       <  2:4734452/94‑1 (MQ=18)
           tGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                                >  2:82292/1‑101 (MQ=255)
           tGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                                <  2:4149620/101‑1 (MQ=255)
           tGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                                >  2:717203/1‑101 (MQ=255)
           tGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cc                                                                                                >  2:1468707/1‑100 (MQ=255)
            gCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                                                                              >  1:233057/1‑101 (MQ=255)
             ccGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc                                                                                             >  1:3550859/1‑101 (MQ=255)
              gcggACTCCCCCCGGG.ATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c                                                                                                 <  1:991578/99‑1 (MQ=18)
              ccattcaccccccGGG.ATATTTTCGGCAAGATGAA...AGCC..C..CCCTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCATAGGGGGGCAG.CG.cc                                                                                             >  2:3809676/9‑101 (MQ=11)
                 gACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccc                                                                                       >  1:4029250/1‑100 (MQ=255)
                 gACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccc                                                                                       <  1:737818/101‑2 (MQ=255)
                  aCT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccct                                                                                      >  2:4380080/1‑99 (MQ=255)
                  aCT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccct                                                                                      <  2:3656213/101‑3 (MQ=255)
                   cT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccctt                                                                                     <  1:3089351/101‑4 (MQ=255)
                    cccccccGGGAATTTTGGAGGAAAGAGGAG...GGGG..G..GTTTTCTTTTTGCAAAAATTACCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccct                                                                                      <  1:3981906/99‑3 (MQ=12)
                      cccccGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.cccttcg                                                                                   <  2:677650/101‑6 (MQ=255)
                      cccccGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGAGAA.TGGCGCA..GCCAGAGGGGGGCAT.CG.CCGC.CC..TTcg                                                                                  >  1:1534838/1‑101 (MQ=25)
                       ccccGGG.ATATTTGCGGCAAGATGAA...GGGGCAT..GATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CC.C.CC..tg                                                                                    >  2:4647477/1‑100 (MQ=0)
                        cccGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTcgcg                                                                                <  2:1213634/101‑1 (MQ=255)
                         ccGGG.ATATTTGAGGCAAGATGAA...GGGG..T..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCgcgg                                                                               <  2:305186/101‑1 (MQ=255)
                          cGGG.ATATTTGCGGCAAGATGAA...GGGG..CATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CC.C.CC..tgccc                                                                                 >  2:3553439/1‑97 (MQ=2)
                          cGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCgcggg                                                                              >  1:851451/1‑101 (MQ=255)
                            gg.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCgcgggcg                                                                            >  1:2113736/1‑101 (MQ=255)
                                tatTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGcc                                                                         >  2:2176177/1‑101 (MQ=255)
                                 atTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGccc                                                                        >  1:106802/1‑100 (MQ=255)
                                 atTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCcg                                                                        >  2:1946907/1‑101 (MQ=255)
                                 atTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCcg                                                                        >  2:2656498/1‑101 (MQ=255)
                                  tttGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATGTCCCGGGGGGTGAA.TGGCGGA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCcgc                                                                       <  1:2774647/101‑1 (MQ=37)
                                    tGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGa                                                                     >  1:170387/1‑101 (MQ=255)
                                     gAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGAc                                                                    <  1:2674363/101‑1 (MQ=255)
                                      gggaaaGATGAG...GGGG..G..GTTTTCATCTTCCCGAAAATACCGGGGGGGGGGAA.GGGCGCA..GCCAGGGGGGGGCAGCCC.CC.C.CC..TTCGCGGGCGGCCGCGAcg                                                                   <  1:3329894/97‑1 (MQ=255)
                                          aaGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGAcgcgcg                                                               >  2:1824702/1‑101 (MQ=255)
                                          aaGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGAcgcgcg                                                               >  2:3307546/1‑101 (MQ=255)
                                           aGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACgcgcg.c                                                             >  1:3305551/1‑101 (MQ=255)
                                            gATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGcg.cc                                                            <  1:4351254/101‑1 (MQ=255)
                                            gATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGcg.cc                                                            <  2:2847709/101‑1 (MQ=255)
                                             aTGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGcg.ccg                                                           >  2:3922272/1‑101 (MQ=255)
                                             aTGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGcg.ccg                                                           <  2:1350217/101‑1 (MQ=255)
                                              tGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGcc.ccac                                                          >  2:3969880/1‑99 (MQ=255)
                                                 a...GGGG..G..TCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAAc                                                       <  1:637770/101‑1 (MQ=255)
                                                     gggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACa                                                      <  1:2747345/101‑1 (MQ=255)
                                                     gggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACa                                                      >  1:2461999/1‑101 (MQ=255)
                                                        g..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CC.CGCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGt                                                   >  2:243587/1‑101 (MQ=255)
                                                           g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTc                                                  <  2:4536517/101‑1 (MQ=255)
                                                                tttCTTCTGGCGAAAAATCCCCGGGGGGGTGAA.TGGCGCC..GCCAGGGGGGGGCCG.CC.CC.C.CCCTTTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGt                                                <  2:4029008/101‑1 (MQ=255)
                                                                 ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGc                                              >  2:1795035/1‑101 (MQ=255)
                                                                  tCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGcc                                             >  2:1070557/1‑101 (MQ=255)
                                                                  tCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCGG.CC.CC.C.CC..CTCGCGGGGCGCCCCGACGCGGGCCCGAAACAGGCGGTGcc                                             >  2:2356170/1‑101 (MQ=16)
                                                                   cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGccc                                            <  2:358880/101‑1 (MQ=255)
                                                                    aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGcccc                                           >  2:4121141/1‑101 (MQ=255)
                                                                      cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCtt                                         >  2:3402257/1‑101 (MQ=255)
                                                                      cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCtt                                         <  1:1135972/101‑1 (MQ=255)
                                                                       ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCttt                                        <  2:505808/101‑1 (MQ=255)
                                                                         gCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGt                                      <  1:4658713/101‑1 (MQ=255)
                                                                           cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCa                                    >  1:3733537/1‑101 (MQ=255)
                                                                            gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCAt                                   >  2:2895895/1‑101 (MQ=255)
                                                                             aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCAtt                                  <  1:3643399/101‑1 (MQ=255)
                                                                               aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCt                                >  1:4276132/1‑101 (MQ=255)
                                                                                atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTg                               >  1:2056425/1‑101 (MQ=255)
                                                                                atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTg                               >  1:4587163/1‑101 (MQ=255)
                                                                                atatCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTGGTGCCCCTTTGTCCTTCTg                               >  1:422611/1‑101 (MQ=37)
                                                                                  atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                             <  2:3122330/101‑1 (MQ=255)
                                                                                  atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                             >  2:1144011/1‑101 (MQ=255)
                                                                                    cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCCACCGGTCGTGCCCCGTTGTCATTCTGTCGc                           >  2:4243977/1‑101 (MQ=38)
                                                                                    cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGc                           >  1:1764809/1‑101 (MQ=255)
                                                                                     ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                          <  1:2215507/101‑1 (MQ=255)
                                                                                     ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGCCATTCTGTCGCa                          >  1:3894795/1‑101 (MQ=255)
                                                                                        gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                       >  2:916540/1‑101 (MQ=255)
                                                                                        gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                       >  1:2543109/1‑101 (MQ=255)
                                                                                        gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                       >  1:4103246/1‑101 (MQ=255)
                                                                                         ggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                      <  2:2285914/101‑1 (MQ=255)
                                                                                          gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                     >  2:4208839/1‑101 (MQ=255)
                                                                                            gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                   >  1:1820113/1‑101 (MQ=255)
                                                                                            gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                   <  1:4294995/101‑1 (MQ=255)
                                                                                            gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                   >  2:1314970/1‑101 (MQ=255)
                                                                                              gAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                 >  1:2440840/1‑101 (MQ=255)
                                                                                               aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                <  1:3586526/101‑1 (MQ=255)
                                                                                                a.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa               <  2:4355028/101‑1 (MQ=255)
                                                                                                  tGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa              >  1:1448251/1‑101 (MQ=255)
                                                                                                      gcA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg          >  1:685165/1‑101 (MQ=255)
                                                                                                             cAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc     >  1:1917808/1‑101 (MQ=255)
                                                                                                             cAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc     <  1:1315111/101‑1 (MQ=255)
                                                                                                              agagGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa    >  2:260490/1‑101 (MQ=255)
                                                                                                              agagGGGGGCAG.CG.CC.CGCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGGCATAAGCTGGCCGCa    >  2:3863459/1‑101 (MQ=255)
                                                                                                               gagGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt   <  1:2151801/101‑1 (MQ=255)
                                                                                                                agGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCGGCATc  >  1:4153838/1‑101 (MQ=255)
                                                                                                                |                                                                                                          
GCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CC.C.CC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATC  >  NC_014034/2249280‑2249479

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG