BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,863	C→A	6.9%	K96N (AAG→AAT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,863	0	C	A	6.9%	2962.1 / 136.4	842	K96N (AAG→AAT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (405/376); new base A (29/31); total (434/408)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 880 total aligned reads displayed.

CACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTT  >  CP009273/4106796‑4106930

 aCGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCt                                                                    >  1:8489234/1‑68 (MQ=255)

 aCGTTTTAAATCCTCGCAGATTTCTGAGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGc                                                                     <  1:39088482/67‑1 (MQ=255)

   gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCtt                                                                   >  1:35099349/1‑67 (MQ=255)

   gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCtt                                                                   >  1:35794800/1‑67 (MQ=255)

     tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGc                                                                 <  1:45400589/67‑1 (MQ=255)

        aaaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGt                                                             <  1:37912170/68‑1 (MQ=255)

         aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:11069881/1‑68 (MQ=255)

         aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:26432426/1‑68 (MQ=255)

         aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:34640237/1‑68 (MQ=255)

         aaTGCTCGCAGATTTCTGCGGCACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:48359792/1‑68 (MQ=255)

           tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:30673481/1‑66 (MQ=255)

           tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:28031999/1‑66 (MQ=255)

           tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:17130171/1‑66 (MQ=255)

           tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:12854039/1‑66 (MQ=255)

                gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCttttt                                                     <  1:7308688/68‑1 (MQ=255)

                  aGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTc                                                    <  1:41285727/67‑1 (MQ=255)

                  aGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTc                                                    >  1:23354916/1‑67 (MQ=255)

                     tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTccc                                                  <  1:40526358/66‑1 (MQ=255)

                     tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAg                                                <  1:13583815/68‑1 (MQ=255)

                     tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAg                                                <  1:23699931/68‑1 (MQ=255)

                     tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAg                                                <  1:21545041/68‑1 (MQ=255)

                     tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAg                                                <  1:29852794/68‑1 (MQ=255)

                     ttacTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAg                                                <  1:16574139/65‑1 (MQ=255)

                        cTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACa                                             <  1:446742/68‑1 (MQ=255)

                        cTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACa                                             <  1:11825957/68‑1 (MQ=255)

                        cTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACa                                             <  1:10601952/68‑1 (MQ=255)

                        cTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACa                                             <  1:39978756/68‑1 (MQ=255)

                          gCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAAt                                           <  1:30374406/68‑1 (MQ=255)

                          gCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAAt                                           >  1:6744903/1‑68 (MQ=255)

                          gCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAAt                                           <  1:16213804/68‑1 (MQ=255)

                            ggTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATa                                          <  1:4634236/67‑1 (MQ=255)

                            ggTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATa                                          <  1:48990132/67‑1 (MQ=255)

                               acgacgGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTgg                                      >  1:42356796/1‑68 (MQ=255)

                                cgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                     <  1:47904339/68‑1 (MQ=255)

                                cgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTGCTTTTTCCCAGACAATAGTGGt                                     <  1:5949521/68‑1 (MQ=255)

                                 gacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGtt                                    >  1:8967554/1‑68 (MQ=255)

                                    ggCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTc                                  >  1:32606101/1‑67 (MQ=255)

                                    ggCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTc                                  >  1:45867819/1‑67 (MQ=255)

                                    ggCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTc                                  >  1:18155681/1‑67 (MQ=255)

                                       aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACg                               <  1:33297390/67‑1 (MQ=255)

                                        cTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCt                             <  1:39590021/68‑1 (MQ=255)

                                          gCCAGACAATGGCGTTATAGATAGGCTTGCCTGTTTCTTTTTCCCAGAAAATAGTGGTTTCACGCTgg                           <  1:38801322/68‑1 (MQ=255)

                                          gCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                           <  1:15900973/68‑1 (MQ=255)

                                          gCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                           <  1:32116425/68‑1 (MQ=255)

                                          gCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                           <  1:5701963/68‑1 (MQ=255)

                                          gCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                           <  1:6229461/68‑1 (MQ=255)

                                          gCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                           <  1:18596811/68‑1 (MQ=255)

                                               aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      >  1:22676902/1‑68 (MQ=255)

                                               aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      >  1:51316313/1‑68 (MQ=255)

                                               aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      >  1:21109892/1‑68 (MQ=255)

                                               aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      >  1:2036038/1‑68 (MQ=255)

                                               aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      >  1:13332995/1‑68 (MQ=255)

                                               aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      >  1:48646878/1‑68 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:36175739/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:40838993/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:33883245/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:29491164/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:4974641/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:27106431/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:7912052/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:22953632/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:19719721/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:49279245/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:17688052/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:14350833/65‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      <  1:28795496/67‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTa                     <  1:28397534/68‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTa                     <  1:12389309/68‑1 (MQ=255)

                                                cAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:22025261/65‑1 (MQ=255)

                                                 aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTaa                    >  1:36896406/1‑68 (MQ=255)

                                                  aTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAAt                   >  1:29191197/1‑68 (MQ=255)

                                                   tGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATa                  <  1:50319124/68‑1 (MQ=255)

                                                    ggCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc                 <  1:19000814/68‑1 (MQ=255)

                                                     gCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAcc                >  1:45047938/1‑68 (MQ=255)

                                                     gCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAcc                >  1:14038430/1‑68 (MQ=255)

                                                       gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGa              >  1:10004117/1‑68 (MQ=255)

                                                       gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTATCACGCTGGTTCGTAATACCGa              >  1:23948715/1‑68 (MQ=255)

                                                            agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATa            >  1:38230115/1‑65 (MQ=255)

                                                            agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATa            >  1:4356966/1‑65 (MQ=255)

                                                            agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATa            >  1:26812954/1‑65 (MQ=255)

                                                            agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATa            >  1:24724715/1‑65 (MQ=255)

                                                              atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:9494426/1‑68 (MQ=255)

                                                              atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:8162590/1‑68 (MQ=255)

                                                              atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:27690569/1‑68 (MQ=255)

                                                              atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:10226613/1‑68 (MQ=255)

                                                              atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:26148611/1‑68 (MQ=255)

                                                              atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:47330739/1‑68 (MQ=255)

                                                              atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:41851434/1‑68 (MQ=255)

                                                              atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:44099690/1‑68 (MQ=255)

                                                               tagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCa      <  1:4298019/68‑1 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:50676545/1‑68 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:37147756/1‑68 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:25700506/1‑68 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:7707865/1‑68 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:25442795/1‑68 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:14949703/1‑68 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:46800294/1‑68 (MQ=255)

                                                                  gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTc                        <  1:11518545/47‑1 (MQ=255)

                                                                  gCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAtt   >  1:15463102/1‑68 (MQ=255)

                                                                   cTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAttt  >  1:3119996/1‑68 (MQ=255)

CACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTT  >  CP009273/4106796‑4106930

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