Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	500,029	G→A	47.2%	L227L (CTG→CTA)	gsk →	inosine/guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	500,029	0	G	A	47.2%	0.9 / 48.1	36	L227L (CTG→CTA)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base G (12/7);  new base A (12/5);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

ATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCAC  >  W3110S.gb/499961‑500089                                                                     |                                                             aTGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTggg                                                            <  1:2200023/71‑1 (MQ=255)    ccGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCAc                                                         >  1:1281424/1‑71 (MQ=255)          gCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGtt                                                   >  1:3959723/1‑71 (MQ=255)          gCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGtt                                                   >  1:3412070/1‑71 (MQ=255)             aCCATGAAACCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGt                                                <  1:3492552/71‑1 (MQ=255)              ccATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGttt                                                  <  1:3595556/68‑1 (MQ=255)              ccATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGttt                                                  <  1:2674797/68‑1 (MQ=255)               cATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCa                                              >  1:489183/1‑71 (MQ=255)                aTGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCAt                                             >  1:1379726/1‑71 (MQ=255)                aTGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCAt                                             >  1:1218924/1‑71 (MQ=255)                   aaaGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTg                                           <  1:3925322/70‑1 (MQ=255)                         aTTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCgaga                                     >  1:798947/1‑70 (MQ=255)                         aTTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGaa                                    >  1:113383/1‑71 (MQ=255)                         aTTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGaa                                    >  1:3390188/1‑71 (MQ=255)                          ttGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAAt                                   >  1:2900907/1‑71 (MQ=255)                          ttGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAAt                                   >  1:1843466/1‑71 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATcc                                 >  1:1062699/1‑71 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATcc                                 >  1:2187465/1‑71 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATcc                                 >  1:2796558/1‑71 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATcc                                 >  1:2954228/1‑71 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATcc                                 >  1:642777/1‑71 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATc                                  >  1:736445/1‑70 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATc                                  >  1:3725080/1‑70 (MQ=255)                             aGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCg                                <  1:3481875/71‑1 (MQ=255)                             aGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCg                                <  1:1906028/71‑1 (MQ=255)                                 cgcgAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCAt                                             >  1:3414724/1‑54 (MQ=255)                                 cgcgAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAgt                            <  1:3473704/71‑1 (MQ=255)                                 cgcgAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATCCGCAgt                            <  1:1765312/71‑1 (MQ=255)                                  gcgAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAgtg                           <  1:1618811/71‑1 (MQ=255)                                       gaaATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCAtt                                            <  1:1652917/49‑1 (MQ=255)                                            ataACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCaa                 >  1:2341756/1‑71 (MQ=255)                                            ataACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCaa                 >  1:3947838/1‑71 (MQ=255)                                                    ccGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGgcagc                   <  1:717615/61‑1 (MQ=255)                                                    ccGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGgcagc                   <  1:2480142/61‑1 (MQ=255)                                                         ggtgCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCATTGGTGGCAGCAATTCCTCAAAGAt     >  1:2889143/1‑70 (MQ=255)                                                             cTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCAc  >  1:1676147/1‑69 (MQ=255)                                                             cTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCAc  >  1:1940000/1‑69 (MQ=255)                                                                     |                                                             ATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCAC  >  W3110S.gb/499961‑500089

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