Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	500,029	G→A	43.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	43.2%	19.3 / 52.7	44	L227L (CTG→CTA)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base G (12/13);  new base A (11/8);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAA  >  W3110S.gb/499961‑500083                                                                     |                                                       aTGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTagg                                                      >  1:1847414/1‑71 (MQ=255)    ccGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCAc                                                   >  1:4205769/1‑71 (MQ=255)        aaGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAg                                               >  1:3421249/1‑71 (MQ=255)         aGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGt                                              >  1:1915005/1‑71 (MQ=255)         aGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGt                                              >  1:662302/1‑71 (MQ=255)          gCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGtt                                             >  1:3306631/1‑71 (MQ=255)           cAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGttt                                            <  1:1543379/71‑1 (MQ=255)             aCCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGt                                          <  1:1294041/71‑1 (MQ=255)              ccATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGt                                          >  1:1770859/1‑70 (MQ=255)               cATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCa                                        >  1:2625983/1‑71 (MQ=255)                   aaaGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGc                                    >  1:2544480/1‑71 (MQ=255)                      gCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGt                                          >  1:1158178/1‑62 (MQ=255)                       ccATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCgag                                <  1:2467207/71‑1 (MQ=255)                        cATTGAGTACGCGAAGAAATATAACGTACCGATGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCgag                                <  1:3130225/70‑1 (MQ=255)                         aTTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGaa                              <  1:1655689/71‑1 (MQ=255)                         aTTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTc                                         <  1:3757622/60‑1 (MQ=255)                          ttGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGt                                          <  1:3940329/58‑1 (MQ=255)                          ttGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTga                                    >  1:1868351/1‑63 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGtt                                             >  1:4307474/1‑53 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGtt                                             >  1:880033/1‑53 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATcc                           >  1:3457647/1‑71 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATcc                           >  1:770529/1‑71 (MQ=255)                            gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATc                            >  1:2937227/1‑70 (MQ=255)                             aGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCg                          <  1:3035446/71‑1 (MQ=255)                             aGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATCCg                          <  1:3920553/71‑1 (MQ=255)                                 cgcgAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGt                                          >  1:385293/1‑51 (MQ=255)                                 cgcgAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAgt                      >  1:1983391/1‑71 (MQ=255)                                 cgcgAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGt                                          <  1:1301764/51‑1 (MQ=255)                                 cgcgAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATCCGCAgt                      <  1:1673000/71‑1 (MQ=255)                                 cgcgAAGAAATATAACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATCCGCAgt                      >  1:399335/1‑71 (MQ=255)                                     aagaaATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAg                                               <  1:761518/42‑1 (MQ=255)                                       gaaATATTACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGt                                          <  1:4647313/45‑1 (MQ=255)                                        aaaTATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGt                                          <  1:1870439/44‑1 (MQ=255)                                        aaaTATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGtggtgg                 >  1:4064832/1‑69 (MQ=255)                                          atatAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGgcag              >  1:2609259/1‑70 (MQ=255)                                           tataACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCgaga                               <  1:533157/52‑1 (MQ=255)                                           tataACGTACCGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGgcagca            <  1:2557692/71‑1 (MQ=255)                                             taACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCAtt                                      <  1:3980464/43‑1 (MQ=255)                                             taACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAAt          >  1:1292722/1‑71 (MQ=255)                                               aCGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAAtt         <  1:1466358/70‑1 (MQ=255)                                                    ccGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCa    <  1:3720859/70‑1 (MQ=255)                                                    ccGGTGGTGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCa    <  1:825396/70‑1 (MQ=255)                                                      ggtggtGCTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATc                            >  1:1591045/1‑44 (MQ=255)                                                             cTGACGCTAGGCACCAAGTTTGTCATTGCCGAGAATCCGCAg                       <  1:290411/42‑1 (MQ=255)                                                                aCGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCaaa  >  1:4306214/1‑60 (MQ=255)                                                                     |                                                       ATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAA  >  W3110S.gb/499961‑500083

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