Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,381 G→A 2.4% intergenic (‑177/‑128) 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,3810G→A 2.4% 4.4 0.0 ‑0.0207intergenic (‑177/‑128)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (3/2):  ref base (109/92):  total (112/95)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.52e-01

Only 100 of 241 total aligned reads displayed.

CCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCG  >  NC_014034/2249281‑2249480
                                                                                                              |                                                                                                         
ccgcgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                           >  1:3510041/4‑101 (MQ=255)
ccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATCTCGCGGGGTGTGAA.TGGCGCA..GCCagag                                                                                                           <  2:3586635/101‑1 (MQ=255)
ccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                           <  1:3116734/101‑1 (MQ=255)
 cccgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                                                          >  1:2206632/1‑101 (MQ=255)
   cgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                                        >  2:81113/1‑101 (MQ=255)
   cgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                                        >  1:4450149/1‑101 (MQ=255)
    gtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                                       <  2:1738376/101‑1 (MQ=255)
      cccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                     <  2:1634417/101‑1 (MQ=255)
       ccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                    >  2:4528299/1‑101 (MQ=255)
       ccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                    <  2:3214748/101‑1 (MQ=255)
        cgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   >  1:296729/1‑101 (MQ=255)
        cgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   <  2:285407/101‑1 (MQ=255)
        cgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                   >  1:1718570/1‑101 (MQ=255)
         gcaagcggACTCCCCCCGGG.ATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                                                      <  2:4360232/95‑1 (MQ=18)
          tGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                                <  2:4149620/101‑1 (MQ=255)
          tGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cc                                                                                                >  2:1468707/1‑100 (MQ=255)
           gCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                                                                              >  1:233057/1‑101 (MQ=255)
            ccGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc                                                                                             >  1:3550859/1‑101 (MQ=255)
             gcggACTCCCCCCGGG.ATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c                                                                                                 <  1:991578/99‑1 (MQ=18)
             ccattcaccccccGGG.ATATTTTCGGCAAGATGAA...AGCC..CCCCTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCATAGGGGGGCAG.CG.cc                                                                                             >  2:3809676/9‑101 (MQ=11)
              gggACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.c                                                                                         >  2:2239127/1‑101 (MQ=255)
                gACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccc                                                                                       <  1:737818/101‑2 (MQ=255)
                gACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccc                                                                                       >  1:4029250/1‑100 (MQ=255)
                  cT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccctt                                                                                     <  1:3089351/101‑4 (MQ=255)
                   t.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.cccttc                                                                                    <  2:3059914/101‑5 (MQ=255)
                   cccccccGGGAATTTTGGAGGAAAGAGGAG...GGGG..GGTTTTCTTTTTGCAAAAATTACCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc.c.ccct                                                                                      <  1:3981906/99‑3 (MQ=12)
                     cccccGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGAGAA.TGGCGCA..GCCAGAGGGGGGCAT.CG.CCGC.CC..TTcg                                                                                  >  1:1534838/1‑101 (MQ=25)
                       cccGGG.ATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CC.C.CC..tgc                                                                                   <  2:4544891/101‑3 (MQ=0)
                       cccGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTcgcg                                                                                <  2:1213634/101‑1 (MQ=255)
                        ccGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCgcgg                                                                               <  2:607886/101‑1 (MQ=255)
                         cGGG.ATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CC.C.CC..tgccc                                                                                 >  2:2731666/1‑97 (MQ=0)
                         cGGG.ATATTTGAGGCAAGATGAA...GGGG..GTCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCgcggg                                                                              <  2:180712/101‑1 (MQ=255)
                         cGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCgcggg                                                                              >  1:851451/1‑101 (MQ=255)
                          ggg.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCgcgggc                                                                             >  2:487361/1‑101 (MQ=255)
                           gg.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCgcgggcg                                                                            >  1:2113736/1‑101 (MQ=255)
                               tatTTGAGGCAAGATGAG...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGcc                                                                         <  2:3019066/101‑1 (MQ=255)
                                atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGccc                                                                        >  1:106802/1‑100 (MQ=255)
                                 tttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATGTCCCGGGGGGTGAA.TGGCGGA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCcgc                                                                       <  1:2774647/101‑1 (MQ=37)
                                  ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCcgcg                                                                      <  2:264081/101‑1 (MQ=255)
                                   tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGa                                                                     >  1:170387/1‑101 (MQ=255)
                                    gAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGAc                                                                    <  1:2674363/101‑1 (MQ=255)
                                     gggaaaGATGAG...GGGG..GGTTTTCATCTTCCCGAAAATACCGGGGGGGGGGAA.GGGCGCA..GCCAGGGGGGGGCAGCCC.CC.C.CC..TTCGCGGGCGGCCGCGAcg                                                                   <  1:3329894/97‑1 (MQ=255)
                                          aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACgcgcgc                                                              >  1:3305551/1‑101 (MQ=255)
                                           gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGcgcc                                                             <  1:4351254/101‑1 (MQ=255)
                                            aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGcgccg                                                            <  2:1350217/101‑1 (MQ=255)
                                            aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGcgccg                                                            >  2:3922272/1‑101 (MQ=255)
                                             tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGccccac                                                           >  2:3969880/1‑99 (MQ=255)
                                                a...GGGG..GTCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                        <  1:637770/101‑1 (MQ=255)
                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                        <  2:334307/101‑1 (MQ=255)
                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                        <  2:244620/101‑1 (MQ=255)
                                                    gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                       >  1:2461999/1‑101 (MQ=255)
                                                    gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                       <  1:2747345/101‑1 (MQ=255)
                                                          ggCTTTGATCTTGCCGAAACTAGCCCGGGGGGTGAA.TGGCGCA..GCCAGGGGGGGGCGG.CGCCC.C.CC..TTCGCGGGCGGCCGCGAGGCGCGCCGCAACAGGTc                                                   <  2:4687873/101‑1 (MQ=25)
                                                             tttCTTCTGGCGAAAAATCCCCGGGGGGGTGAA.TGGCGCC..GCCAGGGGGGGGCCG.CC.CC.C.CCCTTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                                 <  2:4029008/101‑1 (MQ=255)
                                                               tCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                              >  2:1070557/1‑101 (MQ=255)
                                                                 aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                            >  2:3437060/1‑101 (MQ=255)
                                                                 aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                            >  2:4121141/1‑101 (MQ=255)
                                                                   cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                          <  1:1135972/101‑1 (MQ=255)
                                                                      gCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                                       <  1:4658713/101‑1 (MQ=255)
                                                                      gCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                                       >  2:230904/1‑101 (MQ=255)
                                                                        cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                     >  2:4600832/1‑101 (MQ=255)
                                                                        cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                     >  1:3733537/1‑101 (MQ=255)
                                                                          aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGATGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                   <  2:1853448/101‑1 (MQ=255)
                                                                          aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                   <  1:3643399/101‑1 (MQ=255)
                                                                            aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                                 >  1:4276132/1‑101 (MQ=255)
                                                                             atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                                >  1:4587163/1‑101 (MQ=255)
                                                                             atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                                >  1:2056425/1‑101 (MQ=255)
                                                                             atatCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTGGTGCCCCTTTGTCCTTCTg                                >  1:422611/1‑101 (MQ=37)
                                                                               atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                              >  2:1144011/1‑101 (MQ=255)
                                                                                tCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCGTTGTCATTCTGttg                             >  2:2566277/1‑99 (MQ=255)
                                                                                 cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCCACCGGTCGTGCCCCGTTGTCATTCTGTCGc                            >  2:4243977/1‑101 (MQ=38)
                                                                                 cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGc                            >  1:1764809/1‑101 (MQ=255)
                                                                                  ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                           <  1:2215507/101‑1 (MQ=255)
                                                                                  ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGCCATTCTGTCGCa                           >  1:3894795/1‑101 (MQ=255)
                                                                                    ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCGGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                         >  2:2810997/1‑101 (MQ=255)
                                                                                     gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                        >  2:916540/1‑101 (MQ=255)
                                                                                     gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                        >  1:2543109/1‑101 (MQ=255)
                                                                                     gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                        >  1:4103246/1‑101 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                      <  2:2206508/101‑1 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                      >  2:4208839/1‑101 (MQ=255)
                                                                                        ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                     >  2:239263/1‑101 (MQ=255)
                                                                                         gggAA.TGGCGCC..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                    <  2:1970990/99‑1 (MQ=255)
                                                                                         gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                    <  1:4294995/101‑1 (MQ=255)
                                                                                         gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                    >  1:1820113/1‑101 (MQ=255)
                                                                                         gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                    >  2:1314970/1‑101 (MQ=255)
                                                                                           gAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                  >  1:2440840/1‑101 (MQ=255)
                                                                                            aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                 <  1:3586526/101‑1 (MQ=255)
                                                                                            aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                 <  2:2042507/101‑1 (MQ=255)
                                                                                            aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                 >  2:68389/1‑101 (MQ=255)
                                                                                               tGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa               <  2:3635237/101‑1 (MQ=255)
                                                                                               tGGCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa               >  1:1448251/1‑101 (MQ=255)
                                                                                                ggCGCA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg              >  2:836456/1‑101 (MQ=255)
                                                                                                   gcA..GCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg           >  1:685165/1‑101 (MQ=255)
                                                                                                        gCCAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGcc        >  2:4331004/1‑101 (MQ=255)
                                                                                                          cAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc      >  1:1917808/1‑101 (MQ=255)
                                                                                                          cAGAGGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc      <  1:1315111/101‑1 (MQ=255)
                                                                                                           agagGGGGGCAG.CG.CC.CGCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGGCATAAGCTGGCCGCa     >  2:3863459/1‑101 (MQ=255)
                                                                                                            gagGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt    <  1:2151801/101‑1 (MQ=255)
                                                                                                             agGGGGGCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCGGCATc   >  1:4153838/1‑101 (MQ=255)
                                                                                                              ggggggCAG.CGCCC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTTATAAGCTGGCCGCATCg  <  2:366361/101‑1 (MQ=255)
                                                                                                              |                                                                                                         
CCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CC.C.CC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCG  >  NC_014034/2249281‑2249480

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 39 ≤  ATCG