BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,106,456	C→T	12.0%	G232D (GGC→GAC) ‡	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,106,456	0	C	T	12.0%	1730.3 / 186.8	609	G232D (GGC→GAC) ‡	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (206/330); new base T (30/43); total (236/373)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 658 total aligned reads displayed.

CTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAAT  >  CP009273/4106389‑4106523

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

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

    ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgtc                                                                 >  1:23172020/1‑66 (MQ=255)

    ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgtc                                                                 >  1:27535068/1‑66 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgc                                                                  >  1:10761959/1‑62 (MQ=255)

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

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

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

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

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

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

             gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTACGCGTGCCGCCTTTGCCGCCAATGt                                                            <  1:32106083/64‑1 (MQ=255)

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

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

             gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAg                                                         <  1:8952371/67‑1 (MQ=255)

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

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

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

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

                  ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                   <  1:28244355/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                   tgctgGTCACCGGCGTTCCCGGAGATTTGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGAc                                                   <  1:2323578/67‑1 (MQ=255)

                   tgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAcc                                                  >  1:1165671/1‑68 (MQ=255)

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

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

                      tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGt                                                <  1:26301140/67‑1 (MQ=255)

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

                       ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGtat                                              >  1:3613503/1‑68 (MQ=255)

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

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

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

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

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

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

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

                              ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCgg                                       >  1:12250345/1‑68 (MQ=255)

                              ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                                       >  1:22909635/1‑68 (MQ=255)

                                cGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTCGga                                      <  1:23505973/67‑1 (MQ=255)

                                   tCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGgaaga                                   >  1:11156398/1‑67 (MQ=255)

                                     ccGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacg                                <  1:21105492/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                aaTACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                                           <  1:15892323/46‑1 (MQ=255)

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

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

                                                  tACGCGTGCCGCCTGTGCTGCCAATGTTCGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgc                   <  1:31256547/68‑1 (MQ=255)

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

                                                    cgcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                  >  1:16707761/1‑67 (MQ=255)

                                                     gcgTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                >  1:13167071/1‑68 (MQ=255)

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

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

                                                     gcgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagc                <  1:175725/68‑1 (MQ=255)

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

                                                      cgTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca               >  1:646822/1‑68 (MQ=255)

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

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

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

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

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

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

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

                                                                ttGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg     <  1:13701491/68‑1 (MQ=255)

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

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

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

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

                                                                  gccgccAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGaa   <  1:2612571/68‑1 (MQ=255)

CTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAAT  >  CP009273/4106389‑4106523

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