Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 638,614 G→T 100% T14K (ACA→AAA)  thiD ← phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730638,6140GT100.0% 74.7 / ‑6.5 26T14K (ACA→AAA) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (15/10);  total (16/10)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGA  >  CP000730/638477‑638740
                                                                                                                                          |                                                                                                                              
cATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttgt                                                                                                                             >  1:110419/1‑141 (MQ=255)
 ctCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                            <  1:288887/140‑1 (MQ=255)
 cgcccgtGGCAATGGTGTAACATCGTGTTACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                            <  1:22541/135‑1 (MQ=255)
 ccccACTGGTAATGGTGTAACATCGTGCTACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGCAACGTTCTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                            <  1:22478/139‑1 (MQ=255)
    cATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCTACCATGCCATAGGTATCGAATTCTTGGAACGTTTTCAAATCTGCTTGCATCCCTGCGCCAGCACTTGTGTCAGa                                                                                                                         >  1:169783/1‑141 (MQ=255)
      ttGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCTGAAc                                                                                                                       >  2:157681/1‑141 (MQ=255)
          tAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATAGCACACGTTACTAATTCTTGGAACGGTTTCAAATCTGCTTGCATACCAGCGCCAGCACTTTTTTAAGAACCGGc                                                                                                                   >  2:230662/1‑141 (MQ=255)
                     cATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTaaaa                                                                                                        >  1:345449/1‑141 (MQ=255)
                     cATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGACCGAACCGGCAATTGTTaaaa                                                                                                        >  1:169494/1‑141 (MQ=255)
                      aTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCTGCAATTGTTAAAAc                                                                                                       <  1:199433/141‑1 (MQ=255)
                             gACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttcttt                                                                                                <  1:118846/141‑1 (MQ=255)
                              aCCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTa                                                                                               >  1:182970/1‑141 (MQ=255)
                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGc                                                                                   >  2:161448/1‑141 (MQ=255)
                                            tCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTCAAACTTTCTTTAAAGCCATTCAGCtt                                                                                 >  2:4702/1‑141 (MQ=255)
                                                 aGTAACGGGGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTTGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTc                                                                            <  2:102895/141‑1 (MQ=255)
                                                 aGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGACAGCTTGCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTc                                                                            <  1:296695/141‑1 (MQ=255)
                                                  gTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTcc                                                                           >  1:155269/1‑141 (MQ=255)
                                                     aCGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCCGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTGCACTCCTAc                                                                        >  2:262060/1‑141 (MQ=255)
                                                      cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTTCTTGCATACCTGCGCCAGCACTTTTGTCAGAATCGGGAATTGTTAAAACTTTCTTTAAAGACATTGAGATTCACTCGTACa                                                                       >  2:156547/1‑141 (MQ=255)
                                                       gATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCTTACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCActactaaaa                                                                      >  2:68234/1‑138 (MQ=255)
                                                            cAGTTAAAGCGACCATGCCATAAGCATTAATGTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                 >  1:91649/1‑141 (MQ=255)
                                                                           tGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                  >  1:372266/1‑141 (MQ=255)
                                                                             ctatccGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatat                                                <  1:122609/136‑1 (MQ=255)
                                                                                               ttGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTg                               <  2:110419/141‑1 (MQ=255)
                                                                                                      gTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTCCATAATAATATTGTATTCATCATTTCATTTTTAACTTAATGGaaaactt                        >  1:97171/1‑138 (MQ=255)
                                                                                                           cAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAg                   >  1:105447/1‑141 (MQ=255)
                                                                                                                            cTGCGCCAGCATGTTTGCTAGACCCGGCTATTGTTATAACTTTCTTTAAAGCCATTGAGTTTCACTACTACATAATAATATTGTATTCATCATATCATTTTTAGGCTGATTGAAAAATATTACGTATTCGATATTtgatga  <  2:112797/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                              
CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGA  >  CP000730/638477‑638740

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