Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,163,578 G→A A107T (GCA→ACA)  thiK → thiamine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,163,5780GA100.0% 129.8 / NA 34A107T (GCA→ACA) thiKthiamine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/20);  total (14/20)

GCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGGCAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCTGCGCAAA  >  NC_000913/1163430‑1163721
                                                                                                                                                    |                                                                                                                                               
gcAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACa                                                                                                                                               <  2:340545/151‑1 (MQ=255)
   gTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGc                                                                                                                                            <  2:240819/151‑1 (MQ=255)
   gTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGc                                                                                                                                            >  2:644895/1‑151 (MQ=255)
              tGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCAGATACCAACGAACTGACAGGCTTGCTGtatt                                                                                                                                  >  2:239039/1‑150 (MQ=255)
                          gggCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATcaaca                                                                                                                     <  1:314466/151‑1 (MQ=255)
                               ttATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCAc                                                                                                                >  1:222442/1‑151 (MQ=255)
                                tATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCAc                                                                                                                <  2:416423/150‑1 (MQ=255)
                                tATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCAc                                                                                                                <  1:632915/150‑1 (MQ=255)
                                tATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCAc                                                                                                                <  1:21380/150‑1 (MQ=255)
                                    acaACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGtttt                                                                                                           >  1:542833/1‑151 (MQ=255)
                                       aCTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTtggc                                                                                                        <  1:577324/151‑1 (MQ=255)
                                          aCCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTtggct                                                                                                       <  2:100789/149‑1 (MQ=255)
                                           cccGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTtggctggc                                                                                                    >  2:228733/1‑151 (MQ=255)
                                             cGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCg                                                                                                   >  1:431201/1‑150 (MQ=255)
                                                atgcatTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAAt                                                                                                >  1:504210/1‑150 (MQ=255)
                                                   catTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACg                                                                                            >  2:463258/1‑151 (MQ=255)
                                                     ttGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCt                                                                                          >  1:242667/1‑151 (MQ=255)
                                                       gCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGt                                                                                        <  1:648234/151‑1 (MQ=255)
                                                                ccGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTAc                                                                               >  1:375689/1‑151 (MQ=255)
                                                                   cATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTgg                                                                            <  2:77275/151‑1 (MQ=255)
                                                                                                         gCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGAc                                       <  1:880/150‑1 (MQ=255)
                                                                                                         gCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGAc                                       <  2:365184/150‑1 (MQ=255)
                                                                                                                       aaaCGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGACGCGGCGGACAGTGGGTTGGCTGCg                        <  2:496124/151‑1 (MQ=255)
                                                                                                                         aCGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGaa                      <  1:433845/151‑1 (MQ=255)
                                                                                                                          cGTATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAAt                     <  2:261585/151‑1 (MQ=255)
                                                                                                                            tATTTGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGt                   >  2:589078/1‑151 (MQ=255)
                                                                                                                                tGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTaa                >  1:474085/1‑150 (MQ=255)
                                                                                                                                tGCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTaa                <  1:234228/150‑1 (MQ=255)
                                                                                                                                 gCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTaaaa              <  1:660691/151‑1 (MQ=255)
                                                                                                                                 gCCGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTaaaa              <  1:43149/151‑1 (MQ=255)
                                                                                                                                   cGGATACCAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACg            <  2:343733/151‑1 (MQ=255)
                                                                                                                                         ccAACGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCTgcg      >  2:550550/1‑151 (MQ=255)
                                                                                                                                            aCGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCTGCGCa    >  2:453913/1‑150 (MQ=255)
                                                                                                                                             cGAACTGACAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCTGCGCaaa  <  1:476131/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGATGGTAGTCGACTATCTGCCCGGCGCGGTAAAAACGTATTTGCCGGATACCAACGAACTGGCAGGCTTGCTGTATTATCTACATCAACAACCACGTTTTGGCTGGCGAATAACGCTGTTGCCGTTACTGGAACTGTACTGGCAGCAAAGCGATCCGGCGCGGCGGACAGTGGGTTGGCTGCGAATGTTAAAACGTCTGCGCAAA  >  NC_000913/1163430‑1163721

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