BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,945,840	T→G	S220S (TCT→TCG)	bcsA →	cellulose synthase, catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,945,840	0	T	G	100.0%	35.8 / NA	12	S220S (TCT→TCG)	bcsA	cellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand): ref base T (0/0); new base G (12/0); total (12/0)

GTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCTACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGAAA  >  W3110S.gb/3945795‑3945901

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:1335554/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:630157/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:615789/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:301411/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:2913213/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:2734938/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:2605739/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:2131327/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:2033631/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:193981/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:1696499/1‑62 (MQ=255)

gTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCGACGTTAAACTGGGACg                                               >  1:1614434/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCtt                     >  1:923818/1‑43 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:2761552/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:943484/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:852316/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:727866/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:632774/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:110155/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:625506/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1389352/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:504799/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1417604/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:2908785/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:2827838/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1736389/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1439673/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1604778/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:2295553/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:2238807/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:2042582/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:162803/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:2024425/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1936809/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1789308/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1788914/1‑62 (MQ=255)

                                             tACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGAGGGCTTATTCTGCTCTTCGCTGaaa  >  1:1403197/1‑62 (MQ=255)

GTCGCTGACCGTTTCTTGCCGTTATATCTGGTGGCGTTACACCTCTACGCTGAACTGGGACGATCCGGTCAGCCTGGTGTGCGGGCTTATTCTGCTCTTCGCTGAAA  >  W3110S.gb/3945795‑3945901

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