Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,594,183	A→G	61.1%	D28D (GAT→GAC)	yihE ←	predicted kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,594,183	0	A	G	61.1%	10.1 / 15.8	18	D28D (GAT→GAC)	yihE	predicted kinase
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (5/6);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CGTCGACGATCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGGTGT  >  W3110S.gb/3594115‑3594234                                                                     |                                                    cGTCGACGATCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATc                                                   <  1:970603/71‑1 (MQ=255)  gTCGACGATCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATc                                                   <  1:977731/70‑1 (MQ=255)     gacgaTCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAGTCCAcc                                               >  1:2731378/1‑71 (MQ=255)              cGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGC‑‑‑GTAAGACCGGAGTCCACCCGGATcc                                        <  1:1066709/66‑1 (MQ=255)               gTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCa                                      <  1:1453508/70‑1 (MQ=255)                        ttGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACa                            <  1:1877979/71‑1 (MQ=255)                                cacGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAGTCCAc                                                >  1:262550/1‑43 (MQ=255)                                cacGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAGTCCAc                                                >  1:51092/1‑43 (MQ=255)                                  cGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCat                  <  1:2547271/71‑1 (MQ=255)                                  cGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATccc                                       >  1:2196306/1‑50 (MQ=255)                                                 ttAAGCGGGGTAAGACCGGAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtg   <  1:1963169/71‑1 (MQ=255)                                                 ttAAGCGGGGGAAGACCGGAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtg   <  1:846703/71‑1 (MQ=255)                                                   aaGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtg   >  1:2179453/1‑69 (MQ=255)                                                           tAAGACCGGAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCatgat               >  1:764740/1‑49 (MQ=255)                                                             aGACCGGAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGt            >  1:757030/1‑50 (MQ=255)                                                               aCCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCatgatg              >  1:1082527/1‑46 (MQ=255)                                                                ccGGAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtgt  <  1:290713/57‑1 (MQ=255)                                                                  ggAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGt            <  1:1933711/45‑1 (MQ=255)                                                                     |                                                    CGTCGACGATCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGGTGT  >  W3110S.gb/3594115‑3594234

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