BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,988	C→T	6.2%	A55T (GCC→ACC) ‡	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,988	0	C	T	6.2%	1265.2 / 49.7	383	A55T (GCC→ACC) ‡	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (174/185); new base T (14/10); total (188/195)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 399 total aligned reads displayed.

TGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCT  >  CP009273/4106923‑4107053

tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                                                   <  1:27354316/66‑1 (MQ=255)

tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                                                   <  1:1174307/66‑1 (MQ=255)

tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                                                   <  1:3333758/66‑1 (MQ=255)

tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGAcc                                                                  <  1:24396495/67‑1 (MQ=255)

   aaTTTGATCGGAACTGATCTCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                              <  1:26844569/68‑1 (MQ=255)

   aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGa                                                              <  1:4922200/68‑1 (MQ=255)

   aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                              <  1:1097082/68‑1 (MQ=255)

   aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                              <  1:25223405/68‑1 (MQ=255)

   aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                              <  1:11644906/68‑1 (MQ=255)

   aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                              <  1:6697618/68‑1 (MQ=255)

   aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                              <  1:14008098/68‑1 (MQ=255)

    aTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAt                                                             >  1:20066169/1‑68 (MQ=255)

    aTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAt                                                             <  1:5561650/68‑1 (MQ=255)

     tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGAtt                                                            >  1:19639421/1‑68 (MQ=255)

     tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAtt                                                            >  1:15768134/1‑68 (MQ=255)

     tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAtt                                                            >  1:23255385/1‑68 (MQ=255)

        gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGATCTTTGGGTGGCCCAGATTTcc                                                         <  1:16351027/68‑1 (MQ=255)

        gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                         <  1:8433521/68‑1 (MQ=255)

        gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                         <  1:3663306/68‑1 (MQ=255)

        gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                         <  1:17746436/68‑1 (MQ=255)

        gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                         <  1:1740737/68‑1 (MQ=255)

          tCGGAACTGATCTCGGCTTTCGCCAGCTCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAt                                                       <  1:19250751/68‑1 (MQ=255)

              aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGTTTTGGGTGGCCCAGATTTCCATTgg                                                    <  1:17917307/67‑1 (MQ=255)

              aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                    <  1:21516459/67‑1 (MQ=255)

              aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                    <  1:21670512/67‑1 (MQ=255)

              aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                    <  1:8218590/67‑1 (MQ=255)

               aCTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg                                                   >  1:4411446/1‑67 (MQ=255)

               aCTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                  <  1:22520566/68‑1 (MQ=255)

               aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg                                                   >  1:19725101/1‑67 (MQ=255)

               aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                  <  1:1557992/68‑1 (MQ=255)

               aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTggg                                                   >  1:7954133/1‑67 (MQ=255)

                  gttatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                               <  1:14578666/66‑1 (MQ=255)

                  gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCgt                                               <  1:16269828/68‑1 (MQ=255)

                  gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                               <  1:23041388/68‑1 (MQ=255)

                  gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                               <  1:21794701/68‑1 (MQ=255)

                  gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                               <  1:12529165/68‑1 (MQ=255)

                  gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt                                               <  1:28864212/68‑1 (MQ=255)

                    tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             >  1:17252708/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCCCTTCTACCAGCGTGGAGC‑TTGGGTGACCCAGATTTCCATTGGGTCgtgt                                             >  1:21923469/1‑67 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCgtgt                                             >  1:5889143/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCgtgt                                             >  1:7248652/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:4563671/68‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:27523714/68‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             >  1:2944075/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:28222637/68‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             >  1:10134156/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:20638985/68‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:9602936/68‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             >  1:6089189/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             >  1:16576277/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:1879498/68‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:20308059/68‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtg                                              <  1:2489780/67‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCgtgt                                             >  1:25814165/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCgtgt                                             <  1:17033787/68‑1 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCgtgt                                             >  1:19916630/1‑68 (MQ=255)

                    tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCgtgt                                             <  1:13022089/68‑1 (MQ=255)

                      tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             >  1:21240093/1‑66 (MQ=255)

                       cGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt                                          <  1:2471571/68‑1 (MQ=255)

                         gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                         >  1:26474808/1‑67 (MQ=255)

                         gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                         >  1:14310610/1‑67 (MQ=255)

                         gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                         >  1:8700991/1‑67 (MQ=255)

                         gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                         >  1:25086062/1‑67 (MQ=255)

                         gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                         >  1:1288507/1‑67 (MQ=255)

                         gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                         >  1:18406588/1‑67 (MQ=255)

                         gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                         >  1:21099912/1‑67 (MQ=255)

                          cTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                        <  1:14731227/67‑1 (MQ=255)

                            ttCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCa                                     <  1:13890067/68‑1 (MQ=255)

                             tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaa                                    <  1:9853226/68‑1 (MQ=255)

                              cgtcAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCGCAGATTTCCATTGGGTCGTGTTCTACCCAAc                                   <  1:27067177/65‑1 (MQ=255)

                              cGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAc                                   <  1:9271286/68‑1 (MQ=255)

                              cGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAc                                   <  1:18128083/68‑1 (MQ=255)

                                ccAGCACTTCTACCAGCGTGGATCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaa                                    >  1:18595634/1‑65 (MQ=255)

                                ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaa                                    >  1:13432649/1‑65 (MQ=255)

                                ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAATTCCATTGGGTCGTGTTCTACCCAAcag                                 >  1:6973222/1‑66 (MQ=255)

                                     aCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGt                              <  1:4679478/66‑1 (MQ=255)

                                     aCTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGttt                            >  1:13819565/1‑68 (MQ=255)

                                        tCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                         <  1:11971292/68‑1 (MQ=255)

                                        tCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                         <  1:10603090/68‑1 (MQ=255)

                                           aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTa                      >  1:9006218/1‑68 (MQ=255)

                                            ccAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAg                     >  1:5319115/1‑68 (MQ=255)

                                            ccAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAg                     >  1:16062702/1‑68 (MQ=255)

                                            ccAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAg                     >  1:182886/1‑68 (MQ=255)

                                            ccAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAg                     >  1:13660011/1‑68 (MQ=255)

                                             cAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                    <  1:23465572/68‑1 (MQ=255)

                                             cAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                    >  1:3905195/1‑68 (MQ=255)

                                                cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAttt                 >  1:500237/1‑68 (MQ=255)

                                                cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAttt                 >  1:4856212/1‑68 (MQ=255)

                                                cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAttt                 >  1:18037203/1‑68 (MQ=255)

                                                 gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGt                              >  1:16828056/1‑54 (MQ=255)

                                                 gTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAttt                 >  1:25518384/1‑67 (MQ=255)

                                                   ggAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc               <  1:14985260/67‑1 (MQ=255)

                                                   ggAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc               <  1:23154758/67‑1 (MQ=255)

                                                   ggAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc               <  1:11508942/67‑1 (MQ=255)

                                                     aGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCa            >  1:1138568/1‑68 (MQ=255)

                                                        tttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa          <  1:22088168/67‑1 (MQ=255)

                                                             gTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc     >  1:22829732/1‑67 (MQ=255)

                                                             gTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc     >  1:6262026/1‑67 (MQ=255)

                                                             gTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc     >  1:23814004/1‑67 (MQ=255)

                                                              tGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgcg    <  1:15221607/67‑1 (MQ=255)

TGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCT  >  CP009273/4106923‑4107053

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