BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,988	C→T	7.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	7.2%	2409.2 / 115.2	739	A55T (GCC→ACC) ‡	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (308/376); new base T (25/28); total (333/406)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 755 total aligned reads displayed.

CTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGC  >  CP009273/4106922‑4107055

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

  gCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGTCCa                                                                  <  1:21266630/68‑1 (MQ=255)

    aaTTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                                <  1:41239778/68‑1 (MQ=255)

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

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

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

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

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

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

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

    aaTTTGATCGGAACTGATATCGGCGTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                                <  1:24399352/68‑1 (MQ=255)

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

        tGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTc                                                            <  1:3396635/68‑1 (MQ=255)

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

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

         gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTcc                                                           >  1:48658516/1‑68 (MQ=255)

          aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                                                          >  1:47163162/1‑68 (MQ=255)

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

                aCTGATCTCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                    <  1:8987731/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                   gATATCGGCTTTCGCCAGCACTGCTCCCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCgt                                                 <  1:48974887/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                     tatCGGCTTTCGCCAGCACTTCTACCAGCGGGGAGCTTTGGGTGGTCCAGATTTCCATTGGGGCgtgt                                               <  1:42344603/68‑1 (MQ=255)

                         ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                           <  1:2589799/68‑1 (MQ=255)

                         ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCTa                                           <  1:19234531/68‑1 (MQ=255)

                          gcgttCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                           >  1:11427418/4‑67 (MQ=255)

                          gCTTTCGCCAGCACTTCTACCAGCGTGGGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                           >  1:18736461/1‑67 (MQ=255)

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

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

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

                          gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCTa                                           >  1:34562869/1‑67 (MQ=255)

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

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

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

                           cTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCt                                            <  1:311733/65‑1 (MQ=255)

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

                               cgtcAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAc                                     <  1:43940107/65‑1 (MQ=255)

                               cgtcAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAc                                     <  1:36636644/65‑1 (MQ=255)

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

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

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

                                      aCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTgg                                 <  1:30829001/65‑1 (MQ=255)

                                      aCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGttt                              >  1:25417746/1‑68 (MQ=255)

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

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

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

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

                                         tCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                           <  1:49736946/68‑1 (MQ=255)

                                            aCCAGCGTGGAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTa                        >  1:51279273/1‑68 (MQ=255)

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

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

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

                                            aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGATCTACCCAACCTGGTTTTGGGTa                        >  1:45092960/1‑68 (MQ=255)

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

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

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

                                              cAGCGTGGGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                      >  1:23450019/1‑68 (MQ=255)

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

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

                                              cAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                      >  1:37864883/1‑68 (MQ=255)

                                               aGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAt                     >  1:40547641/1‑68 (MQ=255)

                                                 cGTGGAGCTTTGGGTGGCCCAGATTTGCATTGGGTCGTGTTCTACCCAACCTGGGTTTGGGTAGAttt                   >  1:35134348/1‑68 (MQ=255)

                                                 cGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAttt                   >  1:19338300/1‑68 (MQ=255)

                                                    ggAGCTTTGGGTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                 <  1:7313864/67‑1 (MQ=255)

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

                                                       gCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaa             >  1:13188398/1‑68 (MQ=255)

                                                       gCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaa             >  1:35437103/1‑68 (MQ=255)

                                                       gCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaa             >  1:39203239/1‑68 (MQ=255)

                                                       gCTATGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaa             >  1:38728399/1‑68 (MQ=255)

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

                                                            gggTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTc         >  1:32151400/1‑67 (MQ=255)

                                                             ggTGGCCAAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAATTGCTCAAATTagc       >  1:31527316/1‑68 (MQ=255)

                                                              gTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgc       >  1:19598079/1‑67 (MQ=255)

                                                              gTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCACCCTGGTTTTGGGTAGATTTGCTCAAATTcgc       >  1:5958949/1‑67 (MQ=255)

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

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

                                                                ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCt    <  1:50443855/68‑1 (MQ=255)

                                                                ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCt    <  1:15796782/68‑1 (MQ=255)

                                                                ggCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGGTTTGGGTAGATTTGCTCAAATTCGCGCt    <  1:12809280/68‑1 (MQ=255)

                                                                ggCCCAGATTTCCATTGGGGCGTGTTCTACCCAACCTGCTTATGGGTAGATTTGCTCAAATTCGCGCt    <  1:6732208/68‑1 (MQ=255)

                                                                 gCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCt    >  1:23086347/1‑67 (MQ=255)

CTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGC  >  CP009273/4106922‑4107055

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