BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,316,170	G→A	32.1%	A496A (GCG→GCA)	glpK →	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,316,170	0	G	A	32.1%	107.9 / 61.7	121	A496A (GCG→GCA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (43/39); new base A (29/10); total (72/49)
Fisher's exact test for biased strand distribution p-value = 2.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GAAACCACTGAGCGTAATTACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCC  >  minE/2316108‑2316234

gAAACCACTGAGCGTAATTACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTggg                                                              >  1:659560/1‑67 (MQ=255)

gAAACCACTGAGCGTAATTACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGga                                                             >  1:1232559/1‑68 (MQ=255)

gAAACCACTGAGCGTAATTACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGga                                                             >  1:1484684/1‑68 (MQ=255)

gAAACCACTGAGCGTAATTACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGga                                                             >  1:2938791/1‑68 (MQ=255)

gAAACCACTGAGCGTAATTACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGga                                                             >  1:1130684/1‑68 (MQ=255)

gAAACCACTGAGCGTAATTACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGga                                                             >  1:1017303/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                     aaaaaGCGGTTAAACGCGCGATGGCGTGGGAAGACCACGACGAATAATGTAAATGCCGAATGAAGCGtt                       >  1:403282/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                              ttAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTa                                         <  1:3200481/42‑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:2034669‑M2/1‑35 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

GAAACCACTGAGCGTAATTACCGTTACGCAGGCTGGAAAAAAGCGGTTAAACGCGCGATGGCGTGGGAAGAACACGACGAATAATGTAAATGCCGAATGAAGCGTTTATGCCGCATCCGGTAGTCCC  >  minE/2316108‑2316234

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