BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,457	C→T	15.0%	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	15.0%	1574.0 / 249.4	609	G232S (GGC→AGC) ‡	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (207/310); new base T (29/62); total (236/373)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 669 total aligned reads displayed.

TGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAATA  >  CP009273/4106390‑4106524

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

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

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

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

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

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

   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgc                                                                   <  1:31814465/67‑1 (MQ=255)

      aaCAGCGGGGCCTGCTGGTCACCGGCGTTCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:23076205/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                               <  1:788325/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                               <  1:16531809/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                               <  1:16699926/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                               <  1:7188704/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                               <  1:6700919/68‑1 (MQ=255)

      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                               <  1:18245769/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

      aaCACCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                               <  1:7043461/68‑1 (MQ=255)

       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATg                                                              >  1:32885767/1‑68 (MQ=255)

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

        cAGCGCGGCCTGCTGGTCACCGGCGATCGCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGt                                                             >  1:3834764/1‑68 (MQ=255)

        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGt                                                             >  1:28502747/1‑68 (MQ=255)

        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAACGt                                                             >  1:31319129/1‑68 (MQ=255)

         aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGtt                                                            >  1:30456138/1‑68 (MQ=255)

         aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGtt                                                            >  1:22132464/1‑68 (MQ=255)

         aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATAGGAATACGCGTGCCGCCTTTGCTGCCAATGtt                                                            >  1:32412953/1‑68 (MQ=255)

            gcggcgTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTCGCCGCCAATGTTAGt                                                         <  1:26906578/68‑1 (MQ=255)

            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGt                                                             <  1:28942446/64‑1 (MQ=255)

            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGt                                                         <  1:10524254/68‑1 (MQ=255)

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

                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATAGGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                    <  1:7579121/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                  tgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAcc                                                   <  1:20002289/68‑1 (MQ=255)

                     tgGTCACCGGCGTTCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGt                                                 <  1:1489548/67‑1 (MQ=255)

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

                      ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGtat                                               >  1:12647998/1‑68 (MQ=255)

                         cACCGGCGATCCCGGCGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                            <  1:22732642/68‑1 (MQ=255)

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

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

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

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

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

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

                          aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACt                                            >  1:24153898/1‑67 (MQ=255)

                          aCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                            >  1:30903487/1‑67 (MQ=255)

                               cGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGGCAATGTTAGTCTGAACGTATACTTCGgaa                                      >  1:18587116/1‑68 (MQ=255)

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

                                   cccGGATATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAac                                  <  1:2234786/68‑1 (MQ=255)

                                    ccGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacg                                 <  1:1293235/68‑1 (MQ=255)

                                     cGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacga                                <  1:34414101/68‑1 (MQ=255)

                                     cGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacga                                <  1:33668210/68‑1 (MQ=255)

                                      ggAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                               <  1:8557643/68‑1 (MQ=255)

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

                                        agaTTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:4232274/1‑68 (MQ=255)

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

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

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

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

                                         gaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:11673310/1‑67 (MQ=255)

                                         gaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                             >  1:21861696/1‑67 (MQ=255)

                                           ttGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                          >  1:24439058/1‑68 (MQ=255)

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

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

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

                                           ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTACACTTCGGAAGAACGACGCAc                           <  1:27918084/67‑1 (MQ=255)

                                              gAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                          >  1:6059299/1‑65 (MQ=255)

                                                aTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGCACGACGCACTTCTg                      >  1:27243304/1‑67 (MQ=255)

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

                                                 tACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgc                    >  1:3619021/1‑68 (MQ=255)

                                                  aCGCGTGCGGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                   >  1:17880803/1‑68 (MQ=255)

                                                    gcgTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGCAAGAACGACGCACTTCTGgcagc                 >  1:23514787/1‑68 (MQ=255)

                                                    gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                 >  1:10032395/1‑68 (MQ=255)

                                                    gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                 >  1:25928399/1‑68 (MQ=255)

                                                    gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                 >  1:29611927/1‑68 (MQ=255)

                                                    gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                 >  1:29915061/1‑68 (MQ=255)

                                                    gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                 >  1:25688457/1‑68 (MQ=255)

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

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

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

                                                        gccgccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAtc              >  1:5715773/1‑67 (MQ=255)

                                                            ccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgc          >  1:2168243/1‑67 (MQ=255)

                                                              tttGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgcgcg       <  1:29304429/68‑1 (MQ=255)

                                                               ttGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg      <  1:31511126/68‑1 (MQ=255)

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

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

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

                                                                tGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGa     <  1:18717671/68‑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: `‑`