Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 802,452 +T coding (601/1071 nt) dhaK ← dihydroxyacetone kinase, N‑terminal domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE802,4491.T96.7% 220.0 / 0.7 61coding (604/1071 nt)dhaKdihydroxyacetone kinase, N‑terminal domain
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (59/0);  total (61/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

AATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAA  >  minE/802402‑802466
                                                |                  
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCATGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:902515/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCggaaggaa  >  1:700260/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:6152/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:220249/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1000382/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:291459/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:334642/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:354514/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:378350/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:461311/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:474994/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:503564/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:539767/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:562469/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:593472/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:603812/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:606889/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:26164/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:641225/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:656285/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:670094/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:671851/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:694746/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:796438/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:81139/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:817389/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:88186/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:922871/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:940471/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:979046/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1549302/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1007481/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1010520/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1021524/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1026811/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1057571/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1067664/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1096177/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1148061/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:115263/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1153364/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1157701/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1168104/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1216758/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1224794/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1236886/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1291829/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1297200/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:130246/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1343707/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1356124/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1392780/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1428251/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1477536/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1537159/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1539884/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:1546433/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGaa  >  1:156985/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTTGCCCGCGGCAGGa   >  1:608164/1‑65 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGt                   >  1:476266/1‑49 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTGTTGCCCGCGGCAGGaa  >  1:1431519/1‑66 (MQ=255)
aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTCTTGCCCGCGGCAGGaa  >  1:430652/1‑66 (MQ=255)
                                                |                  
AATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTTTGCCCGCGGCAGGAA  >  minE/802402‑802466

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: