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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,106,4560CT44.8% 143.8 / inf 484G232D (GGC→GAC) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (103/163);  new base T (93/123);  total (198/286)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 512 total aligned reads displayed.

TGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAA  >  CP009273/4106390‑4106522
                                                                   |                                                                  
   ccAAACAGCGCGGCCTGCTGGTCCCCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                >  1:2716700/1‑68 (MQ=255)
   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgtcgcc                                                                >  1:30162020/1‑68 (MQ=255)
   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgtc                                                                >  1:18981703/1‑66 (MQ=255)
   ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgcc                                                                >  1:6376131/1‑68 (MQ=255)
     aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCaa                                                              >  1:16515414/1‑68 (MQ=255)
     aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCaa                                                              >  1:6113285/1‑68 (MQ=255)
      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                             <  1:15792710/68‑1 (MQ=255)
      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                             <  1:3557974/68‑1 (MQ=255)
      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAAt                                                             <  1:3326954/68‑1 (MQ=255)
      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                             <  1:21054278/68‑1 (MQ=255)
      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                             <  1:22862986/68‑1 (MQ=255)
      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                             <  1:6850595/68‑1 (MQ=255)
      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                             <  1:15549356/68‑1 (MQ=255)
      aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAAt                                                             <  1:4407984/68‑1 (MQ=255)
       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATg                                                            >  1:14514287/1‑68 (MQ=255)
       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATg                                                            >  1:20214592/1‑68 (MQ=255)
       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATg                                                            >  1:9351385/1‑68 (MQ=255)
       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                            >  1:9007025/1‑68 (MQ=255)
       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                            >  1:30504573/1‑68 (MQ=255)
       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                            >  1:27772982/1‑68 (MQ=255)
       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                            >  1:26813339/1‑68 (MQ=255)
       aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATg                                                            >  1:19382972/1‑68 (MQ=255)
        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGt                                                           >  1:29557898/1‑68 (MQ=255)
        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGt                                                           >  1:20322037/1‑68 (MQ=255)
        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGt                                                           >  1:5375288/1‑68 (MQ=255)
        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGt                                                           >  1:7105337/1‑68 (MQ=255)
        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCTGCCAATGt                                                           >  1:23893607/1‑68 (MQ=255)
        cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGt                                                           >  1:16010675/1‑68 (MQ=255)
         aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGtt                                                          >  1:21811144/1‑68 (MQ=255)
         aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGtt                                                          >  1:17336905/1‑68 (MQ=255)
         aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGtt                                                          >  1:23364020/1‑68 (MQ=255)
         aGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGtt                                                          >  1:11337236/1‑68 (MQ=255)
            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTa                                                         <  1:27344700/66‑1 (MQ=255)
            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAg                                                        >  1:10873051/1‑67 (MQ=255)
            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGt                                                       <  1:30303586/68‑1 (MQ=255)
            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGt                                                       <  1:12081370/68‑1 (MQ=255)
            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAg                                                        >  1:9596060/1‑67 (MQ=255)
            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAg                                                        <  1:12827759/67‑1 (MQ=255)
            gcgGCCTGCTGGTCACCGGCGATCCCGGAGATCGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAg                                                        <  1:22028128/67‑1 (MQ=255)
               gCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTg                                                    >  1:17895272/1‑68 (MQ=255)
                 ctgcttGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:28100365/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:14634372/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:1749595/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:25781364/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:4186235/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:19778606/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:23111914/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:19999574/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGAc                                                  <  1:4693/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGAc                                                  <  1:16332569/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                  <  1:31867817/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                  <  1:243710/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                  <  1:24968511/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                  <  1:31035521/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                  <  1:13645908/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                  <  1:20553092/68‑1 (MQ=255)
                 ctgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                                                  <  1:27999435/68‑1 (MQ=255)
                  tgctgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAcc                                                 <  1:8706320/68‑1 (MQ=255)
                     tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                                               <  1:8343502/67‑1 (MQ=255)
                     tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                                               <  1:11131662/67‑1 (MQ=255)
                     tgGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                                               <  1:25941469/67‑1 (MQ=255)
                      ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGtat                                             >  1:6732430/1‑68 (MQ=255)
                      ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                             >  1:17685171/1‑68 (MQ=255)
                      ggTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                             >  1:21259383/1‑68 (MQ=255)
                        tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGtat                                             >  1:31482417/1‑66 (MQ=255)
                        tCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                                             >  1:5587102/1‑66 (MQ=255)
                          aCCGGCGATCCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACt                                          <  1:11723198/68‑1 (MQ=255)
                           ccGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGTCAATGTTAGTCTGACCGTATACTTc                                        <  1:13414161/68‑1 (MQ=255)
                             ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCgg                                      >  1:16977308/1‑68 (MQ=255)
                             ggCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACAGTATACTTCgg                                      >  1:26474708/1‑68 (MQ=255)
                               cGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGga                                     <  1:15209082/67‑1 (MQ=255)
                               cGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGga                                     <  1:3742517/67‑1 (MQ=255)
                                       ggAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac                             <  1:14237237/68‑1 (MQ=255)
                                         agaTTGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                           >  1:28423695/1‑68 (MQ=255)
                                         agaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                           >  1:24132061/1‑68 (MQ=255)
                                         agaTTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc                           >  1:19592382/1‑68 (MQ=255)
                                            ttGGAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                        <  1:13245289/68‑1 (MQ=255)
                                            ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTCTACTTCGGAAGAACGACGCACt                        <  1:21695059/68‑1 (MQ=255)
                                            ttGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                        <  1:8168639/68‑1 (MQ=255)
                                            ttGGAATACGCGTGACGCCTTTGTCGCTAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt                        <  1:14934256/68‑1 (MQ=255)
                                               gAATACTCGTGCCGCCTTTGACGCCAATGTTAGTCTGCCCGTATAATTAGGAAGAACGAAGCACt                        >  1:7902026/1‑65 (MQ=255)
                                               gAATACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACAACGCACt                        >  1:25217042/1‑65 (MQ=255)
                                                 aTACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                    <  1:3122933/67‑1 (MQ=255)
                                                 aTACGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTg                    >  1:4930721/1‑67 (MQ=255)
                                                   aCGCGTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgca                 >  1:10592119/1‑68 (MQ=255)
                                                      cgTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca              >  1:18593774/1‑68 (MQ=255)
                                                      cgTGCCGCCTTTGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca              >  1:10034089/1‑68 (MQ=255)
                                                      cgTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca              >  1:6568936/1‑68 (MQ=255)
                                                      cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca              >  1:2469341/1‑68 (MQ=255)
                                                      cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca              >  1:8871143/1‑68 (MQ=255)
                                                      cgTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGgcagca              >  1:12282979/1‑68 (MQ=255)
                                                       gTGCCGCCTTTGCTGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCAt             <  1:28664801/68‑1 (MQ=255)
                                                           cgccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcg         <  1:28816523/68‑1 (MQ=255)
                                                             ccTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTcgc        >  1:18087566/1‑67 (MQ=255)
                                                                ttGTCGCCACTGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg    <  1:3004867/68‑1 (MQ=255)
                                                                ttGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg    <  1:26263969/68‑1 (MQ=255)
                                                                ttGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg    <  1:25555324/68‑1 (MQ=255)
                                                                ttGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg    <  1:22443151/68‑1 (MQ=255)
                                                                ttGTCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg    <  1:11932071/68‑1 (MQ=255)
                                                                ttGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCgg    <  1:7355652/68‑1 (MQ=255)
                                                                   |                                                                  
TGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTCTGGCAGCATCTCGCGCGGAA  >  CP009273/4106390‑4106522

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: