BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	7,496	C→G	32.7%	G155A (GGC→GCC)	yaaJ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	7,496	0	C	G	32.7%	54.0 / 35.9	52	G155A (GGC→GCC)	yaaJ	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (26/9); new base G (10/7); total (36/16)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCA  >  minE/7428‑7550

ggggAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGc                                                        <  1:2000316/69‑1 (MQ=255)

                 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)

GGGGAAAATCAAAAGAAAAACTCAGGGCGCGGGCAACGGCGTTCGCTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCA  >  minE/7428‑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: `‑`