Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,106,335 C→T 8.1% M272I (ATG→ATA)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,106,3350CT 8.1% 2256.6 / 140.6 718M272I (ATG→ATA) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (225/433);  new base T (20/38);  total (246/472)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 772 total aligned reads displayed.

GCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACA  >  CP009273/4106269‑4106399
                                                                  |                                                                
gCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGc                                                                  <  1:22796500/67‑1 (MQ=255)
gCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGc                                                                  <  1:36034575/67‑1 (MQ=255)
gCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGc                                                                  <  1:20577045/67‑1 (MQ=255)
gCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGc                                                                  <  1:30992129/67‑1 (MQ=255)
gCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGc                                                                  <  1:48822221/67‑1 (MQ=255)
gCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGc                                                                  <  1:5361285/67‑1 (MQ=255)
gCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGc                                                                  <  1:1312192/67‑1 (MQ=255)
   gCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATaa                                                              >  1:7256629/1‑68 (MQ=255)
      ggCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagc                                                           >  1:37247604/1‑68 (MQ=255)
      ggCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCGTTCTCGCCAGTGTTCATCAGCATAAagc                                                           >  1:49913304/1‑68 (MQ=255)
        cGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagca                                                          >  1:43459226/1‑67 (MQ=255)
        cGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagca                                                          >  1:8840439/1‑67 (MQ=255)
         gATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagcagc                                                        <  1:1870140/68‑1 (MQ=255)
         gATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagcag                                                         <  1:28911347/67‑1 (MQ=255)
         gATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagcag                                                         <  1:44707163/67‑1 (MQ=255)
          aTGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGcc                                                       >  1:23626911/1‑68 (MQ=255)
           tggtggTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGcc                                                       <  1:7567019/67‑1 (MQ=255)
           tggtggTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCa                                                      <  1:31640489/68‑1 (MQ=255)
            ggtggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAg                                                     <  1:38120817/68‑1 (MQ=255)
            ggtggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAg                                                     <  1:37729522/68‑1 (MQ=255)
             gtggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCTGCCAGt                                                    >  1:23886091/1‑68 (MQ=255)
              tggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGt                                                    >  1:10240898/1‑67 (MQ=255)
              tggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGt                                                    <  1:44218499/67‑1 (MQ=255)
              tggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGt                                                    <  1:26522299/67‑1 (MQ=255)
              tggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGt                                                    <  1:7864927/67‑1 (MQ=255)
               ggtCTGCAGGCAGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:4367215/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:42085062/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:946138/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:3292609/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:45673747/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:47869175/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:48510906/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:49634212/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:27720926/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:26117591/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:27220825/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:11522454/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTg                                                   <  1:13851054/67‑1 (MQ=255)
               ggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGc                                                  <  1:50824912/68‑1 (MQ=255)
                 tCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCa                                                >  1:25720913/1‑68 (MQ=255)
                 tCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGcc                                                 <  1:23073815/67‑1 (MQ=255)
                 tCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGcc                                                 <  1:22043447/67‑1 (MQ=255)
                 tCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGcc                                                 <  1:49141492/67‑1 (MQ=255)
                  cagcagGCCGTGTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCa                                                <  1:17616073/67‑1 (MQ=255)
                    gcagGCCGTTTTCTGATTTCACCGGTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATa                                              <  1:26893287/67‑1 (MQ=255)
                    gcagGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAg                                             >  1:40748756/1‑68 (MQ=255)
                      agGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGGAATAGgt                                           >  1:34689270/1‑68 (MQ=255)
                      agGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAg                                             <  1:27969042/66‑1 (MQ=255)
                      agGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGgt                                           >  1:3262937/1‑68 (MQ=255)
                      agGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGgt                                           >  1:33573647/1‑68 (MQ=255)
                         ccGTTTTGTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATTAAGCAGCCAGTGCCATAGgt                                           <  1:43115428/65‑1 (MQ=255)
                         ccGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGgtgt                                         <  1:45218884/67‑1 (MQ=255)
                         ccGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGgt                                           <  1:690501/65‑1 (MQ=255)
                           gTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGttct                                      >  1:18022007/1‑68 (MQ=255)
                             tttCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGttcttc                                    <  1:46053331/68‑1 (MQ=255)
                                    tttCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATcc                             <  1:47277279/68‑1 (MQ=255)
                                    tttCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATcc                             <  1:25091905/68‑1 (MQ=255)
                                    tttCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATcc                             <  1:30134956/68‑1 (MQ=255)
                                    attCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATcc                             <  1:6416243/67‑1 (MQ=255)
                                      tCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCct                           <  1:33833133/68‑1 (MQ=255)
                                           gCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCcttctt                       >  1:36477121/1‑67 (MQ=255)
                                              ttCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCAc                   >  1:6070298/1‑68 (MQ=255)
                                              ttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGc                                                  <  1:15168091/37‑1 (MQ=255)
                                              ttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCAc                   <  1:46982158/68‑1 (MQ=255)
                                              ttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCAc                   <  1:28569386/68‑1 (MQ=255)
                                              ttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCAc                   <  1:50492580/68‑1 (MQ=255)
                                              ttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCAc                   <  1:22552531/68‑1 (MQ=255)
                                              ttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCAc                   <  1:35484664/68‑1 (MQ=255)
                                                ctcGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGCGTTCTTCGCCATCCCTTCTTTCACGc                 >  1:41540586/1‑68 (MQ=255)
                                                ctcGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGc                 >  1:41228295/1‑68 (MQ=255)
                                                   gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              >  1:21588998/1‑68 (MQ=255)
                                                   gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              <  1:8319958/68‑1 (MQ=255)
                                                   gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              <  1:1216504/68‑1 (MQ=255)
                                                   gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              <  1:14815918/68‑1 (MQ=255)
                                                   gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              >  1:29652422/1‑68 (MQ=255)
                                                   gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              <  1:3127660/68‑1 (MQ=255)
                                                   gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              >  1:2414884/1‑68 (MQ=255)
                                                   gCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              >  1:39443015/1‑68 (MQ=255)
                                                   cccAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              <  1:1001103/67‑1 (MQ=255)
                                                    ccAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              >  1:19991940/1‑67 (MQ=255)
                                                    ccAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGcaca              >  1:46470470/1‑67 (MQ=255)
                                                    ccAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTAGCCATCCCTTCTTTCACGCACaa             >  1:39888594/1‑68 (MQ=255)
                                                     cAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAAc            <  1:2122304/68‑1 (MQ=255)
                                                     cAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAAc            <  1:40265861/68‑1 (MQ=255)
                                                        tgtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTg          >  1:48189628/1‑67 (MQ=255)
                                                         gtTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGa         <  1:19554974/67‑1 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTg          <  1:15930527/66‑1 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTg          <  1:16513069/66‑1 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTg          <  1:34946652/66‑1 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTg          <  1:17252773/66‑1 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTg          <  1:42575943/66‑1 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTg          <  1:10817482/66‑1 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGa         >  1:32117128/1‑67 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGAc        <  1:39107708/68‑1 (MQ=255)
                                                         gtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGAc        <  1:24635290/68‑1 (MQ=255)
                                                           tcatcaGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCa      >  1:30717887/1‑68 (MQ=255)
                                                              tcaGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAAc   >  1:29384974/1‑68 (MQ=255)
                                                              tcaGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAAc   >  1:9171464/1‑68 (MQ=255)
                                                                aGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACa  >  1:38478787/1‑67 (MQ=255)
                                                                aGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACa  >  1:18912968/1‑67 (MQ=255)
                                                                  |                                                                
GCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACA  >  CP009273/4106269‑4106399

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