Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,588,010	T→G	D64A (GAC→GCC)	ydeP ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,588,010	0	T	G	90.6%	100.8 / ‑0.2	32	D64A (GAC→GCC)	ydeP	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (14/15);  total (15/17)
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.51e-01

AGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATT  >  W3110S.gb/1587942‑1588077                                                                     |                                                                    ataaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:458133/69‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:1625628/71‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:287801/71‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:2824747/71‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:3064058/71‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:2492505/71‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:3807634/71‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:1974826/71‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:1303302/71‑1 (MQ=255) agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:101382/71‑1 (MQ=255)       agGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCg                                                                   <  1:1581023/65‑1 (MQ=255)              tACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACtgt                                                      >  1:1488452/1‑71 (MQ=255)                      tATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGa                                              >  1:1379143/1‑71 (MQ=255)                      tATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAgg                                               >  1:1466309/1‑70 (MQ=255)                          cGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCTAATGACGCACTGTGCTTAGGa                                              <  1:3356527/67‑1 (MQ=255)                          cGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGa                                              <  1:1508050/67‑1 (MQ=255)                          cGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTg                                          <  1:3958401/71‑1 (MQ=255)                          cGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTg                                          <  1:569821/71‑1 (MQ=255)                              aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCa                                      >  1:2666288/1‑71 (MQ=255)                                        cGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAAccc                            >  1:688209/1‑71 (MQ=255)                                        cGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAAccc                            >  1:3352190/1‑71 (MQ=255)                                         gATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCg                           >  1:251992/1‑71 (MQ=255)                                         gATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCg                           >  1:2213605/1‑71 (MQ=255)                                         gATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCg                           >  1:2147841/1‑71 (MQ=255)                                         gATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCg                           >  1:3935939/1‑71 (MQ=255)                                         gATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCg                           >  1:3964337/1‑71 (MQ=255)                                         gATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCg                           >  1:1575162/1‑71 (MQ=255)                                                           ttCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg         >  1:1897489/1‑71 (MQ=255)                                                           ttCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg         >  1:1380309/1‑71 (MQ=255)                                                           ttCACAAATGGCCAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg         >  1:2106514/1‑71 (MQ=255)                                                                  aTGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTAtt  <  1:3148211/71‑1 (MQ=255)                                                                  aTGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTAtt  <  1:1262589/71‑1 (MQ=255)                                                                     |                                                                    AGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATT  >  W3110S.gb/1587942‑1588077

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