Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,061,031	A→G	78.3%	Q509R (CAG→CGG)	yliK →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,061,031	0	A	G	78.3%	45.7 / 7.5	23	Q509R (CAG→CGG)	yliK	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (12/6);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACC  >  W3110S.gb/3060968‑3061086                                                                |                                                        ttCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCAt                                                     <  1:465815/68‑1 (MQ=255) ttCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCAt                                                     <  1:450861/68‑1 (MQ=255) ttCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCAt                                                     <  1:1509438/68‑1 (MQ=255)   cgcgCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGATCGCCCTGACTCACGCCGCACGGCATAACg                                                 >  1:1898948/1‑70 (MQ=255)   cgcgCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACg                                                 >  1:1755776/1‑70 (MQ=255)   cgcgCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACg                                                 >  1:781904/1‑70 (MQ=255)   cgcgCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACg                                                 >  1:165685/1‑70 (MQ=255)   cgcgCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACg                                                 >  1:1453696/1‑70 (MQ=255)   cgcgCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACg                                                 >  1:1372716/1‑70 (MQ=255)   cgcgCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACg                                                 >  1:1281493/1‑70 (MQ=255)     cgcCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGaaa                                              <  1:1427383/71‑1 (MQ=255)     cgcCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGaaa                                              <  1:1190720/71‑1 (MQ=255)                atgatgCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTgctggct                                   >  1:764335/1‑71 (MQ=255)                     gccgccGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGc                              >  1:1103998/1‑71 (MQ=255)                     gccgccGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGc                              >  1:196027/1‑71 (MQ=255)                     gccgccGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGc                              >  1:416091/1‑71 (MQ=255)                        gccgTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCAcatcat                                                     <  1:667045/45‑1 (MQ=255)                         ccgTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGtt                          <  1:1363321/71‑1 (MQ=255)                         ccgTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGtt                          <  1:1119709/71‑1 (MQ=255)                          cgTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGtt                          <  1:232485/70‑1 (MQ=255)                          cgTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGtt                          <  1:362637/70‑1 (MQ=255)                              aCCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGc                     >  1:1380588/1‑71 (MQ=255)                                         gAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTc          >  1:586070/1‑71 (MQ=255)                                                  gACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCAcc  >  1:1536824/1‑70 (MQ=255)                                                                |                                                        TTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACC  >  W3110S.gb/3060968‑3061086

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