Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,756,529 A→T T278S (ACC→TCC)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,756,5290AT98.8% 283.3 / NA 80T278S (ACC→TCC) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (43/36);  minor base G (1/0);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGT  >  NC_000913/1756432‑1756623
                                                                                                  |                                                                                              
cGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTcc                                                                                              <  2:2275241/100‑1 (MQ=255)
    cTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGa                                                                                          <  2:2056295/100‑1 (MQ=255)
    cTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGa                                                                                          <  2:451376/100‑1 (MQ=255)
        gtgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCt                                                                                      >  1:104376/1‑100 (MQ=255)
        gtgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCt                                                                                      >  2:2031873/1‑100 (MQ=255)
         tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTg                                                                                     >  2:244574/1‑100 (MQ=255)
         tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTg                                                                                     >  2:1179647/1‑100 (MQ=255)
         tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTg                                                                                     >  2:704325/1‑100 (MQ=255)
         tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTg                                                                                     <  2:1239861/100‑1 (MQ=255)
         tgtAGAAATCCCGGGAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTAGCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTg                                                                                     <  1:2347593/100‑1 (MQ=255)
            aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAt                                                                                  >  2:361246/1‑100 (MQ=255)
            aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAt                                                                                  >  2:983374/1‑100 (MQ=255)
            aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAt                                                                                  >  2:1488431/1‑100 (MQ=255)
             gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAtt                                                                                 >  2:312528/1‑100 (MQ=255)
             gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAtt                                                                                 >  1:1908993/1‑100 (MQ=255)
             gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAtt                                                                                 >  2:167036/1‑100 (MQ=255)
              aaaTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTc                                                                                >  2:2210541/1‑100 (MQ=255)
                 tCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCAt                                                                             <  1:792752/100‑1 (MQ=255)
                  cccGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATg                                                                            >  1:303477/1‑100 (MQ=255)
                    cGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGAt                                                                          <  2:1495707/100‑1 (MQ=255)
                    cGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGAt                                                                          <  1:2350007/100‑1 (MQ=255)
                     ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATc                                                                         >  2:147228/1‑100 (MQ=255)
                     ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATc                                                                         >  1:40908/1‑100 (MQ=255)
                      gTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCa                                                                        <  2:862311/100‑1 (MQ=255)
                        agaagaAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCaaa                                                                      <  1:1050362/100‑1 (MQ=255)
                           agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAc                                                                   >  1:2262323/1‑100 (MQ=255)
                           agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAc                                                                   <  2:40908/100‑1 (MQ=255)
                            gaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTTTCCGGTCCCGGAAAAGTCTTAACCCTTCGGCCCCGGTGGTGGGTTCCATGATCAAAAAcc                                                                  >  1:1681153/1‑100 (MQ=255)
                               gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCAc                                                               >  2:879048/1‑100 (MQ=255)
                               gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCAc                                                               >  2:2223558/1‑100 (MQ=255)
                                 tATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGc                                                             <  1:869156/100‑1 (MQ=255)
                                 tATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGc                                                             >  1:870644/1‑100 (MQ=255)
                                    cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                          <  1:879047/100‑1 (MQ=255)
                                    cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                          >  2:2017001/1‑100 (MQ=255)
                                    cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                          >  2:1386723/1‑100 (MQ=255)
                                      tCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGAc                                                        <  1:324993/100‑1 (MQ=255)
                                      tCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGAc                                                        >  2:1326621/1‑100 (MQ=255)
                                       cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt                                                       >  1:1746391/1‑100 (MQ=255)
                                       cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt                                                       >  2:1435306/1‑100 (MQ=255)
                                        gCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTc                                                      >  2:757624/1‑100 (MQ=255)
                                         cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcg                                                     <  2:285304/100‑1 (MQ=255)
                                         cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcg                                                     >  2:734236/1‑100 (MQ=255)
                                         cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGAAAAGTCGTTATCACTGCGTCCCATATGCTAGATTCCATGATCAAAAACCCACGCCCGACTcg                                                     >  1:872623/1‑100 (MQ=255)
                                          ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTCAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGAc                                                        <  2:600566/96‑1 (MQ=255)
                                             gaagaTGATGATCGAGAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCa                                                 <  1:1189774/100‑1 (MQ=255)
                                             gaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCa                                                 >  1:395778/1‑100 (MQ=255)
                                                 atgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAg                                             <  2:1746392/100‑1 (MQ=255)
                                                  tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGc                                            >  1:832662/1‑100 (MQ=255)
                                                  tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGc                                            >  2:1170249/1‑100 (MQ=255)
                                                  tgatgatCGAAAAATGTATCCGTGAACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCAACGCCCGACTCGCGCAGAAGc                                            >  2:208008/1‑100 (MQ=255)
                                                   gatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGcc                                           <  1:2210543/100‑1 (MQ=255)
                                                   gatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGcc                                           <  2:1656067/100‑1 (MQ=255)
                                                     tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                         <  2:122890/100‑1 (MQ=255)
                                                      gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                        >  2:1250694/1‑100 (MQ=255)
                                                       atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                       >  1:5969/1‑100 (MQ=255)
                                                        tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                      <  1:70182/100‑1 (MQ=255)
                                                         cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                     <  1:147229/100‑1 (MQ=255)
                                                         cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                     >  1:987631/1‑100 (MQ=255)
                                                            aaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGtt                                  <  2:402308/100‑1 (MQ=255)
                                                             aaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                 <  2:395778/100‑1 (MQ=255)
                                                              aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                >  2:688809/1‑100 (MQ=255)
                                                              aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCACGACTCGCGCAGAAGCCGGTGACGTTGc                                >  2:1077067/1‑100 (MQ=255)
                                                               aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCa                               <  2:1790453/100‑1 (MQ=255)
                                                                 gTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCaaa                             >  1:592787/1‑100 (MQ=255)
                                                                 gTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCaaa                             <  1:2008524/100‑1 (MQ=255)
                                                                    tCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGc                          <  1:312528/100‑1 (MQ=255)
                                                                    tCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGc                          <  1:1326620/100‑1 (MQ=255)
                                                                      cGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCa                        >  2:2146418/1‑100 (MQ=255)
                                                                          gCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                     <  1:361246/100‑1 (MQ=255)
                                                                          gCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                     <  1:2017003/100‑1 (MQ=255)
                                                                           cACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCt                    >  1:1659734/1‑100 (MQ=255)
                                                                             cGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                  >  1:1607228/1‑100 (MQ=255)
                                                                                 aaGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACgg              <  1:1353973/100‑1 (MQ=255)
                                                                                   gTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTa            >  2:2189847/1‑100 (MQ=255)
                                                                                   gTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTa            >  2:666465/1‑100 (MQ=255)
                                                                                      gTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg         <  2:2177071/100‑1 (MQ=255)
                                                                                       ttATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGa        <  1:2189849/100‑1 (MQ=255)
                                                                                        tATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc       >  1:1794250/1‑100 (MQ=255)
                                                                                        tATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc       <  2:832663/100‑1 (MQ=255)
                                                                                             cTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt  <  2:1607229/100‑1 (MQ=255)
                                                                                                  |                                                                                              
CGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGT  >  NC_000913/1756432‑1756623

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: