BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,438	A→G	9.0%	I238T (ATT→ACT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,438	0	A	G	9.0%	2329.5 / 150.0	716	I238T (ATT→ACT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base A (287/359); new base G (28/36); total (317/399)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 753 total aligned reads displayed.

CCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTG  >  CP009273/4106371‑4106504

cccTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGCTTGGaa                                                                    >  1:44303958/1‑68 (MQ=255)

   ttcttTCACGCACAACTGACCAAACAGCGCGGCCCGCTGGTCACCGGCGATCCCGGAGATTGGAATAc                                                                 >  1:3224919/1‑68 (MQ=255)

      tttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTAcgcg                                                              >  1:35128656/1‑68 (MQ=255)

      tttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATAcgcg                                                              >  1:19716785/1‑68 (MQ=255)

      tttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATAcgcg                                                              <  1:6397273/68‑1 (MQ=255)

       ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGt                                                             >  1:13409295/1‑68 (MQ=255)

       ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGt                                                             >  1:39559752/1‑68 (MQ=255)

       ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGt                                                             >  1:36366834/1‑68 (MQ=255)

       ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGt                                                             >  1:24840580/1‑68 (MQ=255)

           cGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTgccg                                                         <  1:34555114/68‑1 (MQ=255)

            gCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTgccg                                                         >  1:26969970/1‑67 (MQ=255)

              acaACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCt                                                      >  1:18458702/1‑68 (MQ=255)

              acaACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATACCGGAGATTGGAATACGCGTGCAGCCt                                                      >  1:9429903/1‑68 (MQ=255)

                 aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCttt                                                    >  1:27815602/1‑67 (MQ=255)

                 aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCttt                                                    >  1:47000912/1‑67 (MQ=255)

                 aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCttt                                                    >  1:514298/1‑67 (MQ=255)

                 aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCttt                                                    >  1:13188660/1‑67 (MQ=255)

                 aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCttt                                                    >  1:29014406/1‑67 (MQ=255)

                 aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCttt                                                    >  1:26021184/1‑67 (MQ=255)

                   tGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgct                                                 <  1:10407188/68‑2 (MQ=255)

                      ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGATTACGCGTGCCGCCTTtgctgc                                               <  1:42326131/67‑1 (MQ=255)

                      ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTtgctgc                                               <  1:47786326/67‑1 (MQ=255)

                      ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgt                                               <  1:15396158/67‑2 (MQ=255)

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

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

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

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

                        aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCaa                                            >  1:21634336/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                          aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATg                                          >  1:37502236/1‑68 (MQ=255)

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

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

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

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

                            aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGtt                                        >  1:44862129/1‑68 (MQ=255)

                               gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGGGCCGCCTTTGCCGCCAATGTTAg                                      <  1:37796482/67‑1 (MQ=255)

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

                               gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAg                                      <  1:34102303/67‑1 (MQ=255)

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

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

                               gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTCAg                                      >  1:20697629/1‑67 (MQ=255)

                                    ctgctgGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                <  1:42755463/68‑1 (MQ=255)

                                    ctgctgGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                <  1:29794977/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                     tgctgGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                <  1:35635428/67‑1 (MQ=255)

                                      gctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCg                              >  1:1714113/1‑68 (MQ=255)

                                        tgGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                             <  1:43524603/67‑1 (MQ=255)

                                        tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                             <  1:2510232/67‑1 (MQ=255)

                                        tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                             <  1:36219150/67‑1 (MQ=255)

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

                                         ggTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                           >  1:4990201/1‑68 (MQ=255)

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

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

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

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

                                            cACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACt                        <  1:48488320/68‑1 (MQ=255)

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

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

                                               cGGCGATCCCGGAGATTGGATTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCg                     <  1:7682447/68‑1 (MQ=255)

                                               cGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCg                     <  1:8863690/68‑1 (MQ=255)

                                                 gCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTCTCCTTCGga                   >  1:28856643/1‑68 (MQ=255)

                                                  cGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaa                  >  1:12212908/1‑68 (MQ=255)

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

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

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

                                                           agaTTGGAATACGCGTGCCGACTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc         >  1:7879432/1‑68 (MQ=255)

                                                              ttGGAGTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt      <  1:13995791/68‑1 (MQ=255)

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

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

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

                                                                   aTACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg  >  1:3173778/1‑67 (MQ=255)

                                                                   aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGCCGCACTTCTg  >  1:17745384/1‑67 (MQ=255)

                                                                   aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg  >  1:695103/1‑67 (MQ=255)

                                                                   aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg  >  1:29627425/1‑67 (MQ=255)

CCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTG  >  CP009273/4106371‑4106504

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