BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,955	C→T	9.3%	A66T (GCG→ACG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,955	0	C	T	9.3%	3067.5 / 229.1	997	A66T (GCG→ACG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (420/484); new base T (51/42); total (471/526)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1009 total aligned reads displayed.

AGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAC  >  CP009273/4106890‑4107020

aGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCa                                                                 >  1:41418720/1‑68 (MQ=255)

 gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCa                                                                 >  1:48263876/1‑67 (MQ=255)

 gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCa                                                                 >  1:17221299/1‑67 (MQ=255)

   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGc                                                               >  1:24253407/1‑67 (MQ=255)

   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGc                                                               >  1:20055249/1‑67 (MQ=255)

   ggTTTCACGCTGGTTCGTAATACCGATAGCAGCAATTTGATCGGAACTGATATCGGCAATCGCCAGc                                                               >  1:10597896/1‑67 (MQ=255)

     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                             <  1:32405580/67‑1 (MQ=255)

     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                             <  1:5870650/67‑1 (MQ=255)

     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                             <  1:36918244/67‑1 (MQ=255)

     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACt                                                            >  1:22771741/1‑68 (MQ=255)

     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                             <  1:13729414/67‑1 (MQ=255)

     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                             <  1:41037930/67‑1 (MQ=255)

     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACt                                                            >  1:23422373/1‑68 (MQ=255)

     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCACCAGCAc                                                             <  1:26389049/67‑1 (MQ=255)

     ttgcACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCAc                                                             <  1:253473/64‑1 (MQ=255)

      ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACtt                                                           <  1:48712373/68‑1 (MQ=255)

           gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAcc                                                      >  1:3514469/1‑68 (MQ=255)

            cTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCa                                                     >  1:4706315/1‑68 (MQ=255)

              ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGc                                                   <  1:28876698/68‑1 (MQ=255)

              ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGc                                                   <  1:28592901/68‑1 (MQ=255)

              ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGc                                                   <  1:14363535/68‑1 (MQ=255)

                ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGt                                                 >  1:50865638/1‑68 (MQ=255)

                 tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTg                                                >  1:38739645/1‑68 (MQ=255)

                   gTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                              >  1:37728736/1‑68 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                              >  1:10051672/1‑66 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATCTCGGCTTTCGCCAGCAAATCTACCAGCGTGGa                                              >  1:18935183/1‑66 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTggg                                              >  1:7241482/1‑65 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                              >  1:18672675/1‑66 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                              >  1:16868328/1‑66 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGa                                              >  1:50079020/1‑66 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAg                                             <  1:49688330/67‑1 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                            <  1:37260354/68‑1 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                            <  1:30133033/68‑1 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                            <  1:36054233/68‑1 (MQ=255)

                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                            <  1:21649598/68‑1 (MQ=255)

                     aaTACCGATAGCGGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGc                                            <  1:15924603/68‑1 (MQ=255)

                      aTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCt                                           >  1:27619989/1‑68 (MQ=255)

                       tACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCtt                                          >  1:1108356/1‑68 (MQ=255)

                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                                         >  1:40030794/1‑68 (MQ=255)

                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                                         >  1:29623379/1‑68 (MQ=255)

                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                                         >  1:3185699/1‑68 (MQ=255)

                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                                         >  1:2817942/1‑68 (MQ=255)

                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCttt                                         >  1:14085849/1‑68 (MQ=255)

                          cGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                       >  1:4292819/1‑68 (MQ=255)

                             tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTg                                    >  1:22354513/1‑68 (MQ=255)

                             tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCACCAGCACTTCTACCAGCGTGGAGCTTTGGGTg                                    >  1:45329015/1‑68 (MQ=255)

                               gCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                   >  1:47453889/1‑67 (MQ=255)

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

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

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

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

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

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

                                    aaTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGa                             <  1:43478521/68‑1 (MQ=255)

                                      tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAtt                           <  1:33665003/68‑1 (MQ=255)

                                         gaacGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                        <  1:6655410/65‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGt                 <  1:19241542/68‑1 (MQ=255)

                                                   gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGTCCCAGATTTCCATTGGGTCgt              <  1:46438181/68‑1 (MQ=255)

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

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

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

                                                   gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCgt              <  1:44180067/68‑1 (MQ=255)

                                                   gATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTCCATTGGGTCgt              <  1:15076579/68‑1 (MQ=255)

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

                                                   gATATCGGCTTTCACCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt              <  1:29229325/68‑1 (MQ=255)

                                                     tctCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt            <  1:30860514/66‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                        cGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt         >  1:39351231/1‑68 (MQ=255)

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

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

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

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

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

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

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

                                                           cTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGTTGGCCCAGATTTCCATTGGGTCGTGTTCTAc       <  1:11537024/67‑1 (MQ=255)

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

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

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

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

AGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAC  >  CP009273/4106890‑4107020

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