BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,586	G→T	84.0%	R189S (CGT→AGT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,586	0	G	T	84.0%	1514.4 / 248.0	548	R189S (CGT→AGT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (45/43); new base T (237/223); total (282/266)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 553 total aligned reads displayed.

AATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCA  >  CP009273/4106521‑4106653

    tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGc                                                               <  1:27947537/68‑1 (MQ=255)

    tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGc                                                               <  1:1664766/68‑1 (MQ=255)

    tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGc                                                               <  1:29082266/68‑1 (MQ=255)

    tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGc                                                               <  1:11303192/68‑1 (MQ=255)

      cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGt                                                             >  1:10235292/1‑68 (MQ=255)

      cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGt                                                             >  1:12531706/1‑68 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTTTGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                                                         <  1:24155044/68‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTgg                                                          <  1:8341387/67‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTgg                                                          <  1:20639830/67‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTgg                                                          <  1:213803/67‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTg                                                           <  1:24617654/66‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTg                                                           <  1:6830282/66‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGgt                                                         <  1:24398060/68‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGgt                                                         <  1:28247649/68‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGgt                                                         >  1:10602128/1‑68 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTgg                                                          <  1:14659486/67‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTAATAGAGGCGTTGgt                                                         >  1:4257649/1‑68 (MQ=255)

          aCTTCCAGCATTTTGTCGACCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGgt                                                         >  1:27479530/1‑68 (MQ=255)

                  cATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCgg                                                 >  1:26552680/1‑68 (MQ=255)

                  cATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTCTAATCgg                                                 >  1:24712088/1‑68 (MQ=255)

                   aTTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCgg                                                 >  1:4083704/1‑67 (MQ=255)

                   aTTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGt                                                <  1:3146383/68‑1 (MQ=255)

                     tttGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTca                                              <  1:5565235/68‑1 (MQ=255)

                     tttGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTca                                              <  1:333340/68‑1 (MQ=255)

                       tggcgtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGGTGGTGTAATCGGTcaca                                            <  1:22691387/65‑1 (MQ=255)

                       tGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTcaca                                            <  1:11798330/68‑1 (MQ=255)

                       tGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTcaca                                            <  1:14982882/68‑1 (MQ=255)

                       tGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTcaca                                            <  1:19037831/68‑1 (MQ=255)

                          cgtcCCAGTCCAGGGTTTGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATgg                                         >  1:17047060/1‑68 (MQ=255)

                             cccAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGac                                       <  1:16380801/67‑1 (MQ=255)

                             cccAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGac                                       <  1:27697990/67‑1 (MQ=255)

                             cccAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGac                                       >  1:18664762/1‑67 (MQ=255)

                                   ccAGGTTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  <  1:26781923/66‑1 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACTAGAGGGGTTGGTGTAATCGGTCACATGGACACGGc                                  <  1:23306668/66‑1 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  >  1:17743823/1‑66 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  >  1:221217/1‑66 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  >  1:7316077/1‑66 (MQ=255)

                                   acAGGGTATGGATGTTGAACAACAAGGTACTAGAGGCGTTTGTGTAATCGGTCACATGGACACGGc                                  <  1:4670149/65‑1 (MQ=255)

                                       ggTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:18805190/1‑68 (MQ=255)

                                       ggTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:23013650/1‑68 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTGCTTGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCGTGAg                            >  1:20144287/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGa                             <  1:10859561/66‑1 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGa                             <  1:3635144/66‑1 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:26645544/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:2361069/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:2617092/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:28434847/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:28811700/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:2927704/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:9504522/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:22553214/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:10161304/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:1205467/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:20243418/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:14148649/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:15421079/1‑67 (MQ=255)

                                            ggATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtt                       >  1:20779287/1‑68 (MQ=255)

                                             gATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                      >  1:9314027/1‑68 (MQ=255)

                                             gATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                      >  1:6306175/1‑68 (MQ=255)

                                             gATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                      >  1:13732611/1‑68 (MQ=255)

                                             gATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                      >  1:16446792/1‑68 (MQ=255)

                                             gATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                      >  1:2862243/1‑68 (MQ=255)

                                              aTGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                      <  1:4970928/67‑1 (MQ=255)

                                              aTGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                      <  1:19925004/67‑1 (MQ=255)

                                               tGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                     >  1:9589276/1‑67 (MQ=255)

                                                  tGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAg                 >  1:26963072/1‑68 (MQ=255)

                                                   gAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGa                <  1:26360863/68‑1 (MQ=255)

                                                     acaacaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa              >  1:17997056/1‑68 (MQ=255)

                                                     acaacaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa              >  1:8793980/1‑68 (MQ=255)

                                                     acaacaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa              >  1:27275098/1‑68 (MQ=255)

                                                     acaacaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCAGAGTCATTTTCCAGATa              >  1:18417362/1‑68 (MQ=255)

                                                       aacaTGGTACTATAGGCGTTGGGGTAATCGGGCACATGGACACGGCCCTGAGTCATTTTCCAGATAAg            <  1:28958250/68‑1 (MQ=255)

                                                       aacaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATaa             >  1:16047381/1‑67 (MQ=255)

                                                        acaTGGTACTAGAGGCGTTTGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGc           <  1:9102853/68‑1 (MQ=255)

                                                        acaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGc           <  1:25558542/68‑1 (MQ=255)

                                                        acaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGc           <  1:13301360/68‑1 (MQ=255)

                                                        acaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGc           <  1:506293/68‑1 (MQ=255)

                                                          aTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGcc          >  1:25785498/1‑67 (MQ=255)

                                                          aTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGcc          >  1:5800862/1‑67 (MQ=255)

                                                          aTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGcc          >  1:19688643/1‑67 (MQ=255)

                                                          aTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGcc          >  1:15576390/1‑67 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:13385098/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:19394620/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:11565424/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:21699469/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:13105886/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:6528881/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:1401341/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:6915357/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:12393323/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:17373633/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTGGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:8051225/68‑1 (MQ=255)

                                                            ggTACTAGAGGCGTCGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:1584541/68‑1 (MQ=255)

                                                            ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:12670024/68‑1 (MQ=255)

                                                            ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:25339251/68‑1 (MQ=255)

                                                            ggTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:540324/68‑1 (MQ=255)

                                                             gTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGt      >  1:7744237/1‑68 (MQ=255)

                                                              tACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGt      <  1:29566318/67‑1 (MQ=255)

                                                              tACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGt      <  1:21883931/67‑1 (MQ=255)

                                                               aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa     >  1:17582118/1‑67 (MQ=255)

AATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCA  >  CP009273/4106521‑4106653

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