BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,795,807	C→T	40.4%	S99S (TCC→TCT)	dgkA →	diacylglycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,795,807	0	C	T	40.4%	54.3 / 81.4	116	S99S (TCC→TCT)	dgkA	diacylglycerol kinase
Reads supporting (aligned to +/- strand): ref base C (24/45); new base T (29/18); total (53/63)
Fisher's exact test for biased strand distribution p-value = 4.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 123 total aligned reads displayed.

GAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGT  >  minE/2795739‑2795865

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                        tctGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTAtc                                    <  1:593750/69‑1 (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:912675/1‑62 (MQ=255)

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

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

                                  aTGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTg                           <  1:3000537/68‑1 (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:347769/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                      gCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTAc                      >  1:2377205/1‑69 (MQ=255)

                                                       aaGATATGGGATCCGCTGCGGTGCTGATTGGCATTATCGTCGCCGTg                           <  1:377848/47‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:1441097/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:889734/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:2283643/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:1238897/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:563909/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:1936169/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:526360/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:459241/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:1280152/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:2394421/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:311242/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:3206080/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:2986461/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:2966693/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:2921616/67‑1 (MQ=255)

                                                         gATATGGGATCTGCCGCAGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATCCTGTTAt     <  1:3381476/67‑1 (MQ=255)

                                                          atatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGt  <  1:1610263/69‑1 (MQ=255)

                                                          atatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGt  <  1:11572/69‑1 (MQ=255)

                                                          atatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGt  <  1:2276781/69‑1 (MQ=255)

                                                          atatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGt  <  1:2850272/69‑1 (MQ=255)

                                                          atatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTCGTGCATTCTGTTATGGt  <  1:3198317/69‑1 (MQ=255)

                                                           tatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGt  <  1:1413769/68‑1 (MQ=255)

                                                             tGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGtt       >  1:1818485/1‑61 (MQ=255)

                                                             tGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGtt       >  1:2585929/1‑61 (MQ=255)

                                                             tGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGtt       >  1:2678533/1‑61 (MQ=255)

                                                             tGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGtt       >  1:2250697/1‑61 (MQ=255)

                                                             tGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGtt       >  1:160411/1‑61 (MQ=255)

GAAGCAGTGGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGT  >  minE/2795739‑2795865

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`