Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,683 Δ1 bp 100% coding (571/2283 nt) pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,6830C.94.1% 128.9 / 2.7 34coding (571/2283 nt)pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base . (32/0);  total (34/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.67e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCT  >  minE/609661‑609744
                      |                                                             
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCGCCGATGATacggcc                  >  1:354339/1‑67 (MQ=39)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCg                             >  1:806363/1‑56 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:27488/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:861464/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:816810/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:815281/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:796189/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:776574/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:707342/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:514137/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:481777/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:436699/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:395169/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:295751/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:295158/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1175347/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:22578/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1410146/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1475555/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:156287/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:167895/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1711883/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1899162/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:190763/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:20188/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:203585/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:2052302/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:2121431/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:233154/1‑67 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggc                   >  1:1328119/1‑66 (MQ=255)
aGTTTGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggc                   >  1:1272578/1‑66 (MQ=255)
aGTTAGTCTTTCATCAGGTAGTGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:278621/1‑67 (MQ=255)
                 gTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc   >  1:1623122/1‑66 (MQ=255)
                 gTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt  >  1:768436/1‑67 (MQ=255)
                      |                                                             
AGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCT  >  minE/609661‑609744

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: