BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,724,675	C→T	73.4%	K165K (AAG→AAA)	yheT ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,724,675	0	C	T	73.4%	45.8 / 15.2	45	K165K (AAG→AAA)	yheT	predicted hydrolase
Reads supporting (aligned to +/- strand): ref base C (8/4); new base T (14/19); total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCCAGCATAAACGGCGCAGAGACAATCACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGC  >  minE/2724614‑2724725

cTTCCAGCATAAACGGCGCAGAGACAATCACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGcc                                              >  1:2167982/1‑68 (MQ=255)

           aaCGGCGCAGAGACAATCACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGcc                                  >  1:1677180/1‑69 (MQ=255)

           aaCGGCGCAGAGACAATCACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGcc                                  >  1:369309/1‑69 (MQ=255)

             cGGCGCAGAGACAATCACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGCCAg                                >  1:1904846/1‑69 (MQ=255)

                cgcAGAGACAATCACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGCCAGCa                              <  1:225574/68‑1 (MQ=255)

                            cACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCaaa                                        <  1:2967743/46‑1 (MQ=255)

                            cACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCaaa                                        <  1:1941549/46‑1 (MQ=255)

                            cACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCaaa                                        <  1:165656/46‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCGAGCATATTACCGCCGAGCGaa             <  1:2617943/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:3234962/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:3002494/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:3014614/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:3052821/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:1217206/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:3276582/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:459510/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:568381/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:617593/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:666198/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:740756/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:2099490/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:1353626/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:1574816/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:1744418/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:1814779/69‑1 (MQ=255)

                                gccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:1925309/69‑1 (MQ=255)

                                 ccgcATCAATCGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:2473313/68‑1 (MQ=255)

                                      tCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCg       >  1:1904993/1‑69 (MQ=255)

                                      tCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCg       >  1:2191007/1‑69 (MQ=255)

                                      tCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCg       >  1:3200097/1‑69 (MQ=255)

                                      tCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCg       >  1:1720284/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:777163/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:1032910/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:705840/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:628881/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:1468558/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:1626001/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:3372575/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:3346832/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:3299184/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:2093844/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:3195315/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:1856849/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCCAGCATATTACCGCAGAGCGAATAGCCGACGGc  >  1:3094047/1‑69 (MQ=255)

                                           cGGGAGATCATTGCCTTCTTTTGCCAGCAAACAGGCAAGCATATTACCGCCGAGCGAATAGCCGACGGc  >  1:1957618/1‑69 (MQ=255)

                                                        ccttctttTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGaa             <  1:3021406/45‑1 (MQ=255)

CTTCCAGCATAAACGGCGCAGAGACAATCACCGCCGCATCAACCGGGAGATCATTGCCTTCCTTTGCCAGCAAACAGGCCAGCATATTACCGCCGAGCGAATAGCCGACGGC  >  minE/2724614‑2724725

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