BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,586	G→T	23.3%	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	23.3%	1236.4 / inf	574	R189S (CGT→AGT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (288/152); new base T (79/54); total (367/207)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 582 total aligned reads displayed.

AATATCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCA  >  CP009273/4106521‑4106653

  tatCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGagag                                                                 <  1:4431891/68‑1 (MQ=255)

    tCCAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGc                                                               <  1:16294412/68‑1 (MQ=255)

     ccAGCTATTCCATCATTTTGTCGTCCCATTCCAGGGTATGGATGTTTAACAACATGGTACGAGAGGCg                                                              <  1:9218491/68‑1 (MQ=255)

      cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGt                                                             >  1:9630293/1‑68 (MQ=255)

      cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGt                                                             >  1:20113693/1‑68 (MQ=255)

      cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGt                                                             >  1:22174198/1‑68 (MQ=255)

      cAGCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTAGAACAACATGGTACTAGAGGCGt                                                             >  1:29140102/1‑68 (MQ=255)

        gCACTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTg                                                           >  1:1359938/1‑68 (MQ=255)

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

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

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

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                                                         >  1:2732405/1‑68 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                                                         <  1:618095/68‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                                                         >  1:8972194/1‑68 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                                                         >  1:1076273/1‑68 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                                                         <  1:20535723/68‑1 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                                                         >  1:13706659/1‑68 (MQ=255)

          aCTTCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGgt                                                         >  1:28510206/1‑68 (MQ=255)

            ttCCAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGgtgt                                                       <  1:1771907/68‑1 (MQ=255)

              ccAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTaa                                                     <  1:27689087/68‑1 (MQ=255)

               cAGCATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTaa                                                     >  1:11593098/1‑67 (MQ=255)

                  cATTTTGTGGTCCCAGTCCAGGGTATGGATGTTGAACAAGATGGTACGAGAGGCGTTGGTGTAATCgg                                                 >  1:5472672/1‑68 (MQ=255)

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

                  cATTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCgg                                                 >  1:24403545/1‑68 (MQ=255)

                   aTTTTGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCgg                                                 <  1:1836357/67‑1 (MQ=255)

                     tttGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTca                                              <  1:16749778/68‑1 (MQ=255)

                     tttGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTca                                              <  1:931146/68‑1 (MQ=255)

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

                       tGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTcaca                                            <  1:4295959/68‑1 (MQ=255)

                       tGTCGTCCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTcaca                                            <  1:8529324/68‑1 (MQ=255)

                        gtcgtcCCAGTCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACAt                                           >  1:4698927/1‑68 (MQ=255)

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

                                  tCCAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  <  1:10171180/67‑1 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  <  1:1571299/66‑1 (MQ=255)

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

                                   ccAGGGTATGGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  <  1:11742307/66‑1 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGcc                                 <  1:27929841/67‑1 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  <  1:19363746/66‑1 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  <  1:25519309/66‑1 (MQ=255)

                                   ccAGGGTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGc                                  <  1:1486887/66‑1 (MQ=255)

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

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

                                       ggTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTcaca                                            <  1:12234822/52‑1 (MQ=255)

                                       ggTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:12860159/1‑68 (MQ=255)

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

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

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

                                        gTATGGATGTTGAACAACATGGTACGAGTGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGa                             <  1:28199653/66‑1 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:19826886/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:13230114/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:15268720/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:27481759/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:26682227/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:12484575/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:24946019/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:17194160/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:6486480/1‑67 (MQ=255)

                                        gTATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGAAATCGGTCACATGGACACGGCCCTGAg                            >  1:15933070/1‑67 (MQ=255)

                                        gTATAGATGTTGAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAg                            >  1:3999164/1‑67 (MQ=255)

                                         tATGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTc                          >  1:7680475/1‑68 (MQ=255)

                                          aTGGATGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTc                          <  1:28789687/67‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                              aTGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAttt                      >  1:12402814/1‑67 (MQ=255)

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

                                               tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                     <  1:2623450/67‑1 (MQ=255)

                                               tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                     >  1:1116692/1‑67 (MQ=255)

                                               tGTTGAACAACATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCAtttt                     >  1:638389/1‑67 (MQ=255)

                                               ggTTGAACACCATGGTACGAGCGGCGTTGGTGTAATCGCTCACAAGGACACGGCCCTGAGTCAtttc                     <  1:4943059/66‑2 (MQ=255)

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

                                                   gAACAACATGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGa                >  1:24768892/1‑68 (MQ=255)

                                                    aacaacATGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGAt               >  1:21228467/1‑68 (MQ=255)

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

                                                     acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa              >  1:20814014/1‑68 (MQ=255)

                                                     acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa              >  1:16550127/1‑68 (MQ=255)

                                                     acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATa              >  1:17975617/1‑68 (MQ=255)

                                                     acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTCCCAGATa              >  1:14802056/1‑68 (MQ=255)

                                                     acaacaTGGTACGAGAGGCGTTGGTGTAATCGGTCACAAGGACACGGCCCTGAGTCATTTTCCAGATa              >  1:3642222/1‑68 (MQ=255)

                                                        agaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGc           <  1:25809578/66‑1 (MQ=255)

                                                        acaTGGTACTAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAg            >  1:23663751/1‑67 (MQ=255)

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

                                                        acaTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAg            >  1:671744/1‑67 (MQ=255)

                                                          aTGGTACGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGcc          >  1:10297527/1‑67 (MQ=255)

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

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

                                                            ggTACTAGAGGCGTAGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACg       <  1:28941638/68‑1 (MQ=255)

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

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

                                                               aCGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTa     <  1:8694810/67‑1 (MQ=255)

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

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

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

                                                                cGAGAGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATc   >  1:21678216/1‑68 (MQ=255)

                                                                 gagaGGCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCa  <  1:19668488/68‑1 (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: `‑`