BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,457	C→T	16.1%	G232S (GGC→AGC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,457	0	C	T	16.1%	1754.4 / inf	702	G232S (GGC→AGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (236/351); new base T (46/67); total (284/418)
Fisher's exact test for biased strand distribution p-value = 9.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

Only 100 of 782 total aligned reads displayed.

TGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATA  >  CP009273/4106390‑4106524

tGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgcc                                                                     <  1:27282303/68‑1 (MQ=255)

 gACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccg                                                                    <  1:23670616/68‑1 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgctgcc                                                                  >  1:24968849/1‑68 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                  >  1:10313978/1‑68 (MQ=255)

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                  >  1:31062048/1‑68 (MQ=255)

     aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCaa                                                                >  1:22391235/1‑68 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:36630895/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAAt                                                               <  1:426896/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAAt                                                               <  1:14881199/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAAt                                                               <  1:48476740/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:46434961/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               >  1:30408134/1‑68 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:616109/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:21787390/68‑1 (MQ=255)

      aaCAGCGCGGCCGGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAAt                                                               <  1:28585375/68‑1 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATg                                                              >  1:51086494/1‑68 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATg                                                              >  1:26679819/1‑68 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                              >  1:47866911/1‑68 (MQ=255)

           cgcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAg                                                          <  1:39830010/68‑1 (MQ=255)

            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTa                                                           <  1:2528725/66‑1 (MQ=255)

            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGt                                                         <  1:36298358/68‑1 (MQ=255)

                 ctgcttGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                    <  1:32270627/68‑1 (MQ=255)

                 ctgctgTTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:15592626/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                    <  1:4088479/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                    <  1:6760850/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                    <  1:9192487/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGAc                                                    <  1:43489260/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGa                                                     <  1:19935069/67‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:28264978/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:31746985/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:47188454/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:20460910/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:20883637/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:29755815/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:6502537/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:34972475/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:28882458/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:35419570/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:5655528/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:7848654/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:26141861/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:11570691/68‑1 (MQ=255)

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:8540626/68‑1 (MQ=255)

                  tgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGAc                                                    <  1:18453045/67‑1 (MQ=255)

                     tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGta                                                >  1:40299683/1‑68 (MQ=255)

                     tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTGGCCGCCAATGTTAGTCTGACCGt                                                 <  1:45468198/67‑1 (MQ=255)

                      ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGtat                                               >  1:23183441/1‑68 (MQ=255)

                      ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                               >  1:21237406/1‑68 (MQ=255)

                      ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                               >  1:49715027/1‑68 (MQ=255)

                        tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                               >  1:17771788/1‑66 (MQ=255)

                         cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACt                                            <  1:32794206/68‑1 (MQ=255)

                         cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACt                                            <  1:33865276/68‑1 (MQ=255)

                         cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                            <  1:10787489/68‑1 (MQ=255)

                          aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACt                                            >  1:38836117/1‑67 (MQ=255)

                          aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                            <  1:44441344/67‑1 (MQ=255)

                             ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCgg                                        >  1:29029677/1‑68 (MQ=255)

                             ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGCCCGTATACTTCgg                                        >  1:19142647/1‑68 (MQ=255)

                               cGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGga                                       <  1:48940290/67‑1 (MQ=255)

                                  tCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaaga                                    <  1:37578582/67‑1 (MQ=255)

                                   cccGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAac                                  <  1:47502700/68‑1 (MQ=255)

                                      ggAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                               <  1:12454409/68‑1 (MQ=255)

                                        agaTTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:25896983/1‑68 (MQ=255)

                                        agaTTGGAATACGCGTGCCGCGTTTGCCGCCAATGTTAGTCTGACCGTATTCTTCGGAAGAACGACGc                             >  1:43049758/1‑68 (MQ=255)

                                        agaTTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:4399676/1‑68 (MQ=255)

                                        agaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:30054062/1‑68 (MQ=255)

                                        agaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:13774103/1‑68 (MQ=255)

                                         gaTAGGAATACGCGTGCCGCCTTTGCCGCCAAAGTAAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:50053370/1‑67 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                          <  1:10005901/68‑1 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                          <  1:41783791/68‑1 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc                           <  1:24604473/67‑1 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc                           <  1:13002061/67‑1 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                          <  1:16268298/68‑1 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                          <  1:36797531/68‑1 (MQ=255)

                                               agtACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                      >  1:50827673/3‑68 (MQ=255)

                                               aaTACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                      >  1:23905491/1‑68 (MQ=255)

                                               aaTACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                      >  1:38016768/1‑68 (MQ=255)

                                               aaTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                      >  1:2763496/1‑68 (MQ=255)

                                                gtACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                      <  1:44970901/66‑1 (MQ=255)

                                                aTACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                      <  1:14219047/67‑1 (MQ=255)

                                                aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGCCCGTATACTTCGGAAGAACGCCGCACTTCTg                      >  1:29510135/1‑67 (MQ=255)

                                                aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                      <  1:41249672/67‑1 (MQ=255)

                                                  aCGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                   >  1:16641865/1‑68 (MQ=255)

                                                    gcgTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                 >  1:42220933/1‑68 (MQ=255)

                                                    gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGCCGCACTTCTGgcagc                 >  1:35822514/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:31351841/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:5238274/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:9501369/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:7688167/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:7048600/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:46854937/1‑68 (MQ=255)

                                                     cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca                >  1:17242795/1‑68 (MQ=255)

                                                      gTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAt               <  1:39306492/68‑1 (MQ=255)

                                                          cgccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtctc            >  1:8118235/1‑67 (MQ=255)

                                                               ttGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg      <  1:18203830/68‑1 (MQ=255)

                                                               ttGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg      <  1:45868773/68‑1 (MQ=255)

                                                               ttGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg      <  1:34396207/68‑1 (MQ=255)

                                                                tGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa     <  1:49274420/68‑1 (MQ=255)

                                                                tGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa     >  1:48208766/1‑68 (MQ=255)

                                                                   cgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAata  >  1:13238055/1‑68 (MQ=255)

                                                                   cgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAat   <  1:16980655/67‑1 (MQ=255)

TGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATA  >  CP009273/4106390‑4106524

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