Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,067,769	A→G	E285E (GAA→GAG)	glgB →	1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,067,769	0	A	G	100.0%	138.1 / NA	49	E285E (GAA→GAG)	glgB	1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (49/0);  total (49/0)

TGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCACCTCGAACTACTGCCC  >  W3110S.gb/4067717‑4067778                                                     |          tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTCCCCATCTCGAGCTACTGcct  >  1:938106/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:755520/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:44840/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1115811/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:57736/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:608633/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:616969/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:620172/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:649309/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:687641/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:692586/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:731687/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:744063/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:482665/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:824433/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:844248/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:873821/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:875366/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:915089/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:943487/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:970794/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:977679/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:99266/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:15693/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1126960/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1136525/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1187747/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1205612/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1325772/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1346233/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1368445/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1404411/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:1444498/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:385633/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:172291/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:173477/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:220594/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:224055/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:279151/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:352448/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:35645/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:356639/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:360430/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcc   >  1:203251/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcc   >  1:648280/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcc   >  1:892450/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcc   >  1:516019/1‑61 (MQ=255) tGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcat  >  1:1408946/1‑60 (MQ=255) tGGCAGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCATCTCGAGCTACTGcct  >  1:973388/1‑61 (MQ=255)                                                     |          TGGCCGATCAACTGGTGCCTTATGCTAAATGGATGGGCTTTACCCACCTCGAACTACTGCCC  >  W3110S.gb/4067717‑4067778

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