BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,863	C→A	41.6%	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	41.6%	396.2 / inf	560	K96N (AAG→AAT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (163/164); new base A (114/119); total (277/283)
Fisher's exact test for biased strand distribution p-value = 8.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 584 total aligned reads displayed.

CACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTT  >  CP009273/4106796‑4106930

       tAAATTCTCGCAGATTTCTGCGGTACTTCGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCgg                                                              <  1:3518562/68‑1 (MQ=255)

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

         aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGtt                                                            <  1:1678709/68‑1 (MQ=255)

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

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

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

         aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGtt                                                            >  1:3939436/1‑68 (MQ=255)

         aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGtt                                                            >  1:13458551/1‑68 (MQ=255)

          aTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCggct                                                            >  1:17837710/1‑65 (MQ=255)

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

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

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

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

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

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

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

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

                      ttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGa                                               <  1:9417955/68‑1 (MQ=255)

                       tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGAc                                              <  1:17448140/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                          gCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATACCATTGCCGGTATCTTTTTCCCAGACAAt                                           >  1:16019352/1‑68 (MQ=255)

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

                               acgacgGCACTGCCAGACAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTgg                                      <  1:12430954/68‑1 (MQ=255)

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

                                   cgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTc                                  <  1:5237034/68‑1 (MQ=255)

                                   cgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTc                                  <  1:2305353/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                            cAGACAATGGCGTTATCGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                           <  1:31789735/66‑1 (MQ=255)

                                            cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGt                          <  1:31147975/67‑1 (MQ=255)

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

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

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

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

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

                                               aCAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGAGGTTTCACGCTGGTTCGt                      >  1:9039264/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                cAATGGCGTTATAGATAGGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTa                     <  1:13915728/68‑1 (MQ=255)

                                                    ggCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATa                  >  1:10538879/1‑67 (MQ=255)

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

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

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

                                                      cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCg               >  1:10923149/1‑68 (MQ=255)

                                                      cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCg               >  1:18690783/1‑68 (MQ=255)

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

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

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

                                                            agatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAAAACCGATa            >  1:11612755/1‑65 (MQ=255)

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

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

                                                            agatagGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCt         <  1:21160479/68‑1 (MQ=255)

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

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

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

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

                                                                agGATTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCaa     <  1:30568142/68‑1 (MQ=255)

                                                                 ggCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCGCGCTGGTTCGTAATACCGATAGCTGCAAt    >  1:23679451/1‑68 (MQ=255)

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

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

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

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

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

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

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

                                                                   cTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAttt  >  1:7962283/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: `‑`