Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,061,031	A→G	73.5%	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	73.5%	57.6 / 20.2	34	Q509R (CAG→CGG)	yliK	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (9/16);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

ATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTGGGTGA  >  W3110S.gb/3060967‑3061094                                                                 |                                                                aTTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCAt                                                             <  1:183138/69‑1 (MQ=255) aTTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCAt                                                             <  1:196494/69‑1 (MQ=255)    cgcgCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATaa                                                           >  1:240967/1‑68 (MQ=255)      cgcCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACg                                                         <  1:728944/68‑1 (MQ=255)      cgcCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACg                                                         <  1:694089/68‑1 (MQ=255)         cACCCGTGATGATGCCGCCTTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGaaaa                                                     <  1:743813/69‑1 (MQ=255)         cACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGaaaa                                                     <  1:618147/69‑1 (MQ=255)         cACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGaaaa                                                     <  1:2424551/69‑1 (MQ=255)            ccGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACct                                                  >  1:2963390/1‑69 (MQ=255)             cGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACctg                                                 >  1:560509/1‑69 (MQ=255)             cGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACct                                                  <  1:2887414/68‑1 (MQ=255)                   gatgCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTgctggct                                           <  1:1133674/69‑1 (MQ=255)                           cgTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGt                                   <  1:2463739/69‑1 (MQ=255)                            gTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGtt                                  <  1:743217/69‑1 (MQ=255)                               aCCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAAt                               >  1:3172804/1‑69 (MQ=255)                               aCCGCCGCGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAAt                               >  1:18288/1‑69 (MQ=255)                                   ccgcGTTGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCg                           <  1:2053718/69‑1 (MQ=255)                                        ttGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGtt                                  <  1:1413734/57‑1 (MQ=255)                                        ttGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGtt                                  <  1:2895392/57‑1 (MQ=255)                                        ttGAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGtt                                  <  1:560685/57‑1 (MQ=255)                                          gAACGCCCTGACTCACGCCGCACGGCATAACGAAAACctgctg                                              >  1:2286615/1‑43 (MQ=255)                                          gAACGCCCTGACTCACGCCGCACGGCATAACGAAAACctgctg                                              >  1:2190700/1‑43 (MQ=255)                                          gAACGCCCTGACTCACGCCGCACGGCATAACGAAAACctgctg                                              >  1:1491163/1‑43 (MQ=255)                                          gAACGCCCTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGt                    >  1:1920605/1‑69 (MQ=255)                                              gCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGtt                                  <  1:503510/51‑1 (MQ=255)                                              gCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTcgc                >  1:2196500/1‑69 (MQ=255)                                              gCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTcgc                >  1:3498270/1‑69 (MQ=255)                                              gCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTcgc                >  1:3292584/1‑69 (MQ=255)                                                ccTGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCgcgc              <  1:2684347/69‑1 (MQ=255)                                                  tGACTCACGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCa            <  1:2442333/69‑1 (MQ=255)                                                      tCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCAcc          <  1:545786/67‑1 (MQ=255)                                                      tCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCAcc          <  1:2110590/67‑1 (MQ=255)                                                        aCGCCGCACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTgg      >  1:42080/1‑69 (MQ=255)                                                            cgcACGGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTGGGTGa  <  1:3342065/69‑1 (MQ=255)                                                                 |                                                                ATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAATGCCGCTCGCGTTCGCGCCACCCTGGGTGA  >  W3110S.gb/3060967‑3061094

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