BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,316,188	A→G	37.1%	E502E (GAA→GAG)	glpK →	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,316,188	0	A	G	37.1%	69.9 / 78.4	121	E502E (GAA→GAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base A (35/41); new base G (35/10); total (70/51)
Fisher's exact test for biased strand distribution p-value = 6.49e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 122 total aligned reads displayed.

ACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCCGAAACGTGCGGGGGCAACCCCG  >  minE/2316127‑2316256

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1851393/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               >  1:400322/1‑69 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1593552/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1436377/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:3320439/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:3192551/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:3069432/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:565520/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1829109/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1158985/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:852884/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               >  1:1148805/1‑69 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:1855286/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               >  1:2058195/1‑69 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:2666088/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:2119803/69‑1 (MQ=255)

aCCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:231877/69‑1 (MQ=255)

      aCGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGc                                                          >  1:101841/1‑68 (MQ=255)

      aCGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGc                                                          >  1:14539/1‑68 (MQ=255)

        gCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCg                                                        <  1:511232/68‑1 (MQ=255)

        gCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCg                                                        <  1:2248352/68‑1 (MQ=255)

        gCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCg                                                        <  1:913707/68‑1 (MQ=255)

             cTGGAAAAAAGGGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:428772/56‑1 (MQ=255)

             cTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATgca                                                  >  1:2520634/1‑67 (MQ=255)

             cTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGaa                                                  >  1:1527475/1‑69 (MQ=255)

             cTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGaa                                                  >  1:521671/1‑69 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:2509065‑M2/69‑20 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:2661354‑M2/69‑20 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1039767‑M2/69‑20 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:2754557‑M2/69‑20 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1998055‑M2/69‑20 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1111567‑M2/69‑20 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:686114‑M2/69‑20 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1733121‑M2/69‑20 (MQ=255)

                gAAAAAAGCGGTCAAACGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgac                                               <  1:1311068‑M2/69‑20 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:486029/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:3057286/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:342234/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:2308456/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2524643/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:3214827/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2530378/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2608624/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2645272/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:3208947/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2655493/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:3120281/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:853991/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2807395/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:3078246/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2934068/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2988868/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:551255/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:781933/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:2365848/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2336938/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1616157/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2233562/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:1509542/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:2196335/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1645191/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1506235/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:1472134/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:1406415/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             >  1:2015855/1‑69 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1744783/69‑1 (MQ=255)

                  aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGtt                                             <  1:1954386/69‑1 (MQ=255)

                            tAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:73130/41‑1 (MQ=255)

                              aaCGCGCGATTGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:152771/39‑1 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGATGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:3431201‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcttcttat                                >  1:3744‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1197676‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1258210‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1789789‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1061513‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:172837‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:2034669‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:2087966‑M2/1‑35 (MQ=255)

                               aCGCGCGATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:2114535‑M2/1‑35 (MQ=255)

                                cgcgcgATGGCGTGGGAAGAACACGACGAATAATGTa                                                               <  1:2099689/37‑1 (MQ=255)

                                cgcgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1015131‑M2/1‑34 (MQ=255)

                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:919467‑M2/1‑32 (MQ=255)

                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:130344‑M2/1‑32 (MQ=255)

                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1025624‑M2/1‑32 (MQ=255)

                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1696489‑M2/1‑32 (MQ=255)

                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:3156531‑M2/1‑32 (MQ=255)

                                  cgcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1737154‑M2/1‑32 (MQ=255)

                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:3197549‑M2/1‑31 (MQ=255)

                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1476526‑M2/1‑31 (MQ=255)

                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:126346‑M2/1‑31 (MQ=255)

                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:3036409‑M2/1‑31 (MQ=255)

                                   gcgATGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:238090‑M2/1‑31 (MQ=255)

                                       tGGCATGGGAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1152198‑M2/1‑27 (MQ=255)

                                        ggCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTccc                        >  1:2223753/1‑68 (MQ=255)

                                              ggAAGAACACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:3043948‑M2/1‑20 (MQ=255)

                                               gaagaaCACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:1631303‑M2/1‑19 (MQ=255)

                                                   aaCACGACGAGTAATaatgattgcctgatgcgacgcttgcgcgtcttat                                >  1:3311120‑M2/1‑15 (MQ=255)

                                                      acgacgAGTAATaatgattgcctgatgcgacgcttgcgcgtcttatcag                             >  1:1946650‑M2/1‑12 (MQ=255)

                                                      acgacgAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCCGAAACGTGCGGGGGc         <  1:2891201/69‑1 (MQ=255)

                                                             ataatGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCCGAAACGTGCGGGGGAAACCCCg  <  1:384538/69‑1 (MQ=255)

ACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCCGAAACGTGCGGGGGCAACCCCG  >  minE/2316127‑2316256

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