Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,749,892	A→G	S199S (TCT→TCC)	tnaA ←	tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,749,892	0	A	G	100.0%	127.8 / NA	64	S199S (TCT→TCC)	tnaA	tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/64);  total (0/64)

TTCGCGATGCTGTACATCGCTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCAGAGT  >  W3110S.gb/3749836‑3749896                                                         |     ttcGCGATGCTGTACATCACTTTTACGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:4979/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTTAAACCGGCTGACCACCTGCGGAGt  <  1:3979/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:49445/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:41770/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:10862/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:42795/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:44246/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:44317/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:45149/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:46104/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:46755/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:47145/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:47659/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:48186/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:48725/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:49120/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:49313/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:8458/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:8566/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:49853/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:49908/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:50169/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:50982/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:52819/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:53368/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:5362/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:54005/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:54482/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:6422/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:6843/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:8066/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:27273/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:12570/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:14134/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:14412/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:14765/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:15985/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:16465/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:19224/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:19889/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:20183/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:20428/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:20546/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:20577/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:22949/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:23074/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:24620/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:34034/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:39083/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:38242/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:36848/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:35372/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:34136/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:34079/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:40298/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:33903/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:32782/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:32562/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:32117/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:29642/61‑1 (MQ=255) ttcGCGATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:27641/61‑1 (MQ=255)      gATGCTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:43983/56‑1 (MQ=255)          cTGTACATCACTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:42297/52‑1 (MQ=39)                        tAAGTTTGCCAGTGAAACCGGCTGACCACCTGCGGAGt  <  1:50577/38‑1 (MQ=255)                                                         |     TTCGCGATGCTGTACATCGCTTTTAAGTTTGCCAGTGAAACCGGCTGACCACCTGCAGAGT  >  W3110S.gb/3749836‑3749896

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