Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,637,667	+A	coding (273/621 nt)	yfgM ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,637,664	1	.	A	76.3%	161.5 / 62.5	76	coding (276/621 nt)	yfgM	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/18);  new base A (0/58);  total (0/76)
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

CAAGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTCCGCCG  >  W3110S.gb/2637657‑2637719         |                                                        cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:757018/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:233297/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:242180/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:31731/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:33988/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:410859/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:559547/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:57629/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:593560/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:628043/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:638767/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:705963/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:729105/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:734891/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:743106/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:755629/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:103909/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:782192/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:792906/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:807727/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:809714/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:836942/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:844652/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:84644/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:892081/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:920875/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:952779/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:953623/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:954788/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:975754/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:998938/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:152436/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1036188/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1049509/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1130064/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1150626/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1160708/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1305882/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1330428/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1375904/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1384950/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:139676/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1403010/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1497223/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1527527/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1536550/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1572198/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1578988/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:193416/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1645512/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1608029/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1601392/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1575723/61‑1 (MQ=255) cAAGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGAAGCAAATTTTTccg     <  1:1568435/61‑1 (MQ=255)  aaGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:98724/60‑1 (MQ=255)  aaGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:676545/60‑1 (MQ=255)  aaGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:1192371/60‑1 (MQ=255)  aaGTTCCAAAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccg     <  1:84059/60‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:842914/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:303772/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:1242028/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:1247763/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:585918/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:1606796/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:834588/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:379582/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:586345/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:41938/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:509789/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:1593733/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:1589000/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:750945/62‑1 (MQ=255)  aaGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTCTccgccg  <  1:973436/62‑1 (MQ=255)   aGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:262301/61‑1 (MQ=255)   aGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:1047683/61‑1 (MQ=255)   aGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTccgccg  <  1:197584/61‑1 (MQ=255)         |                                                        CAAGTTCC‑AAAGAAGCCAGCGCACCATAAGTATTTTTATTTTCAGCAGCAAATTTTTCCGCCG  >  W3110S.gb/2637657‑2637719

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