BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,808,875	G→T	L264L (CTG→CTT)	pfkB →	6‑phosphofructokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,808,875	0	G	T	80.4%	90.7 / 26.4	51	L264L (CTG→CTT)	pfkB	6‑phosphofructokinase II
Reads supporting (aligned to +/- strand): ref base G (9/1); new base T (0/41); total (9/42)
Fisher's exact test for biased strand distribution p-value = 3.29e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGC > W3110S.gb/1808850‑1808925

aCAGCATGGTCGGCGCGATGAGACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2611663/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCATAAAATGCCTCTCTTGAAGAGATGGtc                <  1:395895/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2810585/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2360175/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:238315/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1176295/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2426492/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2433354/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2505440/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2588446/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2679542/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:279150/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2807038/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2404344/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:315492/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:409963/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:442762/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:494842/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:529947/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:558933/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:768423/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:899018/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:933712/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1801759/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1230163/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1267917/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1437998/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1575773/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1719856/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1721321/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1747817/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1750578/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2267542/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1814501/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2120176/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2230741/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2181666/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2138897/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:1958023/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:204381/62‑2 (MQ=255)

aCAGCATGGTCGGCGCGATGACACTTAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGtc                <  1:2023821/62‑2 (MQ=255)

   gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGt             >  1:2036101/1‑62 (MQ=255)

   gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGt             >  1:437254/1‑62 (MQ=255)

   gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGt             >  1:317111/1‑62 (MQ=255)

   gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGt             >  1:165042/1‑62 (MQ=255)

   gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGt             >  1:712199/1‑62 (MQ=255)

   gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGt             >  1:818771/1‑62 (MQ=255)

   gCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGt             >  1:2186588/1‑62 (MQ=255)

            gcgcGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTa    >  1:249570/1‑62 (MQ=255)

            gcgcGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTa    >  1:1948554/1‑62 (MQ=255)

              gcgATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGc  <  1:261688/62‑1 (MQ=255)

ACAGCATGGTCGGCGCGATGACACTGAAACTGGCAGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGC > W3110S.gb/1808850‑1808925

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