Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,961,422	(T)5→6	coding (412/864 nt)	yjjW ←	predicted pyruvate formate lyase activating enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,961,417	1	.	T	100.0%	59.5 / NA	20	N139K (AAC→AAA)	yjjW	predicted pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (0/20);  total (0/20)

GTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑TTTTTGATAGCAGTAAACAGCGCCAC  >  minE/2961358‑2961443                                                             |                           gTTTCGCTCAACATGCCGTTACTGTCCACCAGCCAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:303662/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:1215217/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:69557/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:3174770/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:2781856/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:2757833/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:2754594/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:2686475/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:2435777/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:2272114/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:2039419/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:1771665/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:1540693/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:1511226/1‑62 (MQ=255) gTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑tt                          >  1:1270578/1‑62 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:1358058/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:705544/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:1097256/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:643576/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:437524/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:356258/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:322615/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:1242156/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:2912325/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:2889741/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:279028/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:1813039/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:1383664/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:1514585/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:266618/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:1518779/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:236754/61‑1 (MQ=255)                           caccaGACAGGTGAGATGGAGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:2420327/61‑1 (MQ=255)                            accaGACAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:1958866/60‑1 (MQ=255)                                 aCAGGTGAGATGGCGCAGTTGCGGATCGTTTTTTGATAGCAGTAAACAGCGccac  <  1:2646845/55‑1 (MQ=255)                                                             |                           GTTTCGCTCAACATGCCGTTACTGTCCACCAGACAGGTGAGATGGCGCAGTTGCGGATCG‑TTTTTGATAGCAGTAAACAGCGCCAC  >  minE/2961358‑2961443

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