BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,795,810	T→C	40.5%	A100A (GCT→GCC)	dgkA →	diacylglycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,795,810	0	T	C	40.5%	45.1 / 84.5	116	A100A (GCT→GCC)	dgkA	diacylglycerol kinase
Reads supporting (aligned to +/- strand): ref base T (24/45); new base C (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 = 9.00e-01

Only 100 of 116 total aligned reads displayed.

GGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGT  >  minE/2795747‑2795865

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

       cGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGAtt                                             <  1:2702086/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:1586865/69‑1 (MQ=255)

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

                    aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCAtt                                       >  1:912675/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:1068364/1‑62 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                            taGCTTTCCGGACGCGCTAAAGATATGGGATATGCCGCAGTGCTGATTGc                                           >  1:3258656/2‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

GGTTGACCGAATTGGCTCTGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGT  >  minE/2795747‑2795865

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