Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,952,349	G→C	R369P (CGG→CCG)	bcsC →	cellulose synthase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,952,349	0	G	C	100.0%	123.8 / NA	67	R369P (CGG→CCG)	bcsC	cellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (0/67);  total (0/67)

GGCTGGCGATCCAGCAGGGCGATGCTGCGCTGAAAGCCAATAATCCTGACCGGGCAGAAC  >  W3110S.gb/3952298‑3952357                                                    |         ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATACTCCTGACCCGGCAGAAc  <  1:50374/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGTCAGAAc  <  1:23413/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:64249/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:74549/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:73496/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:73401/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:70783/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:69330/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:66387/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:64967/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:77596/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:61094/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:60676/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:60322/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:56781/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:10938/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:54926/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:55989/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:78896/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:79535/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:80100/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:8028/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:81835/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:81942/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:82149/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:82666/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:83154/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:84095/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:86954/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:86961/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:87661/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:9585/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:36722/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:12512/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:12850/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:14447/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:16497/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:18055/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:24498/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:2499/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:26322/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:26479/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:26748/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:27304/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:27670/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:32236/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:39264/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:50133/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:54267/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:54223/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:53875/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:53363/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:51751/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:50382/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:40705/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:49300/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:46659/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:45/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:43395/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCACAAc  <  1:35925/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGAACCGGCAGAAc  <  1:73051/60‑1 (MQ=255) ggctggcGATCCAGCAGGGCGAAGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:67403/60‑1 (MQ=38)  gctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:80117/59‑1 (MQ=39)  gctggcGATCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:58238/59‑1 (MQ=39)         aTCCAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:86186/52‑1 (MQ=39)           ccAGCAGGGCGATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:87467/50‑1 (MQ=38)                     gATGCCGCGCTGAAAGCCAATAATCCTGACCCGGCAGAAc  <  1:547/40‑1 (MQ=38)                                                    |         GGCTGGCGATCCAGCAGGGCGATGCTGCGCTGAAAGCCAATAATCCTGACCGGGCAGAAC  >  W3110S.gb/3952298‑3952357

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