BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,831,766	T→C	15.4%	A103A (GCA→GCG)	yjeH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,831,766	0	T	C	15.4%	72.2 / 4.7	45	A103A (GCA→GCG)	yjeH	predicted transporter
Reads supporting (aligned to +/- strand): ref base T (14/24); new base C (5/2); total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCA  >  minE/2831699‑2831807

cTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTg                                          >  1:2925885/1‑69 (MQ=255)

cTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTg                                          >  1:346409/1‑69 (MQ=255)

cTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTg                                          >  1:2468097/1‑69 (MQ=255)

cTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTg                                          >  1:1446212/1‑69 (MQ=255)

        aacaGTTGCCCGCTATGCCAGCCAAACATTGCCTGGCCGAATCCGGCAGCAATTTGTAGCGCggcaggc                                  <  1:3156032/69‑1 (MQ=255)

        aacaGTTGCCCGCTATGCCAGCCAAACATTGCCTGGCCGAATCCGGCAGCAATTTGTAGCGCggcaggc                                  <  1:684565/69‑1 (MQ=255)

                   gCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAAcc                             >  1:1718086/1‑63 (MQ=255)

                   gCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACggg                       <  1:3126634/69‑1 (MQ=255)

                   gCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACggg                       <  1:2648071/69‑1 (MQ=255)

                   gCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACggg                       <  1:352997/69‑1 (MQ=255)

                   gCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACggg                       <  1:2006203/69‑1 (MQ=255)

                   gCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACggg                       <  1:1280546/69‑1 (MQ=255)

                   gCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACggg                       <  1:530152/69‑1 (MQ=255)

                   gCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACggg                       <  1:2046868/69‑1 (MQ=255)

                         ccagccaAACATCGCCTTGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGAc                 <  1:1238516/69‑1 (MQ=255)

                         ccagccaAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGa                  >  1:2000002/1‑68 (MQ=255)

                         ccagccaAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGa                  >  1:1047895/1‑68 (MQ=255)

                         ccagccaAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGAc                 <  1:3152102/69‑1 (MQ=255)

                         ccagccaAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGAc                 <  1:1099058/69‑1 (MQ=255)

                          cagccaAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGAc                 <  1:2675693/68‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:1698243/69‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:346552/69‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:974543/69‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:1264692/69‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:2406508/69‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:1256280/69‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:766876/69‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:2706151/69‑1 (MQ=255)

                                aaCATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATaaa          <  1:2757878/69‑1 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagccg  >  1:1761484/1‑68 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcca  >  1:1112262/1‑69 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcca  >  1:3044428/1‑69 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcca  >  1:1096792/1‑69 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcca  >  1:2302704/1‑69 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcc   <  1:518476/68‑1 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcc   >  1:3147359/1‑68 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcc   <  1:2806867/68‑1 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcc   <  1:2605667/68‑1 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAAcagcc   <  1:1168364/68‑1 (MQ=255)

                                        cTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGAGAAAAAcagcca  >  1:2150669/1‑69 (MQ=255)

                                                          aaTTTGTAGCGCGGCAGGCAAACCCACGGGAATGACCGATa            >  1:2344616/1‑41 (MQ=255)

                                                          aaTTTGTAGCGCGGCAGGCAAACCCACGGGAATGACCGATa            >  1:1574796/1‑41 (MQ=255)

                                                          aaTTTGTAGCGCGGCAGGCAAACCCACGGGAATGACCGATa            >  1:1433777/1‑41 (MQ=255)

                                                          aaTTTGTAGCGCGGCAGGCAAACCCACGGGAATGACCGATa            >  1:1219519/1‑41 (MQ=255)

                                                          aaTTTGTAGCGCGGCAGGCAAACCCACGGGAATGACCGATa            >  1:108669/1‑41 (MQ=255)

CTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACAGCCA  >  minE/2831699‑2831807

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