Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,756,910 (A)6→5 coding (1213/1413 nt) pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,756,9050A.100.0% 147.2 / NA 35coding (1208/1413 nt)pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (20/15);  total (20/15)

TGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGG  >  NC_000913/1756772‑1757034
                                                                                                                                     |                                                                                                                                 
tGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGTCACCATCCTGGCACTGACCACCAACGaaaaac                                                                                                                              >  2:446453/1‑138 (MQ=255)
      ggtgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGTCACCATCCTGGCACTGACCACCAACGAAAAACGGCtca                                                                                                                       <  1:446453/139‑1 (MQ=255)
        tgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCtcatc                                                                                                                     <  2:305948/139‑1 (MQ=255)
                  aCTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACt                                                                                                           >  2:73248/1‑139 (MQ=255)
                                              tCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGtt                                                                               <  2:250553/139‑1 (MQ=255)
                                                  ggTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAag                                                                           >  1:897555/1‑139 (MQ=255)
                                                         aCTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAc                                                                    >  1:370899/1‑139 (MQ=255)
                                                             aGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                         <  1:147396/130‑1 (MQ=255)
                                                             aGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                         >  2:147396/1‑130 (MQ=255)
                                                                                cgcAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTc                                             <  1:73248/139‑1 (MQ=255)
                                                                                               cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAGAACGTCTCATCAGTTGGTACTGAGCAAAGGCg                                                                                                <  1:542950/73‑1 (MQ=255)
                                                                                               cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCg                                                                                                >  2:542950/1‑73 (MQ=255)
                                                                                                 tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGc                                                                                    >  2:873847/1‑83 (MQ=255)
                                                                                                 tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGc                                                                                    <  1:873847/83‑1 (MQ=255)
                                                                                                 tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTaaa                                                                            >  2:298268/1‑91 (MQ=255)
                                                                                                 tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTaaa                                                                            <  1:298268/91‑1 (MQ=255)
                                                                                                 tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAc                                                                    <  2:904580/99‑1 (MQ=255)
                                                                                                 tCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAc                                                                    >  1:904580/1‑99 (MQ=255)
                                                                                                  cccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCtca                                                                                                                       >  2:444993/1‑47 (MQ=39)
                                                                                                  cccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCtca                                                                                                                       <  1:444993/47‑1 (MQ=39)
                                                                                                  cccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctc                           >  2:371925/1‑139 (MQ=255)
                                                                                                                    ggCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                      >  1:391081/1‑62 (MQ=255)
                                                                                                                    ggCACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                                                                                      <  2:391081/62‑1 (MQ=255)
                                                                                                                      cACTGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcacac       <  1:378199/139‑1 (MQ=255)
                                                                                                                         tGACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACaaa    >  2:196144/1‑139 (MQ=255)
                                                                                                                          gACCACCAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAg   >  1:770268/1‑139 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                         <  2:488615/68‑1 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                         >  1:488615/1‑68 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAt                                                                       <  1:23914/70‑1 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAt                                                                       >  2:23914/1‑70 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCt                                                   >  1:53254/1‑90 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCt                                                   <  2:53254/90‑1 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg  >  2:62307/1‑139 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg  >  1:833231/1‑139 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg  >  1:351842/1‑139 (MQ=255)
                                                                                                                           accaccAACGAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACAGTCTGGGTAAAGAACTGGCTCTGAAGAGCGGTCTGGCACACAAAgg  >  1:903762/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                 
TGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGG  >  NC_000913/1756772‑1757034

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