Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,875,612	G→A	D126N (GAT→AAT)	iap →	aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,875,612	0	G	A	100.0%	59.5 / NA	20	D126N (GAT→AAT)	iap	aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/8);  total (12/8)

CCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGA  >  W3110S.gb/2875561‑2875661                                                    |                                                  ccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGa                                      <  1:640057/65‑1 (MQ=255) ccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGa                                      <  1:1727448/65‑1 (MQ=255) ccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGa                                      <  1:2501683/65‑1 (MQ=255) ccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGa                                      <  1:1979590/65‑1 (MQ=255) ccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGa                                      <  1:2445856/65‑1 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:2380740/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:992170/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:885246/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:2832368/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:2598297/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:1171579/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:2113472/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:185170/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:1663844/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:1635779/1‑65 (MQ=255)                                 tACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa      >  1:1266985/1‑65 (MQ=255)                                 tACTTACGACCCGCTGAGCAATGCTGACGCCGATGCCAATCTcg                           >  1:1859795/1‑44 (MQ=38)                                     tACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga  <  1:1916913/65‑1 (MQ=255)                                     tACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga  <  1:2522175/65‑1 (MQ=255)                                     tACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGga  <  1:3063056/65‑1 (MQ=255)                                                    |                                                  CCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGA  >  W3110S.gb/2875561‑2875661

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