BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	229,412	G→T	48.3%	T141T (ACG→ACT)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	229,412	0	G	T	48.3%	27.2 / 130.3	120	T141T (ACG→ACT)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand): ref base G (25/37); new base T (16/42); total (41/79)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 120 total aligned reads displayed.

CACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGA  >  minE/229346‑229477

cacCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGc                                                                 <  1:1626390/69‑1 (MQ=255)

cacCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGc                                                                 <  1:2221313/69‑1 (MQ=255)

cacCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGc                                                                 <  1:633669/69‑1 (MQ=255)

 acCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGc                                                                 <  1:2844158/68‑1 (MQ=255)

   ccTGCGAGCGTTGCTCGATAACAATATCGTTGCGGTAATCAATGAGAACGATGCTGTCGCTACGGCaga                                                              <  1:3145865/69‑1 (MQ=255)

   ccTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCaga                                                              <  1:519782/69‑1 (MQ=255)

   ccTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCaga                                                              <  1:2852899/69‑1 (MQ=255)

   ccTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCaga                                                              <  1:337726/69‑1 (MQ=255)

   ccTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCaga                                                              <  1:1534993/69‑1 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGtag                                                   >  1:1245585/1‑66 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:24169/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:1967120/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:3306457/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:1534511/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:1135208/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:2983833/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:2281521/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:1075647/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:1642855/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:1702012/1‑68 (MQ=255)

               gCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCg                                                   >  1:3102856/1‑68 (MQ=255)

                 tCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGc                                                 <  1:555387/68‑1 (MQ=255)

                 tCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGc                                                 <  1:456772/68‑1 (MQ=255)

                 tCGATAACAATATCGTTCCCGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGc                                                 <  1:2667654/68‑1 (MQ=255)

                  cGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCgg                                                  >  1:1687277/1‑66 (MQ=255)

                  cGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcg                                                >  1:2163933/1‑68 (MQ=255)

                      aaCAATATCGTTCCGGTAATTATTGAGAACGATGCTGTCGCTACTGCAGAGATTATGGTCg                                                   <  1:3333587/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:361748/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:844284/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:2731624/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:2161552/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:3181681/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:2257800/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:665447/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:1387314/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:48958/61‑1 (MQ=255)

                      aaCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:1570725/61‑1 (MQ=255)

                       aCAATATCGTTCCGGTAATTAATGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCg                                                   <  1:322861/60‑1 (MQ=255)

                       aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                              <  1:2504043/65‑1 (MQ=255)

                       aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                              <  1:2150652/65‑1 (MQ=255)

                       aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                              <  1:2454528/65‑1 (MQ=255)

                       aCAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgat                                              <  1:162178/65‑1 (MQ=255)

                           taacGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                       <  1:2633523/65‑1 (MQ=255)

                            atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataactata                                       <  1:1102225/67‑1 (MQ=255)

                            atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                       <  1:2064955/67‑1 (MQ=255)

                            atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                       <  1:2049991/67‑1 (MQ=255)

                            atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                       <  1:3242908/67‑1 (MQ=255)

                            atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                       <  1:48655/67‑1 (MQ=255)

                            atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                       <  1:1076240/67‑1 (MQ=255)

                            atCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                       <  1:855063/67‑1 (MQ=255)

                               gTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCtt                                  >  1:346607/1‑69 (MQ=255)

                                     gTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGcgataacgata                                       <  1:2036773/58‑1 (MQ=255)

                                       aaTCAATGAGAACGATGCTGTCGCTACGGCAGAGAt                                                           >  1:1496974/1‑36 (MQ=255)

                                       aaTCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCt                          >  1:1061109/1‑69 (MQ=255)

                                       aaTCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCt                          >  1:2520103/1‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGtgataa                                            >  1:3299447/3‑49 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:570630/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:2141648/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:1974659/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:235567/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:1934235/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:2908842/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:1761417/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:1732580/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:2786599/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:583756/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:1276936/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:1012135/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:2529399/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTgg                        >  1:2566539/3‑69 (MQ=255)

                                         ttaaTGAGAACGATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTg                         >  1:940989/3‑68 (MQ=255)

                                                    gATGCTGTCGCTTCTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:307558/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCTATTCtg             <  1:36493/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:2427776/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:2615371/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:934070/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1076215/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:766102/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1116829/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:2733478/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:2776216/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:281555/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1351477/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1489476/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:245115/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:3160617/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1783840/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1814150/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1834553/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:2926836/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:3258907/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1969324/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:225412/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:3192344/69‑2 (MQ=255)

                                                    gATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtt             >  1:14850/1‑69 (MQ=255)

                                                    gATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtt             >  1:3355163/1‑69 (MQ=255)

                                                    gATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtt             >  1:2507037/1‑69 (MQ=255)

                                                     aTGCTGTCGCTACTGCAGAGATTAAGGTCGGTGATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:1893231/68‑2 (MQ=255)

                                                            cGCTACTGCAGAGATTAAGGTCGGTTATAACGATAACCTTTCTGCGCTGGCGGCGATTCtg             <  1:645252/61‑2 (MQ=255)

                                                                aCGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGa  >  1:22259/1‑68 (MQ=255)

CACCCTGCGAGCGTTGCTCGATAACAATATCGTTCCGGTAATCAATGAGAACGATGCTGTCGCTACGGCAGAGATTAAGGTCGGCGATAACGATAACCTTTCTGCGCTGGCGGCGATTCTTGCGGGTGCCGA  >  minE/229346‑229477

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