Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,614,294	T→A	P224P (CCA→CCT)	yneH ←	predicted glutaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,614,294	0	T	A	92.4%	204.0 / 6.2	66	P224P (CCA→CCT)	yneH	predicted glutaminase
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base A (0/61);  total (5/61)
Fisher's exact test for biased strand distribution p-value = 1.12e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTC  >  W3110S.gb/1614257‑1614318                                      |                         gcgcGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTc  >  1:1056047/1‑62 (MQ=255) gcgcGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTc  >  1:1116073/1‑62 (MQ=255) gcgcGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTc  >  1:1266851/1‑62 (MQ=255) gcgcGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTc  >  1:2772693/1‑62 (MQ=255) gcgcGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTc  >  1:494451/1‑62 (MQ=255)            tGCCGCGGCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:126302/43‑1 (MQ=38)            tGCCGCGCCTGCATTGTCGTCACCACAGGTTCATCAATATGaa          <  1:2903452/43‑1 (MQ=38)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATGAATATGaa          <  1:1128231/43‑1 (MQ=38)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2561038/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2821506/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2799402/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2618809/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:258142/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2580031/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2567934/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2566906/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:367012/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2446330/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2431401/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2402984/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2317768/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2893972/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:327173/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:217621/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:389862/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:395043/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:445798/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:571836/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:633931/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:681944/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:693732/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:70975/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:875868/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:87861/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:926071/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1705365/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1092790/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1112614/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1266100/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1367409/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1370835/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1435160/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1438987/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1460312/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1515839/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1521079/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1532539/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1546034/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1666674/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1704455/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2292624/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1713152/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1751296/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:178709/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1836268/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1908121/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1937224/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:1943089/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2023242/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2089514/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2135765/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2145368/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:2154844/43‑1 (MQ=255)            tGCCGCGCCTGCATTGGCGTCACCACAGGATCATCAATATGaa          <  1:1858569/43‑1 (MQ=38)            agCCGCGCCTGCATTGGCGTCACCACAGGTTCATCAATATGaa          <  1:782504/42‑1 (MQ=38)                    cTGCATTGGCGTCACCACAGGATCATCAATATGaa          <  1:2743028/35‑1 (MQ=37)                                      |                         GCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTC  >  W3110S.gb/1614257‑1614318

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