Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,391 +C 94.0% 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.0% 32.3 ‑2.4 ‑0.2216intergenic (‑187/‑118)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (113/90):  ref base (4/9):  total (117/99)
Fisher's exact test strand distribution p-value = 9.24e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.01e-01

Only 100 of 255 total aligned reads displayed.

TGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCC  >  NC_014034/2249291‑2249486
                                                                                                               |                                                                                                
tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                                   >  2:5475188/1‑101 (MQ=255)
 gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                                                                                 >  1:2606070/1‑101 (MQ=255)
  ccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc                                                                                                >  2:2292392/1‑101 (MQ=255)
  ccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc                                                                                                <  2:171451/101‑1 (MQ=255)
     ggccT.CCCCCGGGATATTGGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                                                                             <  1:1024759/98‑1 (MQ=255)
     ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGCGAA..GCCAGAGGGGGGCTg.cg.c                                                                                                 >  2:4956631/1‑101 (MQ=17)
     ggACT.CCCCGGGGATATTTGAGGCAAGATGAA...GGGG..GCCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                                                                             <  2:5615259/101‑1 (MQ=255)
     ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                                                                             >  1:4269649/1‑101 (MQ=255)
      gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAG.CG.cc                                                                                                <  2:1610315/96‑1 (MQ=17)
         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                                                                         >  2:580828/1‑97 (MQ=255)
         t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                                                                         >  1:3511685/1‑97 (MQ=255)
             cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTcgcg                                                                                      >  2:3663376/1‑101 (MQ=255)
             cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTcgcg                                                                                      >  2:1122441/1‑101 (MQ=255)
              ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCgcgg                                                                                     <  2:1759336/101‑1 (MQ=255)
               cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCC.tgccc                                                                                       >  2:4065635/1‑97 (MQ=2)
                gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCgcgggc                                                                                   >  1:2198455/1‑101 (MQ=255)
                gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCgcgggc                                                                                   >  1:4592779/1‑101 (MQ=255)
                gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCgcgggc                                                                                   <  2:4246316/101‑1 (MQ=255)
                 ggATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCC.tgcaccg                                                                                     <  2:4900346/101‑7 (MQ=0)
                   atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGggcggc                                                                                >  2:5111094/1‑101 (MQ=255)
                    tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGcc                                                                               <  1:5614770/101‑1 (MQ=255)
                         gAGGCAAGATGAA...GGGG..TGCGTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGAc                                                                          <  2:4309113/101‑1 (MQ=25)
                             cAAGATGAA...GGGG..GCCTTTCATCTGGCCGAAAATATGCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGAcgcgc                                                                      <  2:2175851/101‑1 (MQ=37)
                              aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGAcgcgcg                                                                     <  1:553878/101‑1 (MQ=255)
                              aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGAccggcg                                                                     >  1:3044340/1‑96 (MQ=25)
                               aGATGAG...GGGG..TGTTTTAATCTGGCAAAAAATACCCGGGGGGGTGAA.GGGCGCA..GCCAGGGGGGGGCAGCCC.CCCCC.TTCGCGGGCGGCCGCGACgcgcgc                                                                    <  2:3958692/101‑1 (MQ=255)
                               aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACgcgcgc                                                                    >  1:5437517/1‑101 (MQ=255)
                                gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGcgcc                                                                   >  1:4795376/1‑101 (MQ=255)
                                gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGcgcc                                                                   >  1:5389961/1‑101 (MQ=255)
                                gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGcgcc                                                                   <  2:3335113/101‑1 (MQ=255)
                                   gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGA.GGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                               >  2:5670334/1‑101 (MQ=255)
                                    aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                               >  1:3178166/1‑101 (MQ=255)
                                    aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                               >  1:1283618/1‑101 (MQ=255)
                                     a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                              <  1:1692353/101‑1 (MQ=255)
                                     a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                              <  2:1951910/101‑1 (MQ=255)
                                     a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                              >  1:3393478/1‑101 (MQ=255)
                                      gatggtg..ggCTTTCAACTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                                <  2:744295/95‑1 (MQ=25)
                                          ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                            <  2:4467326/101‑1 (MQ=255)
                                           gg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                                           <  2:284442/101‑1 (MQ=255)
                                            g..ggCTTTCATCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGt                                                          <  1:2292518/101‑1 (MQ=255)
                                               ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                                                         <  2:4237840/101‑1 (MQ=255)
                                               ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                                                         <  2:5398942/101‑1 (MQ=255)
                                                  tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                      >  1:5154420/1‑101 (MQ=255)
                                                  tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                      <  1:1366222/101‑1 (MQ=255)
                                                  tttCATCTTGCCAAAAATACCGGGGGGGGTAAA.GGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                      <  1:61000/101‑1 (MQ=37)
                                                   ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                                     >  2:4698769/1‑101 (MQ=255)
                                                    tCATCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                                    <  1:4357591/101‑1 (MQ=255)
                                                     cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGCCGTGccc                                                   >  1:566588/1‑101 (MQ=255)
                                                      aTCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                                  <  1:1243856/101‑1 (MQ=255)
                                                       tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTTCCCCt                                                 >  1:452330/1‑101 (MQ=255)
                                                       tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                                 <  2:2237122/101‑1 (MQ=255)
                                                       tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                                 >  1:4329083/1‑101 (MQ=255)
                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                                >  2:1313696/1‑101 (MQ=255)
                                                        cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                                >  1:1561900/1‑101 (MQ=255)
                                                         ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCGG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCGGCAACAGGTCGTGCCCCttt                                               >  2:2555882/1‑101 (MQ=37)
                                                         ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                               <  2:3581247/101‑1 (MQ=255)
                                                         ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                               >  2:3708063/1‑101 (MQ=255)
                                                         ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                               >  1:5744331/1‑101 (MQ=255)
                                                          tGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                              >  2:5059159/1‑101 (MQ=255)
                                                           caaaaaaaTCTGCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCGGCCC.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                                             <  2:1037156/97‑1 (MQ=255)
                                                              gAAAATGTCCGGGTGGGTGAA.TGGCGCA..GCCAGAGGGGGGCGG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                          <  1:3809560/101‑1 (MQ=255)
                                                              gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                          <  1:3739804/101‑1 (MQ=255)
                                                              gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                          <  2:374572/101‑1 (MQ=255)
                                                              gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                          >  1:269632/1‑101 (MQ=255)
                                                                aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                        >  1:4728132/1‑101 (MQ=255)
                                                                 aaTATCCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                                       <  1:1160701/101‑1 (MQ=38)
                                                                 aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                                       >  2:4098105/1‑101 (MQ=255)
                                                                  atatCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                                      <  1:4113666/101‑1 (MQ=255)
                                                                  atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGACGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                                      >  2:341770/1‑101 (MQ=255)
                                                                   tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                                     <  2:1875113/101‑1 (MQ=255)
                                                                   tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                                     <  1:1828752/101‑1 (MQ=255)
                                                                    atCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                                    <  1:885779/101‑1 (MQ=255)
                                                                        cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                                >  2:2127671/1‑101 (MQ=255)
                                                                         ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                               >  2:1095929/1‑101 (MQ=255)
                                                                         ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGCGCCCCTTTGTCATTCTGTCGCAcc                               >  1:2786991/1‑101 (MQ=255)
                                                                          gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                              >  2:2077836/1‑101 (MQ=255)
                                                                           ggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                             >  1:4904079/1‑101 (MQ=255)
                                                                           ggggGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                             <  2:3228511/101‑1 (MQ=255)
                                                                            gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                            >  1:1050514/1‑101 (MQ=255)
                                                                             ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                           <  2:2605921/101‑1 (MQ=255)
                                                                              gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                          >  2:3496758/1‑101 (MQ=255)
                                                                                 aa.TGGCGCA..GCCAGAGGGGGGCGG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                       >  2:2709452/1‑101 (MQ=255)
                                                                                 aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                       >  1:1424690/1‑101 (MQ=255)
                                                                                    tGGCGCA..GCCAGAGGGGGCCAG.CGCCCCCC.TCCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                     <  1:419555/101‑1 (MQ=255)
                                                                                     ggCGCA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg                    >  2:1006920/1‑101 (MQ=255)
                                                                                        gcA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg                 <  1:5599253/101‑1 (MQ=255)
                                                                                        gcA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg                 >  1:2018738/1‑101 (MQ=255)
                                                                                        gcA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg                 <  2:4537448/101‑1 (MQ=255)
                                                                                         cA..GCCAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg                <  1:3537913/101‑1 (MQ=255)
                                                                                             gcggggggggggCGG.CC.CCCCCTTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGcc              <  1:1976789/95‑1 (MQ=255)
                                                                                             atctgagGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGcc              <  2:4355326/97‑1 (MQ=255)
                                                                                               cAGAGGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc            >  2:3110152/1‑101 (MQ=255)
                                                                                                agagGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa           >  1:759998/1‑101 (MQ=255)
                                                                                                 gagGGGGGCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt          <  2:3778655/101‑1 (MQ=255)
                                                                                                   ggggggCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCg        <  2:3861415/101‑1 (MQ=255)
                                                                                                    gggggCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGa       >  1:2563139/1‑101 (MQ=255)
                                                                                                     ggggCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGAt      <  1:1885047/101‑1 (MQ=255)
                                                                                                     ggggCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTCTCGCACCCATGTTATAAGCTGGCCGCATCGAt      >  2:558408/1‑101 (MQ=255)
                                                                                                      gggCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATg     >  2:4203212/1‑101 (MQ=255)
                                                                                                        gCAG.CGCCCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTc   <  2:2994509/101‑1 (MQ=255)
                                                                                                               |                                                                                                
TGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCC  >  NC_014034/2249291‑2249486

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