BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	327,202:1	+A	100%	coding (679/963 nt)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	327,202	1	.	A	100.0%	63.7 / NA	18	E227K (GAA→AAA)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base A (0/18); total (0/18)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGCTTCCGCATTGGGGATGGCACCGG‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGAC  >  minE/327176‑327267

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTTGCGCTTTGGTCggg                            <  1:1249947/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:920150/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:213483/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:195150/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:1938325/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:19098/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:506132/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:1746191/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:1743069/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:624869/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:721503/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:1515732/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:1337759/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:1325300/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:1284415/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:745190/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:1111394/67‑1 (MQ=255)

tGGCTTCCGCATTGGGGATGGCACCGGAAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCggg                            <  1:762440/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTTGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1054111/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:729/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:6592/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:516055/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:353994/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1901589/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:2059303/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:2040360/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:2019084/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1920292/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1006225/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:167461/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1632158/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1602465/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1525250/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1419046/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1192616/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1066414/67‑1 (MQ=255)

                         gg‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACcgac  <  1:1006338/67‑1 (MQ=255)

TGGCTTCCGCATTGGGGATGGCACCGG‑AAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGAC  >  minE/327176‑327267

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