Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,248,498	G→A	R152R (CGC→CGT)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,248,498	0	G	A	100.0%	263.5 / NA	95	R152R (CGC→CGT)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (95/0);  total (95/0)

TACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCC  >  W3110S.gb/1248486‑1248547             |                                                  taCCTCTCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1234696/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATagag        >  1:990175/1‑56 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATag          >  1:1217029/1‑54 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATTGAGAATccc  >  1:259568/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGTGAATccc  >  1:1289247/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:391877/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:545538/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:538836/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:475985/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:474502/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:452077/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:443481/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:440133/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:43871/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:409814/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:396086/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:577554/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:385002/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:378632/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:378064/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:375883/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:374401/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:338604/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:337396/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:330409/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:322713/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:305711/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:736190/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:991782/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:982275/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:977969/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:971375/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:839822/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:833609/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:814276/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:782600/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:778211/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:776453/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:740278/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:548670/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:704490/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:689046/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:65188/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:649844/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:646567/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:638029/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:621154/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:582055/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:275800/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:575279/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:111030/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:129571/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1272119/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1235873/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1230691/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1227741/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:121575/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:120423/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1190487/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1182997/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1151874/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:11253/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:13134/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1107853/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1093849/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1090710/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1087142/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1085026/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1059947/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1054577/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1043681/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1029284/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1027309/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:151349/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:292169/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1000665/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:275413/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:26276/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:249321/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:243423/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:223945/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:22223/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:207723/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:186656/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:302567/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1467476/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1449345/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1446625/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1442315/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1410625/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1408525/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1380970/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1380230/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:1365526/1‑62 (MQ=255) taCCTCGCGAAAACGTCCGCCGGGCACGACATAAGGTTTTGGCAGCGGTAATAGAGAATccc  >  1:135272/1‑62 (MQ=255)             |                                                  TACCTCGCGAAAGCGTCCGCCCGGCACGACATAAGGTTCCGGCAGCGGTAACAGAGAATCCC  >  W3110S.gb/1248486‑1248547

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