Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,379 G→A 2.7% intergenic (‑175/‑130) parE ← / → RCAP_rcc02080 DNA topoisomerase IV subunit B/HAD superfamily hydrolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,249,3790G→A 2.7% 2.6 ‑1.0 ‑0.0165intergenic (‑175/‑130)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (80/80):  total (84/81)
Fisher's exact test strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.63e-01

Only 100 of 192 total aligned reads displayed.

TGCCCCGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCA  >  NC_014034/2249279‑2249477
                                                                                                              |                                                                                                       
gcccccGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCag                                                                                                         <  2:1149245/99‑1 (MQ=255)
gcccccGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCag                                                                                                         <  1:2604459/99‑1 (MQ=255)
gcccccGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGC.GCA..GCCag                                                                                                         >  1:2243123/3‑101 (MQ=255)
  ccccGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..CCCagag                                                                                                       >  2:1681540/1‑101 (MQ=255)
     cgtcccgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAg.ggg                                                                                                   >  2:3683041/1‑101 (MQ=255)
        ccgatctcGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGc                                                                                                <  2:425147/94‑1 (MQ=255)
        cccgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGc                                                                                                >  2:3872364/1‑101 (MQ=255)
          cgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.g                                                                                            >  1:84139/1‑101 (MQ=255)
          cgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.g                                                                                            >  2:4102191/1‑101 (MQ=255)
          cgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.g                                                                                            <  2:24645/101‑1 (MQ=255)
           gtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.g.c                                                                                          <  1:1001720/101‑1 (MQ=255)
               cGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.ccc                                                                                     >  1:1918090/1‑101 (MQ=255)
               cGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.ccc                                                                                     >  2:2919430/1‑101 (MQ=255)
               cGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.ccc                                                                                     >  1:2850463/1‑101 (MQ=255)
                gggctccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGC.GAA..GCCAGAG.GGGGGC.A.g.cg.c                                                                                       <  1:3059470/95‑1 (MQ=17)
                gggACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.cccc                                                                                    >  1:354411/1‑101 (MQ=255)
                 ggACTCCCCCGGGATATTTGAGGCAAGATGAC...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.ccccc                                                                                   >  1:3808942/1‑101 (MQ=255)
                  gACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.cccccc                                                                                  <  2:3593020/101‑2 (MQ=255)
                   aCTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.TGGGGC.A.G.CG.cccccct                                                                                 <  2:4129505/101‑3 (MQ=255)
                     tCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.ccccccttc                                                                               >  1:1710117/1‑97 (MQ=255)
                     ccccccGGGTTTTTTGGGGGAAGAGGAA...GGGG..GGTTTTCATTTTGCCAAAAAAATCCCGGGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.cccccctt                                                                                <  1:1422384/100‑4 (MQ=21)
                      cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.ccccccttcg                                                                              <  2:3959378/101‑6 (MQ=255)
                      cccccGGGATATTTGAGGCAAGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.ccccccttcg                                                                              <  2:450745/101‑6 (MQ=255)
                       ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTcgc                                                                              <  1:1547977/101‑1 (MQ=255)
                       ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTcgc                                                                              <  2:1633285/101‑1 (MQ=255)
                       ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTcgc                                                                              >  1:2354331/1‑101 (MQ=255)
                        cccGGGATATTTGCAGCAAGGTGAAACCGGTT..GCCTTTCATCTGGCCCGAAATATCCC.GGGGGGTGAA.TGGC.GAA..GCCAGAG.GGGGGC.A.G.CG.ccctcctgc                                                                               <  1:291216/101‑7 (MQ=14)
                         ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCgcgg                                                                            >  1:82437/1‑101 (MQ=255)
                          cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGAGAAT.TGGC.CAAAGGCCAAGA.GGGGGC.A.G.CG.CCCCCtgacc                                                                              >  2:1879142/1‑97 (MQ=2)
                          cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC..AGGGGTGAATTGGC.CAAAGGCCAAGA.GGGGGC.A.G.CG.CCCCCtgccc                                                                              >  2:1663196/1‑97 (MQ=0)
                          cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC..AGGGGTGAATTGGC.CAAAGGCCAAGA.GGGGGC.A.G.CG.CCCCCtgccc                                                                              >  2:1996119/1‑97 (MQ=2)
                          cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCgcggg                                                                           >  1:636164/1‑101 (MQ=255)
                          cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCgcggg                                                                           >  1:433227/1‑101 (MQ=255)
                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCG.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCgcgggc                                                                          <  2:2916758/101‑1 (MQ=255)
                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCgcgggc                                                                          >  1:1763632/1‑101 (MQ=255)
                               tatTTGCGGAAAGAGGAA...GGGGCATGTTTTCATCTTGCCCAAAATATCCC..GGGGGTGAATTGGCAAAA..GGCAAAGAGGGGGC.A.G.CG.CCCCCtgccccgccg                                                                         <  2:3139117/101‑10 (MQ=9)
                               tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGcc                                                                      <  1:1141657/101‑1 (MQ=255)
                                atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCcg                                                                     >  1:3698587/1‑101 (MQ=255)
                                atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTCCGCGGGCGGCcg                                                                     >  1:991321/1‑101 (MQ=255)
                                       gCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcg                                                              <  2:1112260/101‑1 (MQ=255)
                                         aaGATGAG...GGGG..GGCTTTCATCTTGCCGAAAATAGCGC.GGGGGGTGAA.TGGC.GCA..GCGAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                                                            <  2:1066977/101‑1 (MQ=37)
                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                                                            >  2:1837501/1‑101 (MQ=255)
                                          aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACgcgcgc                                                           >  2:1588134/1‑101 (MQ=255)
                                          aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACgcgcgc                                                           >  1:4232831/1‑101 (MQ=255)
                                             tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccgc                                                        >  2:3799270/1‑101 (MQ=255)
                                              gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                       >  1:735256/1‑101 (MQ=255)
                                              gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                       <  2:3521903/101‑1 (MQ=255)
                                              gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                       <  1:1461113/101‑1 (MQ=255)
                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                     <  1:1899413/101‑1 (MQ=255)
                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                     <  2:2262820/101‑1 (MQ=255)
                                                    gggg..ggTTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                    <  2:899651/101‑1 (MQ=255)
                                                    gggg..ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                    >  2:1269296/1‑101 (MQ=255)
                                                          ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                                                >  2:2568943/1‑101 (MQ=255)
                                                           gCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                                               <  1:3207496/101‑1 (MQ=255)
                                                                cATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                          <  1:2527817/101‑1 (MQ=255)
                                                                cATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                          <  2:3179523/101‑1 (MQ=255)
                                                                  tCTTGCCGAAAATATCCG.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                        <  1:2631767/101‑1 (MQ=255)
                                                                  tCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                        <  2:2302495/101‑1 (MQ=255)
                                                                   cTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                       <  1:3529238/101‑1 (MQ=255)
                                                                   cTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGCCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                       >  1:2265057/1‑101 (MQ=255)
                                                                    ttGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                      <  1:1265657/101‑1 (MQ=255)
                                                                     tggcGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                     >  1:2420108/4‑101 (MQ=255)
                                                                     tGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                     >  1:1651336/1‑101 (MQ=255)
                                                                     tGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                     >  2:2885210/1‑101 (MQ=255)
                                                                     ggcagAAAATACCCG.GGGGGGTGAA.GGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                     <  2:811511/97‑1 (MQ=37)
                                                                      gCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                                    >  2:3357811/1‑101 (MQ=255)
                                                                       ccGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                   >  1:3153813/1‑101 (MQ=255)
                                                                        cGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                  <  2:1009263/101‑1 (MQ=255)
                                                                          aaaaTATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGGGCCCCTTTGTCAtt                                >  1:1053140/1‑101 (MQ=255)
                                                                           aaaTATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGGG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                               <  1:3785762/101‑1 (MQ=255)
                                                                             aatgtcggggggggTGAA.TGGC.GCA..GCCAGGG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                              <  2:2978343/93‑1 (MQ=38)
                                                                               ttCCG.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                           <  1:3364653/100‑1 (MQ=255)
                                                                               atCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGCCGTGCCCCTTTGTCATTCTGTc                           >  1:2452392/1‑101 (MQ=255)
                                                                                tcgc.GGGGGGGGAAGGGGC.GCA..GCCGGGG.GGGGGCGC.T.CC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                           <  1:2944405/98‑1 (MQ=255)
                                                                                gccc.GGGGGGCGAG.GGGG.GGA..GCCGGGG.GGGGGC.A.T.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCg                          <  1:1123182/100‑1 (MQ=255)
                                                                                  cc.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGGGCCCCTTTGTCATTCTGTCGCa                        >  1:3286465/1‑101 (MQ=255)
                                                                                   c.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                       <  2:2787310/101‑1 (MQ=255)
                                                                                     ggggggTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                      <  1:208656/101‑1 (MQ=255)
                                                                                     ggggggTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                      >  2:4010255/1‑101 (MQ=255)
                                                                                      gggggTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                     <  2:1917964/101‑1 (MQ=255)
                                                                                       ggggTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                    >  2:3715464/1‑101 (MQ=255)
                                                                                       ggggTGAA.TGGC.GCA..GCCAGAG.GGGGCC.ACC.CC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                    <  2:514208/101‑1 (MQ=255)
                                                                                        gggTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                   <  2:982947/101‑1 (MQ=255)
                                                                                        gggTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                   <  2:1706323/101‑1 (MQ=255)
                                                                                         ggTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                  <  2:3299626/101‑1 (MQ=255)
                                                                                          gTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                 <  2:1777126/101‑1 (MQ=255)
                                                                                           tGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc                >  2:1436569/1‑101 (MQ=255)
                                                                                            gAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGACGCGCCGCGCGCCGCAACAGGCAGTGCCCCTTTGTCATTCTGTCGTCCCCCTGTCa               >  2:1900639/1‑101 (MQ=255)
                                                                                            gAA.GGGC.GCA..GCCAGGG.GGGGGC.A.GCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa               <  1:3683295/101‑1 (MQ=255)
                                                                                             ttccgatctcA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc                <  2:4117340/92‑1 (MQ=255)
                                                                                             aa.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt              >  1:779448/1‑101 (MQ=255)
                                                                                              a.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa             >  1:1882425/1‑101 (MQ=255)
                                                                                                tGGC.GCA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa            <  2:289922/101‑1 (MQ=255)
                                                                                                   c.gcA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt         >  2:3847444/1‑101 (MQ=255)
                                                                                                     gcA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg        <  1:1251423/101‑1 (MQ=255)
                                                                                                      cA..GCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCagatcggaa       >  1:4117601/1‑93 (MQ=255)
                                                                                                          gCCAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGcc     <  2:1289375/101‑1 (MQ=255)
                                                                                                            cAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc   >  1:3135784/1‑101 (MQ=255)
                                                                                                            cAGAG.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGCCATTCTGTCGCACCCATGTCATAAGCTGGCCGc   >  2:2037613/1‑101 (MQ=255)
                                                                                                             agag.GGGGGC.A.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa  >  1:3947271/1‑101 (MQ=255)
                                                                                                              |                                                                                                       
TGCCCCGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGC.GCA..GCCAGAG.GGGGGC.A.G.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCA  >  NC_014034/2249279‑2249477

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG