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,0060GA100.0% 178.3 / NA 56A437T (GCT→ACT) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/28);  total (28/28)

CCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756871‑1757130
                                                                                                                                       |                                                                                                                            
ccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGActc                                                                                                                           >  1:721846/1‑139 (MQ=255)
       ccaccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                    <  2:725498/139‑1 (MQ=255)
              ccTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             >  2:984311/1‑139 (MQ=255)
                        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  1:1276715/1‑139 (MQ=255)
                        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  2:529834/1‑139 (MQ=255)
                        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  1:134683/1‑139 (MQ=255)
                        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  1:143317/1‑139 (MQ=255)
                        accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  2:565040/1‑139 (MQ=255)
                         ccaccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                  <  2:1276715/139‑1 (MQ=255)
                               aCGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                            <  1:764044/139‑1 (MQ=255)
                                cGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGtt                                                                                           <  1:734910/139‑1 (MQ=255)
                                      aaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGtt                                                                                     <  2:134683/139‑1 (MQ=255)
                                      aaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGtt                                                                                     <  1:984311/139‑1 (MQ=255)
                                                   ttGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTg                                                                        <  2:143317/139‑1 (MQ=255)
                                                            aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                      >  2:1298585/1‑132 (MQ=255)
                                                            aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                      <  1:1298585/132‑1 (MQ=255)
                                                            aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAc                                                               >  2:638472/1‑139 (MQ=255)
                                                             gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                              >  2:1300329/1‑139 (MQ=255)
                                                                      ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                            <  2:471382/100‑1 (MQ=255)
                                                                      ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                            >  1:471382/1‑100 (MQ=255)
                                                                              cAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAaca                                             <  1:853641/139‑1 (MQ=255)
                                                                                 cTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCg                                          <  1:529834/139‑1 (MQ=255)
                                                                                     ttAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                      <  1:1300329/139‑1 (MQ=255)
                                                                                          gagaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                 <  2:721846/139‑1 (MQ=255)
                                                                                           agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAacac                                            >  2:656648/1‑127 (MQ=255)
                                                                                           agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAacac                                            <  1:656648/127‑1 (MQ=255)
                                                                                           agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTACAGAGCGGTCTGGCACACaa                                                                                                     <  1:403823/70‑1 (MQ=255)
                                                                                           agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTACAGAGCGGTCTGGCACACaa                                                                                                     >  2:403823/1‑70 (MQ=255)
                                                                                             aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                       >  2:502624/1‑114 (MQ=255)
                                                                                             aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                       <  1:502624/114‑1 (MQ=255)
                                                                                              tCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                      >  1:322784/1‑98 (MQ=255)
                                                                                              tCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                      <  2:322784/98‑1 (MQ=255)
                                                                                               cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                            >  2:501775/1‑75 (MQ=255)
                                                                                               cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGt                                                                                            <  1:501775/75‑1 (MQ=255)
                                                                                                  ttCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  2:1164907/105‑1 (MQ=255)
                                                                                                  ttCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           >  1:1164907/1‑105 (MQ=255)
                                                                                                    cTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTa                       >  1:238058/1‑139 (MQ=255)
                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGTTGACGTTGTAGTTATGGTTTCTgg                                                                          <  2:890977/82‑1 (MQ=255)
                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          >  1:890977/1‑82 (MQ=255)
                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       <  2:407789/117‑1 (MQ=255)
                                                                                                          atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       >  1:407789/1‑117 (MQ=255)
                                                                                                           tgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  2:1309758/96‑1 (MQ=255)
                                                                                                           tgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           >  1:1309758/1‑96 (MQ=255)
                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGt                                                                                         <  2:490686/65‑1 (MQ=255)
                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGt                                                                                         >  1:490686/1‑65 (MQ=255)
                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTg                                                                        <  2:1306914/82‑1 (MQ=255)
                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTg                                                                        >  1:1306914/1‑82 (MQ=255)
                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       >  1:1165953/1‑115 (MQ=255)
                                                                                                            gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCAt                                       <  2:1165953/115‑1 (MQ=255)
                                                                                                               ttCTACCGTCTGGGTAATGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGa                                                            >  1:822094/1‑91 (MQ=255)
                                                                                                               ttCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGa                                                            <  2:822094/91‑1 (MQ=255)
                                                                                                                       tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           >  1:916136/1‑84 (MQ=255)
                                                                                                                       tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  2:916136/84‑1 (MQ=255)
                                                                                                                         tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaat  >  2:231353/1‑139 (MQ=255)
                                                                                                                                aGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGtt                                                                                           <  1:801055/43‑1 (MQ=255)
                                                                                                                                aGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGtt                                                                                           >  2:801055/1‑43 (MQ=255)
                                                                                                                                       |                                                                                                                            
CCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756871‑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: