Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,493,301	+A	coding (851/1080 nt)	frwC ←	predicted enzyme IIC component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,493,298	1	.	A	100.0%	126.4 / NA	36	C285F (TGT→TTT)	frwC	predicted enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (0/36);  total (0/36)

CCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑AAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGG  >  W3110S.gb/3493247‑3493356                                                     |                                                           cgccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:408222/63‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1857007/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:343552/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:34045/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:308111/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:276201/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:268998/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:2393255/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:488411/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:2360807/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:588275/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:2180513/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:2118387/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:2015619/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1900167/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:67335/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:833372/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:174814/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:952991/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1616369/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1588017/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1430771/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1427741/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1415748/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1388980/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:999808/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1213623/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1191324/65‑1 (MQ=255) ccccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1090506/65‑1 (MQ=255)       ccTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:2396660/59‑1 (MQ=39)        cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:1098391/58‑1 (MQ=39)         tGCATAGCTCTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                <  1:133208/57‑1 (MQ=38)          gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:2370549/65‑1 (MQ=255)          gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:2209895/65‑1 (MQ=255)          gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:1844757/65‑1 (MQ=255)           cATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTTAACCGa                                                <  1:367541/55‑1 (MQ=38)                                               caccacaAAACTGAACCGACCTTGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg  >  1:669774/1‑65 (MQ=255)                                               caccacaAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg  >  1:83198/1‑65 (MQ=255)                                               caccacaAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg  >  1:1684004/1‑65 (MQ=255)                                               caccacaAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAg   >  1:184403/1‑64 (MQ=255)                                                     |                                                           CCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑AAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGG  >  W3110S.gb/3493247‑3493356

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