BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,795,792	A→T	40.3%	A94A (GCA→GCT)	dgkA →	diacylglycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,795,792	0	A	T	40.3%	65.2 / 116.8	119	A94A (GCA→GCT)	dgkA	diacylglycerol kinase
Reads supporting (aligned to +/- strand): ref base A (30/41); new base T (29/19); total (59/60)
Fisher's exact test for biased strand distribution p-value = 6.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

Only 100 of 119 total aligned reads displayed.

TAGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTAC  >  minE/2795729‑2795845

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2957704/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2495264/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:3069381/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:586126/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:1512248/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2117609/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2091501/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:1810603/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:793641/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:1268530/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:1209668/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2807841/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:3185270/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2017959/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACACGCTAAAg                                                   <  1:548848/68‑1 (MQ=255)

tAGCGCCATCGAAGCTGTGGGTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCTAAAg                                                   <  1:2166783/68‑1 (MQ=255)

tAGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGa                                                  <  1:3224604/69‑1 (MQ=255)

 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGGGCAAAAGat                                                 >  1:1743336/1‑69 (MQ=255)

 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGat                                                 >  1:642285/1‑69 (MQ=255)

 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGat                                                 >  1:3260999/1‑69 (MQ=255)

 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGat                                                 >  1:3010857/1‑69 (MQ=255)

 aGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGat                                                 >  1:1328533/1‑69 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAg                                                   <  1:1736464/58‑1 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAg                                                   <  1:2375253/58‑1 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAg                                                   <  1:2325062/58‑1 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAg                                                   <  1:1082414/58‑1 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:2057328/1‑69 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:3133073/1‑69 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:3367935/1‑69 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:1630753/1‑69 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:2569417/1‑69 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:927535/1‑69 (MQ=255)

          gAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATcc                                        >  1:1031045/1‑69 (MQ=255)

                  ggTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGt                                >  1:2431518/1‑69 (MQ=255)

                  ggTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGt                                >  1:1387072/1‑69 (MQ=255)

                  ggTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGt                                >  1:2232713/1‑69 (MQ=255)

                   gTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTg                               >  1:414339/1‑69 (MQ=255)

                   gTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTg                               >  1:1911691/1‑69 (MQ=255)

                     tGACCTAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:2032651/69‑1 (MQ=255)

                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:2905930/69‑1 (MQ=255)

                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:2298835/69‑1 (MQ=255)

                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:692907/69‑1 (MQ=255)

                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:675164/69‑1 (MQ=255)

                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:1867340/69‑1 (MQ=255)

                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:428067/69‑1 (MQ=255)

                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:3394617/69‑1 (MQ=255)

                     tGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCt                             <  1:2176530/69‑1 (MQ=255)

                         cGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                         <  1:2596272/69‑1 (MQ=255)

                         cGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                         <  1:2702086/69‑1 (MQ=255)

                         cGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                         <  1:2716873/69‑1 (MQ=255)

                          gAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                         <  1:1669261/68‑1 (MQ=255)

                               ggCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   <  1:3255607/69‑1 (MQ=255)

                               ggCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   <  1:2667484/69‑1 (MQ=255)

                                  tctGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTAtc                <  1:1586865/69‑1 (MQ=255)

                                  tctGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTAtc                <  1:593750/69‑1 (MQ=255)

                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:912675/1‑62 (MQ=255)

                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:1068364/1‑62 (MQ=255)

                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:2909098/1‑62 (MQ=255)

                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:3045541/1‑62 (MQ=255)

                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:3115381/1‑62 (MQ=255)

                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:2726505/1‑62 (MQ=255)

                                      aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                   >  1:2113389/1‑62 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:2594014/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:3228104/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:3232082/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:2430836/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:290007/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:3000537/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:207964/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:128099/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:823060/68‑1 (MQ=255)

                                            aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:347769/68‑1 (MQ=255)

                                             tGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg       <  1:289025/67‑1 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1502563/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1996311/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1925583/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:547318/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1565471/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1521982/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2116988/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1492644/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1367059/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1331839/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:783812/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:1118890/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:3363858/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:3259407/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:231508/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2359348/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2476793/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2548072/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2582087/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2990289/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2780125/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGAtt    >  1:2904703/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCTAAAGATATGGGATCTGCCGCAGTGCTAATTGCCATTATCGTCGCCGTGAtt    >  1:388389/1‑69 (MQ=255)

                                              gAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGAtt    <  1:1155618/69‑1 (MQ=255)

                                               aGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTa   >  1:3146831/1‑69 (MQ=255)

                                               aGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTa   >  1:2480845/1‑69 (MQ=255)

                                                gCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCCGTGCTGATTGCCATTATCGTCGCCGTGATTAc  >  1:2957406/1‑69 (MQ=255)

TAGCGCCATCGAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTAC  >  minE/2795729‑2795845

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