BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	7,498	A→G	31.5%	Y154Y (TAT→TAC)	yaaJ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	7,498	0	A	G	31.5%	40.2 / 23.1	52	Y154Y (TAT→TAC)	yaaJ	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (26/9); new base G (9/7); total (36/16)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCA  >  minE/7445‑7550

aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:2144471/69‑1 (MQ=255)

aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:2231443/69‑1 (MQ=255)

aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:3074767/69‑1 (MQ=255)

aaaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:679474/69‑1 (MQ=255)

 aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:901386/68‑1 (MQ=255)

 aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:2216594/68‑1 (MQ=255)

 aaCTCAGGGCTCGCGCAACGGCGTTCGCCTGAACTCCGCTGAAAATTATGGCGTAGGCGATGAGCaaa                                       <  1:1511347/68‑1 (MQ=255)

   cTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAg                                    <  1:2900651/69‑1 (MQ=255)

     cAGGGCGCTGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc                                  <  1:740452/69‑1 (MQ=255)

     cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc                                  <  1:44685/69‑1 (MQ=255)

     cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc                                  <  1:164911/69‑1 (MQ=255)

     cAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGAc                                  <  1:2372362/69‑1 (MQ=255)

              ggCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAg                         <  1:1409568/69‑1 (MQ=255)

              ggCAACGGCGTTCGCTTGAACTCCCCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAg                         <  1:2131063/69‑1 (MQ=255)

                   cGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGc                    <  1:3112990/69‑1 (MQ=255)

                    ggCGTTCTCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:3166875/1‑69 (MQ=255)

                    ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:305304/1‑69 (MQ=255)

                    ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:845844/1‑69 (MQ=255)

                    ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:1041534/1‑69 (MQ=255)

                    ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:1018455/1‑69 (MQ=255)

                    ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:1287037/1‑69 (MQ=255)

                    ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:1520955/1‑69 (MQ=255)

                    ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:172562/1‑69 (MQ=255)

                    ggCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGc                    >  1:2252914/1‑68 (MQ=255)

                          cGCTTGCACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:439998/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:360526/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:2010953/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:293591/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:2994938/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:868008/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:2027817/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:2586750/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:1252826/1‑68 (MQ=255)

                          cGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATcc              >  1:817177/1‑68 (MQ=255)

                               gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:2032330/1‑58 (MQ=255)

                               gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:638351/1‑58 (MQ=255)

                               gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:754247/1‑58 (MQ=255)

                               gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGcc                   >  1:3028450/1‑58 (MQ=255)

                               gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCa        >  1:20065/1‑69 (MQ=255)

                               gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCa        >  1:2753905/1‑69 (MQ=255)

                               gAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCa        >  1:2722561/1‑69 (MQ=255)

                                aaCTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGc                    <  1:2667871/56‑1 (MQ=255)

                                     cGCTGAAAATTATGGCTTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:87435/1‑69 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATTAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:3320544/1‑69 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCa                 >  1:1653504/1‑54 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:622299/1‑69 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:1534329/1‑69 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:239230/1‑69 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:3117022/1‑69 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:302538/1‑69 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCa  >  1:2610482/1‑69 (MQ=255)

                                     cGCTGAAAATTATGGCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGAGCATCCCCa  >  1:1580875/1‑69 (MQ=255)

AAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCA  >  minE/7445‑7550

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