Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	539,432	T→A	G354G (GGT→GGA)	allB →	allantoinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	539,432	0	T	A	88.6%	125.4 / 5.8	44	G354G (GGT→GGA)	allB	allantoinase
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base A (39/0);  total (44/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAA  >  W3110S.gb/539404‑539470                             |                                       tGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCaa                 >  1:2052158/1‑52 (MQ=255) tGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGc        >  1:951395/1‑61 (MQ=255) tGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCg       >  1:954553/1‑62 (MQ=255) tGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCg       >  1:1505566/1‑62 (MQ=255) tGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCg       >  1:2865117/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTaa            >  1:1722500/1‑52 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:762268/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2531064/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2606424/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2650881/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:279777/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2864712/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2922790/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:338476/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:404314/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:505561/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2168525/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:813434/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:817013/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:91599/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:916127/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:954939/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:99620/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1845549/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:127559/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1300893/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1303504/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1431488/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1465496/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1482498/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1547784/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1661245/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1816269/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2509094/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1867561/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:188260/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1973238/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1978631/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2041822/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2077727/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2114706/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:2150278/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTaa  >  1:1149681/1‑62 (MQ=255)      gATGAAGCGGTACAGAAACGCGGAATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTa   >  1:530763/1‑61 (MQ=255)                             |                                       TGTTCGATGAAGCGGTACAGAAACGCGGTATGTCTCTGCCAATGTTCGGCAAATTAATGGCGACTAA  >  W3110S.gb/539404‑539470

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