Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,391 +C 94.1% 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.→C94.1% 28.1 0.0 ‑1.6203intergenic (‑187/‑118)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (115/76):  ref base (7/5):  total (122/81)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.43e-02

Only 100 of 244 total aligned reads displayed.

TGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG..CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCCGA  >  NC_014034/2249291‑2249488
                                                                                                          |                                                                                                   
tGCCTGGACTCCCCCGGGATATTTGAGGCAAGATCAG...GGGGGGTTTTCATCTT.TCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAg.cg                                                                                                      <  1:2240441/101‑1 (MQ=37)
 gCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAAGGGCGCAGCCAGAGGGGGGCAg.ag..c                                                                                                   >  2:1737242/1‑101 (MQ=255)
 gCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTAATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAg.cg..c                                                                                                   >  2:2231113/1‑101 (MQ=255)
   cGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGGGAATGGCGCAGCCAGAGGGGGGCAG.CG..ccc                                                                                                 >  1:1052733/1‑101 (MQ=255)
     gactccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTT.GCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAGGGGGGCAG.CG..cc                                                                                                  <  1:4232038/96‑1 (MQ=17)
      gACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG..cccccc                                                                                              <  1:2807368/101‑2 (MQ=255)
       aCTCCCCAGGGTTATTTGAGGCAAGATGTA...GGGGGGCTTTCATCTT.GGCGAAAATTTGCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG..cccccct                                                                                             <  2:4101922/101‑3 (MQ=25)
        cTCCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG..cccccctt                                                                                            <  1:4041099/101‑4 (MQ=255)
               gggATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCgcgggc                                                                                     <  1:654632/101‑1 (MQ=255)
                ggATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCgcggggg                                                                                    >  1:3003233/1‑99 (MQ=255)
                 gATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCgcgggcgg                                                                                   >  1:4569746/1‑101 (MQ=255)
                  atatTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTAAATGGCGCAGACAGAGGGGGGCAG.CGC.CCCCC.TTCGCGggcggc                                                                                  >  2:901926/1‑101 (MQ=37)
                  atatTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGCGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGggcggc                                                                                  >  1:74012/1‑101 (MQ=255)
                   tatTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGcc                                                                                 <  2:2757136/101‑1 (MQ=255)
                     tttGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCcgc                                                                               >  1:4305555/1‑101 (MQ=255)
                     tttGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCcgc                                                                               <  2:4880746/101‑1 (MQ=255)
                      ttGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCcgcg                                                                              >  2:1926736/1‑101 (MQ=255)
                      ttGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCcgcg                                                                              >  1:4922108/1‑101 (MQ=255)
                       tGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGa                                                                             <  1:3864265/101‑1 (MQ=255)
                       tGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGCGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGGGGGCGGCCGCGa                                                                             >  1:477324/1‑101 (MQ=37)
                         aggtaggaggAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGAcg                                                                           <  2:1582857/92‑1 (MQ=25)
                             aaGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGAcgcgcg                                                                       >  2:3936860/1‑101 (MQ=255)
                                aTGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGcgccg                                                                    >  1:2085237/1‑101 (MQ=255)
                                 tGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGcgccgc                                                                   >  2:1777897/1‑101 (MQ=255)
                                  gAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                                  >  2:2673307/1‑101 (MQ=255)
                                  gAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                                  >  2:3637792/1‑101 (MQ=255)
                                  gAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                                  <  2:4682148/101‑1 (MQ=255)
                                  gAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                                  >  2:515019/1‑101 (MQ=255)
                                  gAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                                  <  2:540454/101‑1 (MQ=255)
                                   ta...GGGGGGGCTTCATATTCCCCGAAAATAGCGTGGGGGGTAAACGGCGCAGCCAGGGGGGGGCGGCCC..CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                                  <  1:3078535/100‑1 (MQ=255)
                                   aa...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCGGCaa                                                                 >  1:2565231/1‑101 (MQ=255)
                                   aa...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                                 >  2:1052694/1‑101 (MQ=255)
                                   aa...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                                 >  2:249070/1‑101 (MQ=255)
                                   aa...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                                 >  1:1143316/1‑101 (MQ=255)
                                    a...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                                <  2:2403842/101‑1 (MQ=255)
                                    a...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                                >  2:1840746/1‑101 (MQ=255)
                                        ggggggCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                               >  2:3529936/1‑101 (MQ=255)
                                        ggggggCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                               <  1:2124760/101‑1 (MQ=255)
                                              cTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                                         <  2:3826913/101‑1 (MQ=255)
                                               tttCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                        >  1:2365250/1‑101 (MQ=255)
                                               tttCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                        <  2:1218367/101‑1 (MQ=255)
                                                ttCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                                       <  2:4382807/101‑1 (MQ=255)
                                                 tCATCTT.GCCGAAAATATCCGGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                                      <  2:1296295/101‑1 (MQ=255)
                                                 tCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                                      >  1:3201183/1‑101 (MQ=255)
                                                   aTCTT.GCCGAAAATATCCCGGGGGGTGAAGGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                                    <  1:3981074/101‑1 (MQ=255)
                                                    tCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                                   <  1:2859550/101‑1 (MQ=255)
                                                    tCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                                   >  2:2699450/1‑101 (MQ=255)
                                                     cTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                                  <  2:3394470/101‑1 (MQ=255)
                                                     cTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                                  >  1:4349614/1‑101 (MQ=255)
                                                          ccGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                              <  2:4216837/101‑1 (MQ=255)
                                                           cGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                             >  1:2150742/1‑101 (MQ=255)
                                                            gAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                            >  1:1794654/1‑101 (MQ=255)
                                                             aaaaTATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                           >  1:3784659/1‑101 (MQ=255)
                                                             aaaaTATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                           >  2:5019650/1‑101 (MQ=255)
                                                              aaaTATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                          <  2:308107/101‑1 (MQ=255)
                                                              aaaTATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                          >  1:4295762/1‑101 (MQ=255)
                                                               aaTTTCCGGGGGGGTGAATGGCGCAGCCAAAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGCGCCCCTTTGTTATTCt                                         >  2:165061/1‑101 (MQ=37)
                                                               aaTGTCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAC.CGC.CCCCC.TTCGCTGCCGACCGCGACGCGCGCCCCACCTGGTCGTGCCCCTATGTCATTCt                                         >  2:1642620/1‑101 (MQ=255)
                                                                 tatCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCC.TTCGCGGGCGCCCGCACCCCCGGCCGCAAAAGGCCGGGCCCCTTGTTCTTTCtt                                        >  2:588032/1‑100 (MQ=12)
                                                                  atCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                                      <  1:4945884/101‑1 (MQ=255)
                                                                  atCCCGGGGGATGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                                      >  2:432316/1‑101 (MQ=255)
                                                                      gggggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                                  <  2:1009809/100‑1 (MQ=255)
                                                                      cGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                                  >  2:2894041/1‑101 (MQ=255)
                                                                        gggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCTCCccc                                >  2:3170935/1‑101 (MQ=38)
                                                                        gggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                                >  2:1458945/1‑101 (MQ=255)
                                                                        gggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                                <  2:3985285/101‑1 (MQ=255)
                                                                        gggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                                >  1:702558/1‑101 (MQ=255)
                                                                          gggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                              >  2:2503832/1‑101 (MQ=255)
                                                                          gggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                              >  1:3478696/1‑101 (MQ=255)
                                                                          gggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                              >  2:4470104/1‑101 (MQ=255)
                                                                           ggTGAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TGCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                             <  1:4829377/101‑1 (MQ=255)
                                                                             ggAATGGCGCAGCCAGAGGGGGGCAGCCC..CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc                           <  2:3478431/100‑1 (MQ=255)
                                                                              gAATGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                          >  1:4527252/1‑101 (MQ=255)
                                                                               aaTGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                         <  2:2274322/101‑1 (MQ=255)
                                                                               aaTGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGCCATTCTGTCGCACCCATGTCAt                         >  1:2320318/1‑101 (MQ=255)
                                                                                aTGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa                        >  1:1413744/1‑101 (MQ=255)
                                                                                aTGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa                        <  2:2427176/101‑1 (MQ=255)
                                                                                 tGGCGCAGCGGGGGGGGGGCGG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                       <  1:579683/101‑1 (MQ=255)
                                                                                 tGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                       >  2:2808364/1‑101 (MQ=255)
                                                                                 tGGCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                       >  1:4413198/1‑101 (MQ=255)
                                                                                  ggCGCAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg                      <  1:4976029/101‑1 (MQ=255)
                                                                                     gcAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg                   <  1:2936855/101‑1 (MQ=255)
                                                                                      cAGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg                  >  2:3420182/1‑101 (MQ=255)
                                                                                       aGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc                 >  2:4603970/1‑101 (MQ=255)
                                                                                       aGCCAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc                 >  1:3435802/1‑101 (MQ=255)
                                                                                          cAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc              >  1:2070822/1‑101 (MQ=255)
                                                                                          cAGAGGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc              >  1:905595/1‑101 (MQ=255)
                                                                                          cAGAGGGGGCCGGCCC..CCCCC.TCCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc              <  2:2003903/101‑1 (MQ=255)
                                                                                           agagGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa             >  2:3004394/1‑101 (MQ=255)
                                                                                            gagGGGGGCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt            >  1:1576532/1‑101 (MQ=255)
                                                                                              ggggggTAG.CC..CCCCCGTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTTCCCCTTTGTCATTCTGTCGCGCCCATGTCATAAGCTGGCCGCATCg          <  1:2014668/101‑1 (MQ=255)
                                                                                              ggggggCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCg          <  2:3993861/101‑1 (MQ=255)
                                                                                              ggggggCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCg          <  2:875741/101‑1 (MQ=255)
                                                                                               gggggCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGa         <  1:363917/101‑1 (MQ=255)
                                                                                                 gggCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATg       >  2:1152696/1‑101 (MQ=255)
                                                                                                   gCAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTc     <  2:3197998/101‑1 (MQ=255)
                                                                                                    cAG.CGC.CCCCC.TTGGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTcc    <  2:1550676/101‑1 (MQ=255)
                                                                                                    cAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTcc    >  1:4485527/1‑101 (MQ=255)
                                                                                                    cAG.CGC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTcc    >  2:4766799/1‑101 (MQ=255)
                                                                                                      gcgg..cCCCC.TTCGCGGGCGTCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCCGa  <  2:3808867/98‑1 (MQ=255)
                                                                                                          |                                                                                                   
TGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTT.GCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG..CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCCGA  >  NC_014034/2249291‑2249488

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG