Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,512,691	(G)5→4	coding (650/1248 nt)	ypdG ←	predicted enzyme IIC component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,512,687	0	G	.	100.0%	92.1 / NA	23	coding (654/1248 nt)	ypdG	predicted enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (17/6);  total (17/6)

CATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAA  >  W3110S.gb/2512618‑2512747                                                                      |                                                             cATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg                                                             <  1:3312276/71‑1 (MQ=255) cATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg                                                             <  1:2000457/71‑1 (MQ=255) cATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg                                                             <  1:3489305/71‑1 (MQ=255) cATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAgg                                                             <  1:3145460/71‑1 (MQ=255)     aGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAggggt                                                          <  1:3791124/70‑2 (MQ=255)     aGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAggggt                                                          <  1:2771850/70‑2 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:3647977/1‑60 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:3336900/1‑60 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:3228499/1‑60 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:3582250/1‑60 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:2908094/1‑60 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:3764454/1‑60 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:1849578/1‑60 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:1679057/1‑60 (MQ=255)                           cAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:1451192/1‑60 (MQ=255)                             gTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAggggt                                                          >  1:3474599/1‑45 (MQ=255)                                aGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:3329869/1‑55 (MQ=255)                                aGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACgtagta                                             >  1:2249357/1‑55 (MQ=255)                                        aCGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATgg                      >  1:573287/1‑70 (MQ=255)                                        aCGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGc                     >  1:1245588/1‑71 (MQ=255)                                        aCGGATACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGc                     >  1:3810030/1‑71 (MQ=255)                                             tACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACaaaaa                         >  1:2816571/1‑62 (MQ=255)                                             tACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACaaaaa                         >  1:1561501/1‑62 (MQ=255)                                             tACCGCCGTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACaaaaa                         >  1:414401/1‑62 (MQ=255)                                                cgccgTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGaa             <  1:2999739/71‑1 (MQ=255)                                                cgccgTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGaa             <  1:439369/71‑1 (MQ=255)                                                  ccgTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATc           <  1:3431134/71‑1 (MQ=255)                                                  ccgTTGATCCAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATc           <  1:3165954/71‑1 (MQ=255)                                                           cAGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCaaa  <  1:3545169/71‑1 (MQ=255)                                                            aGCCACCAAA‑GGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGaa             <  1:795471/59‑1 (MQ=255)                                                                      |                                                             CATCAGCGCGCCCTTCTCACCTGCGGCAGTCAGCACAGTACGGATACCGCCGTTGATCCAGCCACCAAAGGGGGTGATGACGTAGTACATCGCCAGCATGACAAAAATGGCGGAAAGAATCGGGAGCAAA  >  W3110S.gb/2512618‑2512747

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