Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,379 G→A 2.4% intergenic (‑175/‑130) 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,3790G→A 2.4% 5.0 0.0 ‑1.0209intergenic (‑175/‑130)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (3/2):  ref base (114/90):  total (117/92)
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.14e-02

Only 100 of 240 total aligned reads displayed.

TGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.C.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAT  >  NC_014034/2249279‑2249478
                                                                                                               |                                                                                                          
gcccccGTCCCCTGCCGGGGCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                            <  1:2817521/99‑1 (MQ=255)
   cccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                                                         >  1:2206632/1‑101 (MQ=255)
     cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                                       >  1:4450149/1‑101 (MQ=255)
     cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAgggg                                                                                                       >  2:2556298/1‑101 (MQ=255)
         ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                   >  2:2215383/1‑101 (MQ=255)
         ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGAGGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                   <  2:3214748/101‑1 (MQ=255)
          cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                  <  2:285407/101‑1 (MQ=255)
          cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                  >  2:2785752/1‑101 (MQ=255)
          cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                  >  1:1718570/1‑101 (MQ=255)
          cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                                      <  1:3922950/94‑1 (MQ=18)
           gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                                                     <  2:4360232/95‑1 (MQ=18)
           gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                                                     <  1:3576625/95‑1 (MQ=18)
            tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                               >  2:717203/1‑101 (MQ=255)
            tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                               >  2:82292/1‑101 (MQ=255)
            tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cc                                                                                               >  2:1468707/1‑100 (MQ=255)
             gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                                                                             >  1:233057/1‑101 (MQ=255)
               gcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c                                                                                                <  1:991578/99‑1 (MQ=18)
               ccattcaccccccGGGATATTTTCGGCAAGATGAA...AGCC..C..CCCTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCATAGGGGGGCAG.CG.c.c                                                                                           >  2:3809676/9‑101 (MQ=11)
                  gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.c.cccc                                                                                      <  1:737818/101‑2 (MQ=255)
                    cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTCCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.c.cccctt                                                                                    <  1:4020275/101‑4 (MQ=255)
                       cccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.c.ccccttcg                                                                                  >  1:812051/1‑96 (MQ=255)
                       cccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.c.ccccttcg                                                                                  <  1:4528552/101‑6 (MQ=255)
                       cccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.c.c.ccccttcg                                                                                  <  2:677650/101‑6 (MQ=255)
                       cccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGAGAA.TGGCGCA..GCCAGAGGGGGGCAT.CG.C.CGCCC..TTcg                                                                                 >  1:1534838/1‑101 (MQ=25)
                         cccGGGATATTTGCGGCAAGATGAA...GGGGCAT..GATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.C.C.CCC..tgc                                                                                  <  2:4544891/101‑3 (MQ=0)
                         cccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTcgcg                                                                               <  2:1213634/101‑1 (MQ=255)
                          ccGGGATATTTGAGGCAAGATGAA...GGGG..T..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCgcgg                                                                              <  2:305186/101‑1 (MQ=255)
                          ccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCgcgg                                                                              <  2:607886/101‑1 (MQ=255)
                           cGGGATATTTGCGGCAAGATGAA...GGGGCAT..GATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.C.C.CCC..tgccc                                                                                >  2:2731666/1‑97 (MQ=0)
                           cGGGATATTTGCGGCAAGATGAA...GGGG..CATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.C.C.CCC..tgccc                                                                                >  2:3553439/1‑97 (MQ=2)
                            gggATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCgcgggc                                                                            >  2:487361/1‑101 (MQ=255)
                             ggATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCgcgggcg                                                                           >  1:2113736/1‑101 (MQ=255)
                               atatTTGAGGCAAGATGAA...GGGG..G..GCTTTCACCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGggcggc                                                                         >  1:4206156/1‑101 (MQ=255)
                                tatTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGcc                                                                        >  2:2176177/1‑101 (MQ=255)
                                tatTTGAGGCAAGAGGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGcc                                                                        <  1:3167913/101‑1 (MQ=255)
                                 atTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGccc                                                                       >  1:106802/1‑100 (MQ=255)
                                 atTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCcg                                                                       >  2:1946907/1‑101 (MQ=255)
                                  tttGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATGTCCCGGGGGGTGAA.TGGCGGA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCcgc                                                                      <  1:2774647/101‑1 (MQ=37)
                                    tGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGa                                                                    >  1:170387/1‑101 (MQ=255)
                                    tGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGa                                                                    >  1:3362466/1‑101 (MQ=255)
                                     gAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGAc                                                                   <  1:2674363/101‑1 (MQ=255)
                                          aaGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGAcgcgcg                                                              >  2:3307546/1‑101 (MQ=255)
                                           aGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACgcgcg.c                                                            >  2:3871958/1‑101 (MQ=255)
                                           aGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACgcgcg.c                                                            >  2:1297249/1‑101 (MQ=255)
                                            gATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGcg.cc                                                           <  1:4351254/101‑1 (MQ=255)
                                             aTGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGcg.ccg                                                          >  2:3922272/1‑101 (MQ=255)
                                              tGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGcc.ccac                                                         >  2:3969880/1‑99 (MQ=255)
                                                 a...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAAc                                                      <  2:244620/101‑1 (MQ=255)
                                                     gggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACa                                                     >  1:2461999/1‑101 (MQ=255)
                                                     gggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACa                                                     <  1:2747345/101‑1 (MQ=255)
                                                       gg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAgg                                                   >  2:4531718/1‑101 (MQ=255)
                                                           g..gCTTTGATCTTGCCGAAACTAGCCCGGGGGGTGAA.TGGCGCA..GCCAGGGGGGGGCGG.CGCC.C.CCC..TTCGCGGGCGGCCGCGAGGCGCG.CCGCAACAGGTc                                                 <  2:4687873/101‑1 (MQ=25)
                                                           g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTc                                                 <  1:3663229/101‑1 (MQ=255)
                                                              gCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.GGCC.C.CCC..TTCCCGGGCGGCCGCGCGGCGCG.CCGCacaaggccg                                                >  2:1752640/1‑93 (MQ=25)
                                                                tttCTTCTGGCGAAAAATCCCCGGGGGGGTGAA.TGGCGCC..GCCAGGGGGGGGCCG.CC.C.C.CCCCTTTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGt                                               <  2:4029008/101‑1 (MQ=255)
                                                                 ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGc                                             >  2:1795035/1‑101 (MQ=255)
                                                                  tCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCGG.CC.C.C.CCC..CTCGCGGGGCGCCCCGACGCGGGCCCGAAACAGGCGGTGcc                                            >  2:2356170/1‑101 (MQ=16)
                                                                   cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGccc                                           <  2:358880/101‑1 (MQ=255)
                                                                    aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGcccc                                          >  2:3437060/1‑101 (MQ=255)
                                                                    aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGcccc                                          >  2:4121141/1‑101 (MQ=255)
                                                                      cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCtt                                        <  1:1135972/101‑1 (MQ=255)
                                                                       ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCAAGAGGGGGGAAG.AG.CAA.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCttt                                       >  1:607164/1‑101 (MQ=38)
                                                                         gCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGt                                     <  1:4658713/101‑1 (MQ=255)
                                                                         gCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGt                                     >  2:230904/1‑101 (MQ=255)
                                                                          cgaaaaaTATCGCGGGGGGTGAA.GGGCGCA..GCCAGGGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTc                                    <  2:1688591/98‑1 (MQ=38)
                                                                           cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCa                                   >  2:4600832/1‑101 (MQ=255)
                                                                            gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCAt                                  <  1:3002827/101‑1 (MQ=255)
                                                                            gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCAt                                  >  2:2895895/1‑101 (MQ=255)
                                                                             aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGATGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCAtt                                 <  2:1853448/101‑1 (MQ=255)
                                                                             aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCAtt                                 >  1:9936/1‑101 (MQ=255)
                                                                              aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGCCGTGCCCCTTTGTCATTc                                >  1:3641255/1‑101 (MQ=255)
                                                                               aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCt                               >  1:4276132/1‑101 (MQ=255)
                                                                                atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTg                              >  1:2056425/1‑101 (MQ=255)
                                                                                  atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                            >  2:1144011/1‑101 (MQ=255)
                                                                                    cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCCACCGGTCGTGCCCCGTTGTCATTCTGTCGc                          >  2:4243977/1‑101 (MQ=38)
                                                                                    cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGc                          >  1:1764809/1‑101 (MQ=255)
                                                                                     ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                         <  1:2215507/101‑1 (MQ=255)
                                                                                        gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                      >  2:916540/1‑101 (MQ=255)
                                                                                        gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                      >  1:2543109/1‑101 (MQ=255)
                                                                                         ggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                     >  2:3721327/1‑101 (MQ=255)
                                                                                          gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                    >  2:4208839/1‑101 (MQ=255)
                                                                                          gggTGAA.GGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                    <  1:382215/101‑1 (MQ=255)
                                                                                           ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                   >  1:3321197/1‑101 (MQ=255)
                                                                                            gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                  <  1:4294995/101‑1 (MQ=255)
                                                                                            gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                  >  1:1820113/1‑101 (MQ=255)
                                                                                             tGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc                 >  2:1388287/1‑101 (MQ=255)
                                                                                              gAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                >  1:2440840/1‑101 (MQ=255)
                                                                                               aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt               <  2:2042507/101‑1 (MQ=255)
                                                                                                aagggCGCT..GCCAGGGGGGGGCAGCCC.C.C.CCC..TGCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCTCCCATGTCAt               <  1:4102709/98‑1 (MQ=255)
                                                                                                  tGGCGCA..GCCAGAGGGGGGCAGCCC.C.C.CCC..GTCGCGGGCGGCCACGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa             <  1:310059/101‑1 (MQ=255)
                                                                                                  tGGCGCA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa             >  1:1448251/1‑101 (MQ=255)
                                                                                                     cgcG..GCCGGAGGGGGGCAGCCC.C.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt          <  1:4136294/101‑1 (MQ=255)
                                                                                                      gcA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg         >  1:685165/1‑101 (MQ=255)
                                                                                                       cA..GCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg        <  2:4157753/101‑1 (MQ=255)
                                                                                                           gCCAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGcc      >  2:4331004/1‑101 (MQ=255)
                                                                                                             cAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc    >  1:1917808/1‑101 (MQ=255)
                                                                                                             cAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc    >  1:3520520/1‑101 (MQ=255)
                                                                                                             cAGAGGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc    <  1:1315111/101‑1 (MQ=255)
                                                                                                              agagGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa   >  2:260490/1‑101 (MQ=255)
                                                                                                               gagGGGGGCAG.CGCC.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt  <  1:2151801/101‑1 (MQ=255)
                                                                                                               |                                                                                                          
TGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.C.C.CCC..TTCGCGGGCGGCCGCGACGCGCG.CCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAT  >  NC_014034/2249279‑2249478

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