BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,863	C→A	22.8%	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	22.8%	1466.7 / inf	664	K96N (AAG→AAT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (268/244); new base A (67/85); total (335/329)
Fisher's exact test for biased strand distribution p-value = 7.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 687 total aligned reads displayed.

CACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTT  >  CP009273/4106796‑4106930

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

     tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGc                                                                 <  1:5848459/67‑1 (MQ=255)

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

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

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

          aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            >  1:10673774/1‑67 (MQ=255)

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

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

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

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

           tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGtt                                                            >  1:10421550/1‑66 (MQ=255)

            gCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCt                                                         <  1:13505658/68‑1 (MQ=255)

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

                gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCttttt                                                     <  1:14017151/68‑1 (MQ=255)

                    aTTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCa                                                 >  1:26594537/1‑68 (MQ=255)

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

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

                     tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAg                                                <  1:10163221/68‑1 (MQ=255)

                      ttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAg                                                >  1:29341701/1‑67 (MQ=255)

                      ttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGa                                               <  1:7012376/68‑1 (MQ=255)

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

                        cTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACa                                             <  1:13222150/68‑1 (MQ=255)

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

                            ggTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAg                                         >  1:4140177/1‑68 (MQ=255)

                            ggTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAg                                         >  1:26218761/1‑68 (MQ=255)

                            ggTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAg                                         >  1:20780548/1‑68 (MQ=255)

                              tACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGt                                        <  1:21790240/67‑1 (MQ=255)

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

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

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

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

                                cgacgGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                     <  1:27802808/68‑1 (MQ=255)

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

                                    ggCACTGCCATACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCa                                 <  1:15909333/68‑1 (MQ=255)

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

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

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

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

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

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

                                      cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCAc                                <  1:30015794/67‑1 (MQ=255)

                                       aCTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTCGGTTCAc                                >  1:5154615/1‑66 (MQ=255)

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

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

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

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

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

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

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

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

                                            cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGt                          >  1:2691337/1‑67 (MQ=255)

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

                                               aCAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      >  1:21134855/1‑68 (MQ=255)

                                               aCAATGGCGGTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      >  1:5596289/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                cAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGt                      <  1:27497381/67‑1 (MQ=255)

                                                    ggCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc                 <  1:3557385/68‑1 (MQ=255)

                                                    ggCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATa                  >  1:17732139/1‑67 (MQ=255)

                                                    ggCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATa                  >  1:3344122/1‑67 (MQ=255)

                                                        ttATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc                 <  1:18665206/64‑1 (MQ=255)

                                                           tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATa            >  1:25886224/1‑66 (MQ=255)

                                                           tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGc          <  1:22871940/68‑1 (MQ=255)

                                                           tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGc          <  1:6544699/68‑1 (MQ=255)

                                                           tagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGc          <  1:9138437/68‑1 (MQ=255)

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

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

                                                            agatagGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATa            >  1:22538514/1‑65 (MQ=255)

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

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

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

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

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

                                                              atagGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGc       >  1:31885739/1‑68 (MQ=255)

                                                                agGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCaa     <  1:28544161/68‑1 (MQ=255)

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

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

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

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

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

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

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATACCTGCAAt    >  1:4386606/1‑68 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCACAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:15564213/1‑68 (MQ=255)

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

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

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

CACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTT  >  CP009273/4106796‑4106930

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`