Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,106,955 C→T 100% A66T (GCG→ACG)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,106,9550CT96.1% 2114.6 / 41.9 583A66T (GCG→ACG) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base T (293/265);  total (302/281)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 585 total aligned reads displayed.

AGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAA  >  CP009273/4106890‑4107019
                                                                 |                                                                
 gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCa                                                                >  1:21137841/1‑67 (MQ=255)
   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGc                                                              >  1:17760297/1‑67 (MQ=255)
   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGc                                                              >  1:10541410/1‑67 (MQ=255)
   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGc                                                              >  1:20111991/1‑67 (MQ=255)
   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGc                                                              >  1:22329088/1‑67 (MQ=255)
   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGc                                                              >  1:25182640/1‑67 (MQ=255)
   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGc                                                              >  1:27178209/1‑67 (MQ=255)
   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCa                                                             <  1:5921025/68‑1 (MQ=255)
   ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                             <  1:25815430/68‑1 (MQ=255)
   ggTATCACGCTCGTTCGTAAAACCGATAGCTGAAATATCATCGGAAATGATAACCGCTTTCCTCAGc                                                              >  1:27710209/1‑67 (MQ=255)
    gTTTCACGCTGGTTCGTAATAGCGATAGGTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                            <  1:4351631/68‑1 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                            <  1:2545044/67‑1 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                            <  1:18459087/67‑1 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                            <  1:2444250/67‑1 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                            <  1:28542508/67‑1 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                            <  1:15197408/67‑1 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACt                                                           >  1:27378296/1‑68 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACt                                                           >  1:28879525/1‑68 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACt                                                           >  1:20799935/1‑68 (MQ=255)
     tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACt                                                           >  1:15337331/1‑68 (MQ=255)
      ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACtt                                                          <  1:11274/68‑1 (MQ=255)
        cACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCt                                                        <  1:23427030/68‑1 (MQ=255)
         aCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCt                                                        >  1:32153907/1‑67 (MQ=255)
          cGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTa                                                       >  1:15783285/1‑67 (MQ=255)
           gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTAcc                                                     >  1:5745416/1‑68 (MQ=255)
            cTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCa                                                    >  1:9644760/1‑68 (MQ=255)
             tGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAg                                                   <  1:3164822/68‑1 (MQ=255)
              ggTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGc                                                  <  1:3543207/68‑1 (MQ=255)
                  cGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCAc                                                            <  1:17545410/54‑1 (MQ=255)
                     aaTACCTATAGCTGCAATTTGATCGGAACTGATATCGTCTTTCGTCAGCACTTCTACCAGCGGGGAGc                                           <  1:33272887/68‑1 (MQ=255)
                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGa                                             >  1:6317868/1‑66 (MQ=255)
                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGa                                             >  1:24850589/1‑66 (MQ=255)
                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAg                                            <  1:14479728/67‑1 (MQ=255)
                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAg                                            <  1:24255660/67‑1 (MQ=255)
                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGc                                           <  1:33461754/68‑1 (MQ=255)
                     aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGc                                           <  1:30394106/68‑1 (MQ=255)
                       tACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCtt                                         >  1:5305266/1‑68 (MQ=255)
                       tACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCtt                                         >  1:3395787/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:34251917/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:17167283/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:8760466/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:30670304/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:29321511/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:12575245/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:26837135/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:20428999/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCttt                                        >  1:14942757/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTCGAGCttt                                        >  1:29984542/1‑68 (MQ=255)
                        aCCGATAGCTGCAATTTGATAGGAACTGATATCGGCTTTCGTCAGCACTTCTAAAAGCGTGGAGCttt                                        >  1:12270615/1‑68 (MQ=255)
                          cGATAGCTGCAATTTGATCGGAACTGATATCGGCATTCGTCAGCCCTTCAACCAGCGTGGAGCTTTgg                                      >  1:8984482/1‑68 (MQ=255)
                             tAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTg                                   >  1:16138766/1‑68 (MQ=255)
                               gCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                  >  1:26693457/1‑67 (MQ=255)
                               gCTGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCAACCAGCGTGGAGCTTTGGGTgg                                  >  1:5499575/1‑67 (MQ=255)
                                 tGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                 <  1:15596997/66‑1 (MQ=255)
                                 tGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                 <  1:1196687/66‑1 (MQ=255)
                                 tGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                 <  1:19593793/66‑1 (MQ=255)
                                 tGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                 <  1:6537931/66‑1 (MQ=255)
                                 tGCAATTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                 <  1:17956334/66‑1 (MQ=255)
                                    aaTTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                            <  1:27638242/68‑1 (MQ=255)
                                    aaTTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                            <  1:10162977/68‑1 (MQ=255)
                                    aaTTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                            <  1:8294561/68‑1 (MQ=255)
                                    aaTTTGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGa                            <  1:14663334/68‑1 (MQ=255)
                                      tttGATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAtt                          >  1:10741253/1‑68 (MQ=255)
                                         gATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                       <  1:12774004/68‑1 (MQ=255)
                                         gATCGGAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                       <  1:2313856/68‑1 (MQ=255)
                                         gATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGACCCAGATTTcc                       <  1:32840226/68‑1 (MQ=255)
                                             ggAACTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTg                   <  1:28740568/68‑1 (MQ=255)
                                                aCTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg                 >  1:21614510/1‑67 (MQ=255)
                                                aCTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTggg                 >  1:9168455/1‑67 (MQ=255)
                                                aCTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                <  1:13926474/68‑1 (MQ=255)
                                                aCTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                <  1:2371503/68‑1 (MQ=255)
                                                aCTGATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                <  1:34699215/68‑1 (MQ=255)
                                                aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGt                <  1:4007748/68‑1 (MQ=255)
                                                   gATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt             <  1:4496537/68‑1 (MQ=255)
                                                   gATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt             <  1:13154940/68‑1 (MQ=255)
                                                   gATATCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgt             <  1:1652916/68‑1 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           <  1:7745975/68‑1 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:16865491/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:6767267/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:7600486/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:4956812/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           <  1:14231024/68‑1 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:22779649/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:10024932/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:18712481/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:11041029/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:19031595/1‑68 (MQ=255)
                                                     tatCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:34396011/1‑68 (MQ=255)
                                                       tCGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt           >  1:31042991/1‑66 (MQ=255)
                                                        cGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt        <  1:7446976/68‑1 (MQ=255)
                                                        cGGCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt        >  1:19280961/1‑68 (MQ=255)
                                                         ggCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCt        >  1:22606100/1‑67 (MQ=255)
                                                         ggCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa       <  1:22040131/68‑1 (MQ=255)
                                                          gCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa       >  1:13419934/1‑67 (MQ=255)
                                                          gCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa       >  1:28185267/1‑67 (MQ=255)
                                                          gCTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa       >  1:31332983/1‑67 (MQ=255)
                                                           cTTTCGTCAGCACTTCTACCTGCGGGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc      <  1:30214010/67‑1 (MQ=255)
                                                           cTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc      <  1:26344677/67‑1 (MQ=255)
                                                           cTTTCGTCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc      <  1:7096974/67‑1 (MQ=255)
                                                                 ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaa  >  1:11744286/1‑65 (MQ=255)
                                                                 |                                                                
AGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAA  >  CP009273/4106890‑4107019

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: