Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,306,896 G→A P302S (CCA→TCA)  alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,8960GA100.0% 42.7 / NA 18P302S (CCA→TCA) alsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/7);  total (11/7)

CGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGG  >  NC_000913/4306665‑4307105
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                 
ngTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGACCGCGCCATTTTCTTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGcc                                                                                                                                                                                                >  1:60613/2‑251 (MQ=255)
              ngAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCAcc                                                                                                                                                                                  >  1:192285/2‑251 (MQ=255)
                      cGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCCTTTGTGCTGCAGCGGGAAAATTACAGTACAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAgcgc                                                                                                                                                                          >  1:159219/1‑251 (MQ=255)
                                     gAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATg                                                                                                                                                                                                    >  1:106334/1‑210 (MQ=255)
                                            gAGCGGAACGCGCCATTTTATTAGGCATTTGTGATGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGGTTTATGTGTTTATTATATATATTGTTCCTGTCGTTGGCACAGCGGGACATGGCGCTACTGTGCCAACTTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAATCGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAggc                                                                                                                                                    <  1:157923/251‑1 (MQ=255)
                                                                               gAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAGGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGAc                                                                                                                 <  1:305122/251‑1 (MQ=255)
                                                                                      naaaaTTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTCTTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGAGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGgcgc                                                                                                            >  1:33927/2‑249 (MQ=255)
                                                                                       aaaaTTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGGTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAgg                                                                                                         <  1:148526/251‑1 (MQ=255)
                                                                                                  naaaCTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGt                                                                                              >  1:32385/2‑251 (MQ=255)
                                                                                                   aaaCTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTCGGCCCAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTc                                                                                             >  1:280848/1‑251 (MQ=255)
                                                                                                                   aaaTAATATGAGTCGTTTTATGTGTTTATTATATATATTGTTCCTGTCGTTGGCAAAGCGGGGAATTGCGCTCCTGTGACACCGTTTTACTCATGGGGCTTTAGCACAGCAATGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGCGTCATGGCAACGagaggct                                                                             >  1:287785/1‑248 (MQ=255)
                                                                                                                                                     atnnNNNACCNGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAg                                                                                                                                                      <  2:192285/144‑1 (MQ=255)
                                                                                                                                                             ccTGTCGTTGGCACACCGGGACATTGCGCTACTGTGCCCGCGTTTTACTCATGGGGCTTTAGCAAAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGCAGCTGAGGCGGCAATTAAGCGCACGACATGATGCGGCAGTGGACGGCGCCGGTACTTTTGGGTCATGGCAACGCGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACg                                                >  1:45486/1‑238 (MQ=255)
                                                                                                                                                                                     tGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCAccc                                                                                                                                                                                 >  1:13458/1‑85 (MQ=255)
                                                                                                                                                                                     tGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCAccc                                                                                                                                                                                 <  2:13458/85‑1 (MQ=255)
                                                                                                                                                                                              gtgCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCtgg  <  1:119288/251‑1 (MQ=255)
                                                                                                                                                                                                                      tttAGCCCAGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACtttt                                                                                                  <  1:43782/131‑1 (MQ=255)
                                                                                                                                                                                                                             caGAACTGTGACAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAACGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCg                                                               >  1:147960/1‑159 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                 
CGTTATCATAATTGGGAAATTGCGACTTTAATAAGTGGAAGTGTGAGCGGAACGCGCCATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGG  >  NC_000913/4306665‑4307105

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