Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,918,445	A→G	R158R (CGA→CGG)	yebS →	conserved inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,918,445	0	A	G	100.0%	175.3 / NA	63	R158R (CGA→CGG)	yebS	conserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/63);  total (0/63)

TTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGG  >  W3110S.gb/1918416‑1918457                              |             ttACGCCCGGTACTTCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:331460/42‑1 (MQ=25) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:732221/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:469650/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:498584/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:512415/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:515031/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:530188/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:531532/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:543696/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:596467/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:599016/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:649266/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:667452/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:680796/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:723289/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:728702/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:424534/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:773122/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:789728/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:861533/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:863648/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:864812/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:906833/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:91399/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:933173/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:935068/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:936722/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:967031/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:96971/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:973407/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:982830/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:42030/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1004527/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1004986/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1057726/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1059398/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1092168/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1110391/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1132027/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1151829/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:116110/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1244024/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1256649/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1262937/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1266376/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1306491/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:191402/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1003747/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:383175/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:366551/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:308285/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:297970/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:287648/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:225552/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1307511/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1389780/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1339618/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1338793/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:133829/42‑1 (MQ=37) ttACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:13133/42‑1 (MQ=37)  tACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:1371034/41‑1 (MQ=37)  tACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:953518/41‑1 (MQ=37)  tACGCCCGGTACTGCTGATGCTTGAGCGGCTGAAAGAGTgg  <  1:748444/41‑1 (MQ=37)                              |             TTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGG  >  W3110S.gb/1918416‑1918457

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: ‑