Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	486,827	G→C	L356L (CTG→CTC)	kefA →	fused mechanosensitive channel proteins

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	486,827	0	G	C	84.5%	151.2 / 20.1	58	L356L (CTG→CTC)	kefA	fused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base G (0/9);  new base C (49/0);  total (49/9)
Fisher's exact test for biased strand distribution p-value = 9.39e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGA  >  W3110S.gb/486786‑486855                                          |                             ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGa                >  1:1855020/1‑56 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2425898/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2089688/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2095060/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:587652/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2161700/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2175270/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:542518/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2188077/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2190367/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2205327/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2278927/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2385515/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:665059/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2601543/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2724623/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:508112/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2759978/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:503862/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:2884081/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:364425/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:386113/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:392385/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:420020/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1871057/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1184895/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1188783/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1213648/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1217898/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1227499/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:948094/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:911840/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:904239/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1798240/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:845169/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1882546/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1894827/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1910216/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1927762/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1937820/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1946007/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:195024/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1963184/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1965338/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTgg           >  1:1997249/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTg            >  1:976194/1‑60 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTg            >  1:1867712/1‑60 (MQ=255) ctgctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTAGGCAGATGAACTgg           >  1:540106/1‑61 (MQ=255) ctgctgTTGTCACGTATCCTCTACCAGCAACAACAAACACTCCCCTCGGCAGATGAAc              >  1:20110/1‑58 (MQ=255)         gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:2787516/62‑1 (MQ=255)         gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:2187547/62‑1 (MQ=255)         gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:21047/62‑1 (MQ=255)         gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:69363/62‑1 (MQ=255)         gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:1564645/62‑1 (MQ=255)         gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:1556384/62‑1 (MQ=255)         gtCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:100084/62‑1 (MQ=255)          tctcGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:2742204/61‑1 (MQ=255)          tctcGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGa  <  1:1261116/61‑1 (MQ=255)                                          |                             CTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGA  >  W3110S.gb/486786‑486855

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