Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,892,999 A→C 7.4% N297T (AAC→ACC)  pyk → pyruvate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,892,9990AC 7.4% 159.3 / 3.9 54N297T (AAC→ACC) pykpyruvate kinase
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base C (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

ACCGGCCGAGAGCGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGAACGGCCCGGACTACCAGGCCCAGCTCGACGTCG  >  NC_002947/4892853‑4893142
                                                                                                                                                  |                                                                                                                                               
aCCGGCCGAGAGCGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCACcg                                                                                                                                               >  7:66292/1‑149 (MQ=255)
   ggCCGAGAGCGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCgg                                                                                                                                            <  6:250895/149‑1 (MQ=255)
         gagCGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTg                                                                                                                                          >  3:99038/1‑145 (MQ=255)
         gagCGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTg                                                                                                                                          <  4:99038/145‑1 (MQ=255)
           gCGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGaa                                                                                                                                    <  2:81375/149‑1 (MQ=255)
           gCGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGaa                                                                                                                                    <  1:120585/149‑1 (MQ=255)
            cGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCACCGCGTTGGccacgt                                                                                                                                   <  8:66292/149‑5 (MQ=255)
              tGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGt                                                                                                                                 <  7:251408/149‑1 (MQ=255)
                  gCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCGCGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCg                                                                                                                             <  7:44836/149‑1 (MQ=255)
                  gCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGcc                                                                                                                              >  4:204945/1‑148 (MQ=255)
                      aTCCAGAAGCGCATCATCGGCACCTGGCGCCTGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGc                                                                                                                         <  2:252796/149‑1 (MQ=255)
                      aTCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATg                                                                                                                          >  5:389529/1‑148 (MQ=255)
                            aaGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGAt                                                                                                                   <  3:204945/149‑1 (MQ=255)
                                gcATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCAGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTg                                                                                                               >  8:44416/1‑149 (MQ=255)
                                  atcatcGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTc                                                                                                             >  3:81585/1‑149 (MQ=255)
                                     atcGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGc                                                                                                          <  4:81585/149‑1 (MQ=255)
                                      tcGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTTGCCGAAGGTGCCGATGCAGTGATGCTGTCGGct                                                                                                         >  5:258852/1‑148 (MQ=255)
                                        ggCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCg                                                                                                        <  6:389529/148‑1 (MQ=255)
                                                   ccAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCgg                                                                                            <  8:281336/149‑1 (MQ=255)
                                                   ccAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCgg                                                                                            <  1:281346/149‑1 (MQ=255)
                                                        tAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGAt                                                                                       >  8:246825/1‑149 (MQ=255)
                                                          ggCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGAt                                                                                       <  7:246825/147‑1 (MQ=255)
                                                           gCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTAc                                                                                    >  4:303489/1‑149 (MQ=255)
                                                               aCCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGc                                                                                <  3:303489/149‑1 (MQ=255)
                                                                 ccGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTg                                                                              >  1:222291/1‑149 (MQ=255)
                                                                  cGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGTGTGATTACCCGCTg                                                                              >  1:140440/1‑148 (MQ=255)
                                                                    tggtggtggCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGaa                                                                           <  6:279178/149‑1 (MQ=255)
                                                                     ggtggtggCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCg                                                                                                            >  2:167977/1‑115 (MQ=255)
                                                                     ggtggtggCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCg                                                                                                            <  1:167977/115‑1 (MQ=255)
                                                                       tggtggCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGcc                                                                        <  2:222291/149‑1 (MQ=255)
                                                                           ggCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGc                                                                    <  7:226814/149‑1 (MQ=255)
                                                                              cACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTg                                                                     >  8:56184/1‑145 (MQ=255)
                                                                                cGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAgatg                                                               <  5:390781/149‑1 (MQ=255)
                                                                                    gATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAgatgat                                                             >  8:22967/1‑147 (MQ=255)
                                                                                         tCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTAGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGc                                                         >  4:87309/1‑146 (MQ=255)
                                                                                                gATGAGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCTTGCAGATGATGAGCACGATCATcc                                               >  4:20312/1‑149 (MQ=255)
                                                                                                   gCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGc                                             >  6:32340/1‑148 (MQ=255)
                                                                                                         ctcGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTTGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGtt                                       <  3:87309/148‑1 (MQ=255)
                                                                                                         ctcGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGtt                                       <  7:22967/148‑1 (MQ=255)
                                                                                                         ctcGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGtt                                       <  7:295153/148‑1 (MQ=255)
                                                                                                            gCCAGCGCCAACCCGCGCCGACGTCACGGACGTGGCAAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCGGTGCAGATGATGAGCAAGATCATCCGCCAGGTTga                                     <  3:20312/147‑1 (MQ=255)
                                                                                                            gCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTgag                                    >  1:149532/1‑148 (MQ=255)
                                                                                                             ccAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCaa                                                       >  6:348655/1‑128 (MQ=255)
                                                                                                             ccAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCaa                                                       <  5:348655/128‑1 (MQ=255)
                                                                                                             ccAGCGCCAACCCGCGCAGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGaa                                  <  7:44416/149‑1 (MQ=255)
                                                                                                                    cAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAgatg                                                               <  5:178394/113‑1 (MQ=255)
                                                                                                                    cAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAgatg                                                               >  6:178394/1‑113 (MQ=255)
                                                                                                                    cAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGAACGGCCCg                           >  6:173770/1‑149 (MQ=255)
                                                                                                                        ccGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGAACGGCCCGGAc                        >  2:153519/1‑148 (MQ=255)
                                                                                                                           cgcCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCACGTTGAGAACGGCCCGGACTAcc                    <  7:56184/149‑1 (MQ=255)
                                                                                                                                         cttcggCCACCGCGTTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGAACGGCCCGGACTACCAgg                 >  2:214585/5‑138 (MQ=255)
                                                                                                                                         cttcggCCACCGCGTTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGAACGGCCCGGACTACCAgg                 <  1:214585/134‑1 (MQ=255)
                                                                                                                                         aCGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGAACGGCCCGGACTACCAGGCCCAGCTCGAc      >  5:251998/1‑149 (MQ=255)
                                                                                                                                              gCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGAACGGCCCGGACTACCAGGCCCAGCTCGACGTcg  <  5:32340/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACCGGCCGAGAGCGTGCCGCAGATCCAGAAGCGCATCATCGGCACCTGCCGCCAGCTAGGCAAACCCGTGGTGGTGGCCACGCAGATGCTCGAGTCGATGCGTTTCTCGCCAGCGCCAACCCGCGCCGAAGTCACGGACGTGGCCAACGCGGTGGCCGAAGGTGCCGATGCAGTGATGCTGTCGGCCGAAACCGCCTCGGGTGATTACCCGCTGGAAGCCGTGCAGATGATGAGCAAGATCATCCGCCAGGTTGAGAACGGCCCGGACTACCAGGCCCAGCTCGACGTCG  >  NC_002947/4892853‑4893142

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