Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,181,531	A→G	I183I (ATT→ATC)	gspE ←	general secretory pathway component, cryptic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,181,531	0	A	G	100.0%	127.7 / NA	66	I183I (ATT→ATC)	gspE	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/66);  total (0/66)

GCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTAATTCTTCCATCC  >  W3110S.gb/4181482‑4181543                                                  |             gCAGTGTGGATACGCTGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:53044/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:6934/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:5118/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:53888/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:54724/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:58313/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:61525/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:62350/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:62748/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:63895/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:64134/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:64421/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:6529/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:67890/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:67907/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:13230/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:75069/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:75090/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:77162/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:79696/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:79880/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:79933/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:80070/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:83918/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:84308/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:85393/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:85505/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:86234/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:87174/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:88574/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:48409/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:17063/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:17103/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:17344/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:27425/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:29237/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:31778/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:32958/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:33610/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:34164/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:35221/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:39672/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:40447/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:40481/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:41543/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:42082/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:48328/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:47436/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:46401/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:44357/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:44336/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:43938/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:43647/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:43043/62‑1 (MQ=255) gCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGACTCTTCCATcg  <  1:43329/62‑2 (MQ=255) gCAGTGTGGATACGCGGACATCCATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:80840/62‑1 (MQ=255) gCAGTGTGGATACCCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:28069/62‑1 (MQ=255)  cAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:63713/61‑1 (MQ=255)  cAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:5394/61‑1 (MQ=255)  cAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:17104/61‑1 (MQ=255)  cAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:48708/61‑1 (MQ=255)  cAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:11946/61‑1 (MQ=255)  cAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:9136/61‑1 (MQ=255)   aGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:47800/60‑1 (MQ=255)                     tCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:74604/42‑1 (MQ=255)                     tCTATGTTACGTCGCCCGATACGCAAACTGATTCTTCCATcc  <  1:65684/42‑1 (MQ=255)                                                  |             GCAGTGTGGATACGCGGACATCTATGTTACGTCGCCCGATACGCAAACTAATTCTTCCATCC  >  W3110S.gb/4181482‑4181543

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