Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,391 +C 90.3% intergenic (‑187/‑118) 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,3911.→C90.3% 38.6 ‑1.8 ‑1.3144intergenic (‑187/‑118)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (75/55):  ref base (5/9):  total (80/64)
Fisher's exact test strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.74e-02

Only 100 of 187 total aligned reads displayed.

TGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCCGAACC  >  NC_014034/2249291‑2249491
                                                                                                                |                                                                                                      
 gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c.g..c                                                                                                      >  1:3413164/1‑101 (MQ=255)
 gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c.g..c                                                                                                      >  2:594752/1‑101 (MQ=255)
    gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..cc.cc                                                                                                  <  2:4247913/101‑1 (MQ=255)
    cggccT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..cc.cc                                                                                                  <  2:3154963/97‑1 (MQ=255)
     ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c.g..c                                                                                                      <  1:2451327/101‑1 (MQ=17)
      gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAG.C.C..cc                                                                                                     <  1:3450227/96‑1 (MQ=17)
       aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..cc.cccct                                                                                               >  1:3062505/1‑99 (MQ=255)
         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..cc.ccccttc                                                                                             >  2:157834/1‑97 (MQ=255)
         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..cc.ccccttc                                                                                             <  1:1780179/101‑5 (MQ=255)
         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..cc.ccccttc                                                                                             <  2:741661/101‑5 (MQ=255)
           cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..cc.ccccttcg                                                                                            <  2:1288833/101‑6 (MQ=255)
           cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGGGCA..GCCAGAGGGGGGCAG.C.G..cc.ccccttcg                                                                                            >  2:3336467/1‑96 (MQ=255)
             cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G..CC.CCCtgc                                                                                              >  2:1905000/1‑99 (MQ=0)
             cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTcgcg                                                                                           >  2:4227192/1‑101 (MQ=255)
                gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCgcgggc                                                                                        >  1:1799212/1‑101 (MQ=255)
                gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCgcgggc                                                                                        >  1:3520456/1‑101 (MQ=255)
                 ggATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G..CC.CCCtgccccg                                                                                          <  2:3054626/101‑7 (MQ=0)
                 ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCgcgggcg                                                                                       >  1:2624461/1‑101 (MQ=255)
                   atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGggcggc                                                                                     <  2:1624353/101‑1 (MQ=255)
                    tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGcc                                                                                    >  2:4603816/1‑101 (MQ=255)
                     atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCcg                                                                                   >  1:4445668/1‑101 (MQ=255)
                     atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCcg                                                                                   >  1:3942320/1‑101 (MQ=255)
                       ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCcgcg                                                                                 >  2:3916653/1‑101 (MQ=255)
                        ttaGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGa                                                                                <  2:2773756/99‑1 (MQ=255)
                        tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGa                                                                                <  2:3821064/101‑1 (MQ=255)
                         gAGACAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGAc                                                                               <  2:4170505/101‑1 (MQ=255)
                           ggCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGAcgc                                                                             <  1:1910438/101‑1 (MQ=255)
                              aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGAcgcgcc                                                                          >  1:516387/1‑100 (MQ=255)
                               aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACgcgcgc                                                                         >  1:3308320/1‑101 (MQ=255)
                               aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACgcgcgc                                                                         <  1:6732/101‑1 (MQ=255)
                                gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGcgcc                                                                        >  2:1985429/1‑101 (MQ=255)
                                gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGcgcc                                                                        <  1:2975977/101‑1 (MQ=255)
                                  tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGcgccgc                                                                      >  2:1996427/1‑101 (MQ=255)
                                  tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGTGAA.AGGCGCA..GCCAGAAGGGGGGAG.CGG..CC.CCCCTCGCCGGCGGCCGCGGCGCGCGCCGCa                                                                     >  2:4505292/1‑101 (MQ=25)
                                    aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGAGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                                    >  1:3556136/1‑101 (MQ=255)
                                     a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGAAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                                   <  1:4115703/101‑1 (MQ=255)
                                         gggg..ggCTTTCATCTTGCCGAAAATATCCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                                  <  1:3371318/101‑1 (MQ=37)
                                         gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAG.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                                  >  1:4271465/1‑101 (MQ=255)
                                         gggg..ggCTTTCATCTTGCCGAAAATAGCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                                  <  1:4622264/101‑1 (MQ=255)
                                          ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                                 >  1:200484/1‑101 (MQ=255)
                                          ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                                 <  2:2208392/101‑1 (MQ=255)
                                          ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                                 >  1:124686/1‑101 (MQ=255)
                                            g..ggCTTTCATCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGt                                                               <  1:3767092/101‑1 (MQ=255)
                                                  tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                           <  1:321237/101‑1 (MQ=255)
                                                  tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                           >  2:112372/1‑101 (MQ=255)
                                                  tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                           >  1:2049931/1‑101 (MQ=255)
                                                   ttCATCTTGCCGAAAATATCCCGGGTGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGGCGTGc                                                          >  1:2608560/1‑101 (MQ=37)
                                                   ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                                          >  1:1385125/1‑101 (MQ=255)
                                                   ttCATCATGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                                          <  1:768392/101‑1 (MQ=255)
                                                     cATCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                                        <  1:4206859/101‑1 (MQ=255)
                                                     cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                                        >  1:2855236/1‑101 (MQ=255)
                                                      aTCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCTCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACATGTCGTGcccc                                                       <  1:4032399/101‑1 (MQ=37)
                                                      aTCTTGCCGAAAATATCCCGGGGGGGGAA.TGGGGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGGCGGGcccc                                                       >  2:1368781/1‑101 (MQ=25)
                                                        cTTGGCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCATAGGGTGGCAG.C.GC.CC.CCCTGCGCGGGCGGCCGCGACTCTCGCCGCACCAGGTCGTGCCCCtt                                                     <  2:2329171/101‑1 (MQ=21)
                                                        cTTGCGGAAAATATCCGGGGGGGTGAA.GGGCGCA..GCCGGGGGGGGGCAGCC.C..CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                                     <  1:997103/101‑1 (MQ=255)
                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                                     >  1:379997/1‑101 (MQ=255)
                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                                     >  2:2626868/1‑101 (MQ=255)
                                                         ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                                    <  1:3628035/101‑1 (MQ=255)
                                                         ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                                    >  2:2635338/1‑101 (MQ=255)
                                                             cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                                <  2:3694959/101‑1 (MQ=255)
                                                             cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                                >  1:1066059/1‑101 (MQ=255)
                                                              gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                               <  1:2321959/101‑1 (MQ=255)
                                                              gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                               >  2:2994249/1‑101 (MQ=255)
                                                               aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                              >  1:2161595/1‑101 (MQ=255)
                                                                aaaTATCGCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                             <  2:2461328/101‑1 (MQ=255)
                                                                 aagcccgtggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                                            <  1:594813/93‑1 (MQ=38)
                                                                 aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                                            >  1:445804/1‑101 (MQ=255)
                                                                  atatCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGCCACGCGCCGCAACAGGTCGTGCCCCTTGGTCAGTCTg                                           >  2:2276934/1‑101 (MQ=37)
                                                                   tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                                          >  2:3475903/1‑101 (MQ=255)
                                                                   tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                                          >  2:2572905/1‑101 (MQ=255)
                                                                       ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                                      >  1:927373/1‑101 (MQ=255)
                                                                        cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                                     <  2:3898244/101‑1 (MQ=255)
                                                                        cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGAGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                                     >  2:1807724/1‑101 (MQ=255)
                                                                          gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                                   >  2:3150915/1‑101 (MQ=255)
                                                                          gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                                   <  1:644304/101‑1 (MQ=255)
                                                                           ggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                                  >  1:4292033/1‑101 (MQ=255)
                                                                           ggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                                  <  2:4086457/101‑1 (MQ=255)
                                                                            gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                                 <  2:4291230/101‑1 (MQ=255)
                                                                             ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                                <  2:1147735/101‑1 (MQ=255)
                                                                               tGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc                              >  2:1055919/1‑101 (MQ=255)
                                                                                gAA.TGGCGCA..GCCAGGGGGGGGCAGCC.C..CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                             <  1:4066306/101‑1 (MQ=255)
                                                                                    tGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                          <  2:3062261/101‑1 (MQ=255)
                                                                                    tGGCGCA..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                          <  2:649149/101‑1 (MQ=255)
                                                                                    gggCGCA..GCGAGGGGGGGGCAGCC.C..CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                          <  1:157848/100‑1 (MQ=255)
                                                                                          a..GCCAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc                    >  2:2905234/1‑101 (MQ=255)
                                                                                             gCCAGGGGGGGGCCC.C.GCCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc                    <  1:1658215/101‑1 (MQ=255)
                                                                                              ccggagGGGGGCAGCC.C..CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg                  <  1:1049738/98‑1 (MQ=255)
                                                                                              ccAGAGGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg                  >  2:612521/1‑101 (MQ=255)
                                                                                                agagGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa                >  2:3520522/1‑101 (MQ=255)
                                                                                                agagGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa                >  2:897599/1‑101 (MQ=255)
                                                                                                 gagGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt               <  1:349240/101‑1 (MQ=255)
                                                                                                  agGGGGGCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATc              <  1:3819061/101‑1 (MQ=255)
                                                                                                    gggggCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGa            >  1:526913/1‑101 (MQ=255)
                                                                                                     ggggCAGCC.C..CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGAt           <  2:254508/101‑1 (MQ=255)
                                                                                                     ggggCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGAt           >  1:690054/1‑101 (MQ=255)
                                                                                                     ggggCAG.C.G..CCACCCTTCGCGGGCGGCCGCGACGTGCGCGGCAACAGGTCGTGCCCCTTTGACATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGAt           >  1:3851508/1‑101 (MQ=255)
                                                                                                     gcggCAG.C.GC.CC.CCCTGCGCGGGCGGCCCCGCCGCGCGCCGCACCAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGAt           <  2:516339/99‑1 (MQ=255)
                                                                                                      gggCAG.C.CC.CC.CCCTCCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATg          <  1:1736739/101‑1 (MQ=255)
                                                                                                        gCAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTc        <  2:3412910/101‑1 (MQ=255)
                                                                                                         cAG.C.GC.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTcc       >  2:3985544/1‑101 (MQ=255)
                                                                                                                |                                                                                                      
TGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G..CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCCGAACC  >  NC_014034/2249291‑2249491

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