Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,380 A→G 2.4% intergenic (‑176/‑129) 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,3800A→G 2.4% 13.4 ‑1.5 ‑1.3453intergenic (‑176/‑129)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (8/3):  ref base (226/216):  total (234/219)
Fisher's exact test strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.09e-02

Only 100 of 509 total aligned reads displayed.

GCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.CC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATC  >  NC_014034/2249280‑2249479
                                                                                                             |                                                                                                       
ccccgGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCG.AAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                        <  1:3658952/100‑1 (MQ=255)
cccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                                                                         >  1:7602815/2‑101 (MQ=255)
 ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                                        <  2:5535179/101‑1 (MQ=255)
  cccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                                                       <  1:1521022/101‑1 (MQ=255)
   ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggg                                                                                                      <  1:3531674/101‑1 (MQ=255)
     gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                                    >  2:9158144/1‑101 (MQ=255)
      tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                                   <  1:4053185/101‑1 (MQ=255)
       cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.c                                                                                                 >  1:1388815/1‑101 (MQ=255)
         cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg                                                                                               <  2:8502363/101‑1 (MQ=255)
           tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGG.CAgcg                                                                                             >  2:3533680/1‑101 (MQ=255)
           caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGG.c                                                                                                 <  2:6110797/96‑1 (MQ=18)
           caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGG.c                                                                                                 >  2:6516338/6‑101 (MQ=18)
            aagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggggg                                                                                                 >  2:7429631/5‑99 (MQ=18)
              cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c                                                                                        <  1:7664704/101‑1 (MQ=255)
               gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.c                                                                                      <  2:1435021/101‑1 (MQ=255)
                 gACT.CCCCCGGGATATTTGACGCAAGATGAN...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.ccc                                                                                    <  1:7209702/101‑2 (MQ=255)
                  actccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGG.CAGCG.cc.c                                                                                        <  1:1177201/97‑1 (MQ=17)
                    t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.cccttc                                                                                 >  1:94288/1‑97 (MQ=255)
                    t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.cccttc                                                                                 <  2:108225/101‑5 (MQ=255)
                    t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTGGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.cccttc                                                                                 <  2:1735425/101‑5 (MQ=255)
                       ccccGGGATATTTGAGGCAAGATGAA...GGTG..TCCTTTCATCTTGCCGAAAATATCCCAGGCGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTcgc                                                                                <  2:7325280/101‑1 (MQ=21)
                       ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTcgc                                                                                <  1:8837916/101‑1 (MQ=255)
                        cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGG.CAGCG.CC.C.CCtgc                                                                                  >  2:8121891/1‑99 (MQ=2)
                        cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTCTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGTGGGG.CAGCGCCC.C.CCTTcgcg                                                                               >  1:6446734/1‑101 (MQ=25)
                         ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCgcgg                                                                              >  1:7989119/1‑101 (MQ=255)
                              atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGggcggc                                                                         >  2:3121539/1‑101 (MQ=255)
                               tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGcc                                                                        >  2:3705213/1‑101 (MQ=255)
                               tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGcc                                                                        >  1:201674/1‑101 (MQ=255)
                                atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCcg                                                                       >  2:4167831/1‑101 (MQ=255)
                                   tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGa                                                                    <  1:7902126/101‑1 (MQ=255)
                                   tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGa                                                                    >  2:6930111/1‑101 (MQ=255)
                                    gAGGCAAGATGAA...GGGG..GTCTTTCATCTTGCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAc                                                                   <  2:9023061/101‑1 (MQ=37)
                                    gAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAc                                                                   <  2:5682043/101‑1 (MQ=255)
                                    gAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAc                                                                   <  2:2203705/101‑1 (MQ=255)
                                    gAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAc                                                                   >  2:4908152/1‑101 (MQ=255)
                                       gCAAGATGAA...GGGG..GTCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAcgcg                                                                <  2:2026013/101‑1 (MQ=255)
                                        cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAcgcgc                                                               >  1:3376280/1‑101 (MQ=255)
                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGGcgcgcg                                                              >  2:6260781/1‑101 (MQ=255)
                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGCcgcgcg                                                              >  1:2233317/1‑101 (MQ=255)
                                            aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGcgccg                                                           >  2:5111174/1‑101 (MQ=255)
                                              gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                         <  2:8915125/101‑1 (MQ=255)
                                              gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                         >  1:6942826/1‑101 (MQ=255)
                                              gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                         >  2:7758247/1‑101 (MQ=255)
                                               aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                        >  2:8296706/1‑101 (MQ=255)
                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                       <  2:4716052/101‑1 (MQ=255)
                                                     ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                     >  1:2566544/1‑101 (MQ=255)
                                                     ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                     <  2:2704485/101‑1 (MQ=255)
                                                     ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                     <  2:7194864/101‑1 (MQ=255)
                                                       g..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAggc                                                   >  1:5012030/1‑100 (MQ=255)
                                                       g..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAAGCGCGCCGCAACAggc                                                   >  1:8282087/1‑100 (MQ=255)
                                                          ggCGTTCATCTGCCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCGAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                                                  <  2:1852102/101‑1 (MQ=37)
                                                            cTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                                <  2:2398581/101‑1 (MQ=255)
                                                             tttCATGTTGCCGAAAGTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                               <  2:7624217/101‑1 (MQ=37)
                                                             tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                               >  2:46374/1‑101 (MQ=255)
                                                             tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                               <  2:2085355/101‑1 (MQ=255)
                                                              ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                              >  1:470451/1‑101 (MQ=255)
                                                                 aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                           >  2:5970751/1‑101 (MQ=255)
                                                                  tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGGCGTGCCCCt                                          >  1:6734588/1‑101 (MQ=255)
                                                                   cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                         >  2:1207273/1‑101 (MQ=255)
                                                                   cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                         >  2:3901468/1‑101 (MQ=255)
                                                                    ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                        <  1:9108057/101‑1 (MQ=255)
                                                                      gCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                                      <  1:6010541/101‑1 (MQ=255)
                                                                       ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCTTGCCCCTTTGTc                                     <  1:6630177/101‑1 (MQ=255)
                                                                       ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGCGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                     >  2:545694/1‑101 (MQ=255)
                                                                        cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                    >  2:6984678/1‑101 (MQ=255)
                                                                          aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                  <  1:7425348/101‑1 (MQ=255)
                                                                            aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGAGGGCG.CC.C.CTTCCGCGGG.GGGCGCGACGCGCGCCGCAACAGGTCCGGCCCCTTTGTCAttttt                               >  2:4449441/1‑98 (MQ=25)
                                                                            aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                                <  1:168561/101‑1 (MQ=255)
                                                                             atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.CC.CACCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGGGCCCCTTTGTCATTGTg                               >  1:4225264/1‑101 (MQ=38)
                                                                             atatCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGACGTGCCCCTTTGTCATTCTg                               >  2:1300017/1‑101 (MQ=38)
                                                                             atagCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                               <  1:7729471/101‑1 (MQ=255)
                                                                               atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                             <  2:2864885/101‑1 (MQ=255)
                                                                               atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                             >  2:4768171/1‑101 (MQ=255)
                                                                                 cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGGCGc                           >  1:2953672/1‑101 (MQ=255)
                                                                                    ggggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                        <  1:3820811/101‑1 (MQ=255)
                                                                                    ggggggTGAA.GGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                        <  1:3222386/101‑1 (MQ=255)
                                                                                     gggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                       >  1:5520198/1‑101 (MQ=255)
                                                                                       gtgTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                     <  2:5871198/99‑1 (MQ=255)
                                                                                        ggTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                    <  2:1523736/101‑1 (MQ=255)
                                                                                         gTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                   <  1:8624604/101‑1 (MQ=255)
                                                                                         gTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                   <  2:6665599/101‑1 (MQ=255)
                                                                                          ttggactgcgcA..TCGAGCGGGAGG.GAGCG.CCGC.TCGTCGGTGGGGGCCGGGAGGCGCGCCGCAACAGGTCGTTCCCCTTTTTCAATCTGTCGCGCCCATGt                   <  2:7961388/94‑1 (MQ=255)
                                                                                           gAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                 >  1:4823683/1‑101 (MQ=255)
                                                                                            atctGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                 <  2:5295001/98‑1 (MQ=255)
                                                                                            aa.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                <  1:5397756/101‑1 (MQ=255)
                                                                                             a.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa               <  2:6597307/101‑1 (MQ=255)
                                                                                             a.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa               >  1:5697934/1‑101 (MQ=255)
                                                                                              agggCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa               <  1:6240183/99‑1 (MQ=255)
                                                                                               tGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa              >  1:4361087/1‑101 (MQ=255)
                                                                                               tGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa              <  1:6282056/101‑1 (MQ=255)
                                                                                                ggCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCaagat             >  1:5295185/1‑97 (MQ=255)
                                                                                                  cgcA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt           <  2:2604963/101‑1 (MQ=255)
                                                                                                   gcA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAgaga          >  1:2346513/1‑98 (MQ=255)
                                                                                                        gCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGcc       <  1:3142289/101‑1 (MQ=255)
                                                                                                           agagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa    >  2:8654391/1‑101 (MQ=255)
                                                                                                           agagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa    >  1:1882767/1‑101 (MQ=255)
                                                                                                           agagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa    >  1:8482025/1‑101 (MQ=255)
                                                                                                            gagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt   >  1:2849322/1‑101 (MQ=255)
                                                                                                             agGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATc  >  1:1024218/1‑101 (MQ=255)
                                                                                                             agGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATc  <  2:3333783/101‑1 (MQ=255)
                                                                                                             |                                                                                                       
GCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.CC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATC  >  NC_014034/2249280‑2249479

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG