Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,757,006 G→A A437T (GCT→ACT)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,757,0060GA97.1% 205.0 / ‑2.0 69A437T (GCT→ACT) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (32/35);  total (33/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 = 7.79e-01

CGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756872‑1757130
                                                                                                                                      |                                                                                                                            
cGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGActct                                                                                                                          >  2:393226/1‑139 (MQ=255)
      ccaccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                    <  1:428003/139‑1 (MQ=255)
      ccaccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                    <  1:983406/139‑1 (MQ=255)
          caTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGt                                                                                                                <  2:1450345/139‑1 (MQ=255)
             ccTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             >  1:95766/1‑139 (MQ=255)
                ggCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGACATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGca                                                                                                          <  1:917806/139‑1 (MQ=255)
                     tGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACaa                                                                                                     >  1:729607/1‑139 (MQ=255)
                       accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             <  2:133304/129‑1 (MQ=255)
                       accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             >  1:133304/1‑129 (MQ=255)
                       accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  1:275862/1‑139 (MQ=255)
                           ccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGa                                                                                               >  1:865409/1‑139 (MQ=255)
                            caACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGActct                                                                                                                          <  2:804108/111‑1 (MQ=255)
                            caACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGActct                                                                                                                          >  1:804108/1‑111 (MQ=255)
                                gAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTg                                                                                          <  1:234686/139‑1 (MQ=255)
                                          cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcacac                                                                                                       <  1:1343807/116‑1 (MQ=255)
                                          cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcacac                                                                                                       >  2:1343807/1‑116 (MQ=255)
                                             atcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcacaca                                                                                                      >  2:172546/1‑114 (MQ=255)
                                             atcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcacaca                                                                                                      <  1:172546/114‑1 (MQ=255)
                                                aGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          <  2:729607/139‑1 (MQ=255)
                                                  ttGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTg                                                                        <  2:865409/139‑1 (MQ=255)
                                                          gAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACaa                                                                                                     >  1:811543/1‑102 (MQ=255)
                                                          gAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACaa                                                                                                     <  2:811543/102‑1 (MQ=255)
                                                               aaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAGAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  1:1150236/139‑1 (MQ=255)
                                                                     ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCac                                                     >  1:1029896/1‑139 (MQ=255)
                                                                             cAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGa                                                                                               >  1:332433/1‑89 (MQ=255)
                                                                             cAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGa                                                                                               <  2:332433/89‑1 (MQ=255)
                                                                                  tGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                    <  2:843448/63‑1 (MQ=255)
                                                                                  tGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                    >  1:843448/1‑63 (MQ=255)
                                                                                  tGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCa                                        >  2:1379656/1‑139 (MQ=255)
                                                                                          agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          <  1:1314875/97‑1 (MQ=255)
                                                                                          agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          >  2:1314875/1‑97 (MQ=255)
                                                                                            aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAacac                                            <  2:1368800/125‑1 (MQ=255)
                                                                                            aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAacac                                            >  1:1368800/1‑125 (MQ=255)
                                                                                            aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       >  2:225262/1‑130 (MQ=255)
                                                                                            aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       <  1:225262/130‑1 (MQ=255)
                                                                                            aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              >  2:1319043/1‑139 (MQ=255)
                                                                                            aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              >  1:861485/1‑139 (MQ=255)
                                                                                            aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTCTCTGGTGCACTGGTACCGAGCGGCACTACTAACTCCGCATCTGTTCACg                              <  2:1155702/139‑1 (MQ=255)
                                                                                              cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                  <  2:1151453/69‑1 (MQ=255)
                                                                                              cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                  >  1:1151453/1‑69 (MQ=255)
                                                                                              cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  2:1095978/108‑1 (MQ=255)
                                                                                              cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           >  1:1095978/1‑108 (MQ=255)
                                                                                                      cTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCa                                <  1:663220/127‑1 (MQ=255)
                                                                                                      cTGATGATTTCTACCGTCTGGGCAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCa                                >  2:663220/1‑127 (MQ=255)
                                                                                                        gatgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTgg                                                                                                            <  1:1526431/49‑1 (MQ=255)
                                                                                                        gatgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTgg                                                                                                            >  2:1526431/1‑49 (MQ=255)
                                                                                                         atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       <  1:1408481/117‑1 (MQ=255)
                                                                                                         atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       >  2:1408481/1‑117 (MQ=255)
                                                                                                           gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAacac                                            >  1:198587/1‑110 (MQ=255)
                                                                                                           gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAacac                                            <  2:198587/110‑1 (MQ=255)
                                                                                                              ttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                <  2:1014803/71‑1 (MQ=255)
                                                                                                              ttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                >  1:1014803/1‑71 (MQ=255)
                                                                                                               tCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTtg           <  1:393226/139‑1 (MQ=255)
                                                                                                                   ccGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaa       <  1:1319043/139‑1 (MQ=255)
                                                                                                                     gTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcacac                                                                                                       <  2:1610176/41‑1 (MQ=255)
                                                                                                                     gTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGcacac                                                                                                       >  1:1610176/1‑41 (MQ=255)
                                                                                                                      tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCactact                                                 >  2:981669/1‑94 (MQ=255)
                                                                                                                      tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCactact                                                 <  1:981669/94‑1 (MQ=255)
                                                                                                                      tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              <  2:1409109/113‑1 (MQ=255)
                                                                                                                      tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              >  1:1409109/1‑113 (MQ=255)
                                                                                                                       cTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaa   <  2:95766/139‑1 (MQ=255)
                                                                                                                        tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       <  2:1451952/102‑1 (MQ=255)
                                                                                                                        tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       >  1:1451952/1‑102 (MQ=255)
                                                                                                                        tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaat  <  2:1719837/139‑1 (MQ=255)
                                                                                                                        tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaat  >  1:1719837/1‑139 (MQ=255)
                                                                                                                          ggTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                >  1:1222761/1‑59 (MQ=255)
                                                                                                                          ggTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGt                                                                                <  2:1222761/59‑1 (MQ=255)
                                                                                                                          ggTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                      <  2:249626/101‑1 (MQ=255)
                                                                                                                          ggTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                      >  1:249626/1‑101 (MQ=255)
                                                                                                                                      |                                                                                                                            
CGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756872‑1757130

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: