BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	163,666	C→T	72.4%	E77K (GAA→AAA)	yaeI ←	predicted phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	163,666	0	C	T	72.4%	198.2 / 91.1	116	E77K (GAA→AAA)	yaeI	predicted phosphatase
Reads supporting (aligned to +/- strand): ref base C (22/10); new base T (29/55); total (51/65)
Fisher's exact test for biased strand distribution p-value = 1.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

CGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGT  >  minE/163610‑163719

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:282096/67‑1 (MQ=255)

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:1493745/67‑1 (MQ=255)

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:1587692/67‑1 (MQ=255)

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:1193551/67‑1 (MQ=255)

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:1195094/67‑1 (MQ=255)

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:842875/67‑1 (MQ=255)

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:783797/67‑1 (MQ=255)

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:1372713/67‑1 (MQ=255)

cGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAaga                                             <  1:1293531/67‑1 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                      >  1:161761/1‑48 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:720754/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:526984/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:250082/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:59258/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:1597742/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:1610651/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:2128314/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:1640909/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:846688/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:44075/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:1142081/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:1834930/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:1665148/1‑59 (MQ=255)

          aataatACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGc                                           >  1:292139/1‑59 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:457639/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:520230/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:1570479/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:1566901/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:1978996/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:2150027/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:766008/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:33987/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:937618/1‑66 (MQ=255)

                       cgccgcCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                       >  1:1076466/1‑66 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:273760/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:392025/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:1766787/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:1808804/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:535537/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:164095/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:155559/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:192473/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:679070/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:810493/64‑1 (MQ=255)

                                      tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaaa          <  1:1086557/64‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:2099524/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:2109479/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1051201/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:316824/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:412405/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:420405/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:533873/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:567439/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:681328/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:690331/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:695886/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:947400/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:96109/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1805285/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1484376/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1506378/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1536915/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1562729/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1660667/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1481950/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1481868/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1709808/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1481099/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1398762/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:133271/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1279724/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1275472/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1198483/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1191644/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1148242/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:2000076/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:2002816/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:2007692/60‑1 (MQ=255)

                                                cTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:2091563/60‑1 (MQ=255)

                                                cTTATGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:504932/60‑1 (MQ=255)

                                                 ttttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:790098/59‑1 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAgg                >  1:2060256/1‑46 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAgg                >  1:2450/1‑46 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACgaga    <  1:1086088/58‑1 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1222360/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1027173/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:951039/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1949343/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:260042/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:121382/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1905246/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1694635/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:186020/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:747185/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1826506/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:705524/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1818756/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1729917/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1729315/1‑60 (MQ=255)

                                                  tttGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGt  >  1:1704407/1‑60 (MQ=255)

CGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGT  >  minE/163610‑163719

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