Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,757,109 A→G *471* (TAA→TGA)  pykF → pyruvate kinase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,757,1090AG100.0% 75.4 / NA 22*471* (TAA→TGA) pykFpyruvate kinase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/15);  total (7/15)

AATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCAAAAATAAGTGCCTTCCCATCAAAAAAATA  >  NC_000913/1756863‑1757346
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
aaTACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGatat                                                                                                                                                                                                                                           <  1:80217/251‑1 (MQ=255)
            aTGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGaa                                                                                                                                                                                                                               <  1:374370/251‑1 (MQ=255)
                    aTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTg                                                                                                                                                                                                                       <  1:324012/251‑1 (MQ=255)
                                accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAgcg                                                                                                                                                                                                           <  1:35334/251‑1 (MQ=255)
                                 ccaccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAgcgc                                                                                                                                                                                                          >  1:131174/1‑251 (MQ=255)
                                   cccnccGAAAAAACGGCTCATCAGTTGGTACTGAGCTAAGGCTTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTTGGTAAAGAACTGGCTCTGCAGAGCGTTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTTGTGCACTGTTACCGAGCGGCACTACTAACTCCGCATCTCTTCACTGCCTGTGATATTGCTTTTGTGAATTATTTTGTATAGCGAATCGccc                                                                                                                                                                                                        <  2:369325/246‑1 (MQ=255)
                                                  gCTCATCAGTTGGTACTGAGCAAAGGCTTTTTGCCGGAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCGCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTTGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGGCCTGATTTGCGCttttt                                                                                                                                                                                         <  2:54647/251‑1 (MQ=255)
                                                         aGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGCTCCCTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCGCTTCAGAGCGGTCTTGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTCCTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGCAGCGCCCTGATGGGCGCTTTTtttattta                                                                                                                                                                                  <  2:137921/251‑1 (MQ=255)
                                                                       aaaGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGaa                                                                                                                                                                    <  1:42733/251‑1 (MQ=255)
                                                                                               aaaGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGgc                                                                                                                                                                                                <  2:70194/199‑1 (MQ=255)
                                                                                               aaaGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGgc                                                                                                                                                                                                >  1:70194/1‑199 (MQ=255)
                                                                                                                                   ggTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGAACCAAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCg                                                                                                         >  1:90995/1‑250 (MQ=255)
                                                                                                                                                           ggTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCaa                                                                                <  2:303998/251‑1 (MQ=255)
                                                                                                                                                                ggCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAAttttt                                                                           <  2:5467/251‑1 (MQ=255)
                                                                                                                                                                                                      cACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGa                                      >  2:109681/1‑250 (MQ=255)
                                                                                                                                                                                                             aCCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCaa                              <  2:185573/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggCACTACTAACACCGCATCTGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCAAAAATAAg                       <  2:44417/251‑1 (MQ=255)
                                                                                                                                                                                                                       actactAACACCGCATATGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCAAAAATAAGTGc                    >  1:363412/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tGTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGAAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCAAAAATAAGTGCCTTCCCATCAAAAAAat   >  1:33777/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gTTCACGTCCTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAGCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCTCTTTATCTAAGATGAATCCGATTGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCAAAAATAAGTGCCTTCTCATCAAAAAAata  <  2:189798/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCaaaaa                           <  1:56087/217‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTGTGATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCaaaaa                           >  2:56087/1‑217 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
AATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGATTTCACTATATAATCTCACTTTATCTAAGATGAATCCGATGGAAGCATCCTGTTTTCTCTCAATTTTTTTATCTAAAACCCAGCGTTCGATGCTTCTTTGAGCGAACGATCAAAAATAAGTGCCTTCCCATCAAAAAAATA  >  NC_000913/1756863‑1757346

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