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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,106,4380AG98.4% 2879.9 / 17.6 788I238T (ATT→ACT) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base A (8/5);  new base G (344/429);  total (353/435)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 788 total aligned reads displayed.

CTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTC  >  CP009273/4106373‑4106502
                                                                 |                                                                  
 ttcttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑c                                                               >  1:52450368/1‑68 (MQ=255)
   cttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑cgc                                                             >  1:31257172/1‑68 (MQ=255)
   cttTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑cgc                                                             >  1:46774332/1‑68 (MQ=255)
     ttCACGCACAACTGACCAAACCGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGt                                                           >  1:13631514/1‑68 (MQ=255)
     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGt                                                           >  1:23446337/1‑68 (MQ=255)
     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGt                                                           >  1:10078822/1‑68 (MQ=255)
     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGt                                                           >  1:31864352/1‑68 (MQ=255)
     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGt                                                           >  1:52780593/1‑68 (MQ=255)
     ttCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGt                                                           >  1:51570447/1‑68 (MQ=255)
       cACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTg                                                          >  1:43611487/1‑67 (MQ=255)
             caACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCtt                                                   <  1:12184786/68‑1 (MQ=255)
               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCttt                                                  >  1:26790592/1‑67 (MQ=255)
               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCttt                                                  >  1:46377040/1‑67 (MQ=255)
               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCttt                                                  >  1:44212433/1‑67 (MQ=255)
               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCttt                                                  >  1:18183083/1‑67 (MQ=255)
               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCttt                                                  >  1:16666575/1‑67 (MQ=255)
               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCttt                                                  >  1:40533710/1‑67 (MQ=255)
               aCTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATA‑‑CGCGTGCCGCCttt                                                  >  1:22411781/1‑67 (MQ=255)
                 tGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTgcc                                               <  1:2021811/68‑1 (MQ=255)
                    ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTgccgcc                                            >  1:17816628/1‑68 (MQ=255)
                    ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTgccgcc                                            >  1:19440142/1‑68 (MQ=255)
                    ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTgccgcc                                            >  1:14931871/1‑68 (MQ=255)
                    ccAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTgccgc                                             <  1:15524782/67‑1 (MQ=255)
                    ccAAACAGCGCGGCCTGCTCGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCCTGCCGCCTTTgccgcc                                            >  1:33622012/1‑68 (MQ=255)
                      aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCaa                                          >  1:9641959/1‑68 (MQ=255)
                      aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCaa                                          >  1:30636742/1‑68 (MQ=255)
                      aaaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCaa                                          >  1:54727225/1‑68 (MQ=255)
                       accAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         <  1:14622710/66‑1 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGGGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         <  1:38017401/68‑1 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         >  1:35178706/1‑68 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         <  1:3307215/68‑1 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         <  1:38387237/68‑1 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         <  1:41995144/68‑1 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         <  1:42608196/68‑1 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         >  1:17400866/1‑68 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         <  1:51027554/68‑1 (MQ=255)
                       aaCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAt                                         >  1:25758018/1‑68 (MQ=255)
                        aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATg                                        >  1:32779915/1‑68 (MQ=255)
                        aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATg                                        >  1:44770889/1‑68 (MQ=255)
                        aCAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATg                                        >  1:21863427/1‑68 (MQ=255)
                         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGt                                       >  1:48434459/1‑68 (MQ=255)
                         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGt                                       >  1:55368791/1‑68 (MQ=255)
                         cAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCAGCCAATGt                                       >  1:39297785/1‑68 (MQ=255)
                             gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTa                                     <  1:36259587/66‑1 (MQ=255)
                             gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGt                                   <  1:57802646/68‑1 (MQ=255)
                             gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGt                                   <  1:732985/68‑1 (MQ=255)
                             gcgGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGt                                   <  1:35809512/68‑1 (MQ=255)
                                  ctggtgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:36804255/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCTATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCCTTGCCGCCAATGTTAGTCTGAc                              <  1:53686234/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCGGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:54332871/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGa                               <  1:8765286/67‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:18428608/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:47923189/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:49323101/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:49741548/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:45040126/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:24772322/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:44519910/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:26514268/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:50832552/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:20601437/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:51901103/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:5332998/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:27843736/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:1419607/68‑1 (MQ=255)
                                  ctgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAc                              <  1:38760693/68‑1 (MQ=255)
                                   tgctgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAcc                             <  1:25114146/68‑1 (MQ=255)
                                      tgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                           <  1:6483201/67‑1 (MQ=255)
                                      tgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                           <  1:24098934/67‑1 (MQ=255)
                                      tgGTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGt                           <  1:34441503/67‑1 (MQ=255)
                                       ggTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                         >  1:34217770/1‑68 (MQ=255)
                                       ggTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                         >  1:57378548/1‑68 (MQ=255)
                                       ggTCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                         >  1:10971537/1‑68 (MQ=255)
                                         tCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGtat                         >  1:56638068/1‑66 (MQ=255)
                                         tCACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAAGTTAGTCTGACCGtata                        <  1:50246696/67‑1 (MQ=255)
                                          cACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATTTTAGTCTGACCGTATACt                      <  1:11493458/68‑1 (MQ=255)
                                          cACCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                      <  1:37396218/68‑1 (MQ=255)
                                           aCCGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                      >  1:27532301/1‑67 (MQ=255)
                                            ccGGCGATCCCGGAGATTGGAGTACGCGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACt                      <  1:16123452/68‑1 (MQ=255)
                                             cGGCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCg                   <  1:27101843/68‑1 (MQ=255)
                                              ggCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGCCCGTTTACTTCgg                  >  1:55994432/1‑68 (MQ=255)
                                              ggCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                  >  1:22891410/1‑68 (MQ=255)
                                              ggCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                  >  1:32508702/1‑68 (MQ=255)
                                              ggCGATCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCgg                  >  1:12758207/1‑68 (MQ=255)
                                                   tCCCGGAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGgaaga              <  1:30049524/67‑1 (MQ=255)
                                                       ggAGATTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac         <  1:29195949/68‑1 (MQ=255)
                                                       cgagaATGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgac         <  1:28370695/67‑1 (MQ=255)
                                                         agaTTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAacgacg        <  1:6007539/67‑1 (MQ=255)
                                                         agaTTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGc       >  1:39644182/1‑68 (MQ=255)
                                                         agaTTGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAAGGTTAGTCTGACCGTATACTTCGGAAGAACGACGc       >  1:10655628/1‑68 (MQ=255)
                                                            ttGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc     <  1:41563296/67‑1 (MQ=255)
                                                            ttGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc     <  1:18889970/67‑1 (MQ=255)
                                                            ttGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCAc     <  1:42958320/67‑1 (MQ=255)
                                                            ttGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt    <  1:21191703/68‑1 (MQ=255)
                                                            ttGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt    <  1:45690231/68‑1 (MQ=255)
                                                            ttGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt    <  1:9229320/68‑1 (MQ=255)
                                                             tGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt   <  1:41085609/68‑1 (MQ=255)
                                                             tGGAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACtt   <  1:8090707/68‑1 (MQ=255)
                                                              ggAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACt    >  1:697560/1‑66 (MQ=255)
                                                              ggAGTA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTc  <  1:39974273/68‑1 (MQ=255)
                                                                 |                                                                  
CTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATA‑‑CGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTTCGGAAGAACGACGCACTTC  >  CP009273/4106373‑4106502

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: