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
*	minE	609,683	0	C	.	94.1%	128.9 / 2.7	34	coding (571/2283 nt)	pflB	pyruvate 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                       |                                                              aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCGCCGATGATacggcc                  >  1:354339/1‑67 (MQ=39) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCg                             >  1:806363/1‑56 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:27488/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:861464/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:816810/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:815281/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:796189/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:776574/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:707342/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:514137/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:481777/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:436699/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:395169/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:295751/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:295158/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1175347/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:22578/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1410146/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1475555/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:156287/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:167895/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1711883/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:1899162/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:190763/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:20188/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:203585/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:2052302/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:2121431/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  1:233154/1‑67 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggc                   >  1:1328119/1‑66 (MQ=255) aGTTTGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggc                   >  1:1272578/1‑66 (MQ=255) aGTTAGTCTTTCATCAGGTAGT‑GATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                  >  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: ‑