BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,987	G→A	6.9%	A55V (GCC→GTC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,987	0	G	A	6.9%	1680.2 / 79.1	511	A55V (GCC→GTC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (196/279); new base A (15/20); total (212/299)
Fisher's exact test for biased strand distribution p-value = 8.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 523 total aligned reads displayed.

AGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCT  >  CP009273/4106920‑4107053

 gCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                    >  1:4474665/1‑67 (MQ=255)

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

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

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

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

  cTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                                                   >  1:29469247/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

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

      aaTTTGAGCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                                                              <  1:28188503/68‑1 (MQ=255)

        tttGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAtt                                                            <  1:19586570/68‑1 (MQ=255)

        tttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGAtt                                                            <  1:32757486/68‑1 (MQ=255)

         ttGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAtt                                                            >  1:16970422/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

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

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

              cGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt                                                      >  1:13268691/1‑68 (MQ=255)

               ggAGCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTg                                                     >  1:35970930/1‑68 (MQ=255)

                gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                    >  1:3226771/1‑68 (MQ=255)

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

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

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

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

                  aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg                                                   <  1:23521364/67‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                       tttCGGGTATCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:30815003/66‑1 (MQ=255)

                       tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:26519469/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                       tatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTGGGGTGGCCCAGATTTCCATTGGGTCgtgt                                             <  1:639314/68‑1 (MQ=255)

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

                          cGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATATCCATTGGGTCGTGttat                                          >  1:13035118/1‑66 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                 cGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGGGTTCTACCCAAc                                   <  1:31789922/68‑1 (MQ=255)

                                      gCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGt                              >  1:22746972/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

                                               ccAGCGTGGATCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAg                     >  1:5556041/1‑68 (MQ=255)

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

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

                                                cAGCGTGGAGCTTTGGGTGGCCCAGATTTCAATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGa                    >  1:28758547/1‑68 (MQ=255)

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

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

                                                  gCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAtt                  >  1:10606645/1‑68 (MQ=255)

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

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

                                                         gCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGtt                                            >  1:4103924/1‑35 (MQ=255)

                                                           tttGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAt         >  1:7360528/1‑68 (MQ=255)

                                                           tttGGGTGGCCCAGATTTCCATTGGGTCGTGTTATACCCAACCTGGTTTTGGGTAGATTTGCTCAAAt         >  1:22936523/1‑68 (MQ=255)

                                                           tttGGGTCGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa          <  1:34630412/67‑1 (MQ=255)

AGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCT  >  CP009273/4106920‑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: `‑`