Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,132,425	G→T	L63L (CTC→CTA)	yghS ←	hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,132,425	0	G	T	86.7%	28.0 / 0.9	15	L63L (CTC→CTA)	yghS	hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (0/13);  total (2/13)
Fisher's exact test for biased strand distribution p-value = 9.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGTGCACATGTGCCGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGC  >  W3110S.gb/3132364‑3132464                                                              |                                        tttCGTGCACATGTGCCGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAggg                                         >  1:1844258/1‑62 (MQ=255) tttCGTGCACATGTGCCGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAggg                                         >  1:485064/1‑62 (MQ=255)                 cGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCAtt                         >  1:1645373/1‑62 (MQ=255)                 cGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCAtt                         >  1:5281/1‑62 (MQ=255)                                        cAAAAGGTGCGCCAATAACAGGTAGCTGTGAATTCCATTCGCCAATTCGCCCGGACGATTGc  <  1:946678/62‑1 (MQ=255)                                        cAAAAGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:2488462/62‑1 (MQ=255)                                        cAAAAGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:1057694/62‑1 (MQ=255)                                        cAAAAGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:390000/62‑1 (MQ=255)                                        cAAAAGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:3060081/62‑1 (MQ=255)                                        cAAAAGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:351824/62‑1 (MQ=255)                                         aaaaGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:2500803/61‑1 (MQ=255)                                         aaaaGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:720890/61‑1 (MQ=255)                                         aaaaGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:2588104/61‑1 (MQ=255)                                         aaaaGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:239312/61‑1 (MQ=255)                                         aaaaGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:1742724/61‑1 (MQ=255)                                         aaaaGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:1733177/61‑1 (MQ=255)                                         aaaaGGTGCGCCAATAACAGGTAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGc  <  1:1565827/61‑1 (MQ=255)                                                              |                                        TTTCGTGCACATGTGCCGCTTTACTTCGCAGATAACGCCCAAAAGGTGCGCCAATAACAGGGAGCTGTGAAATCCATTCGCCAATTCGCCCGGACGATTGC  >  W3110S.gb/3132364‑3132464

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