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.

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGA  >  CP000730/638477‑638740
                                                                                                                                          |                                                                                                                              
cATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttgt                                                                                                                             >  1:110419/1‑141 (MQ=255)
 ctCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                            <  1:288887/140‑1 (MQ=255)
 cgcccgtGGCAATGGTGTAACATCGTGTTACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                            <  1:22541/135‑1 (MQ=255)
 ccccACTGGTAATGGTGTAACATCGTGCTACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGCAACGTTCTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTc                                                                                                                            <  1:22478/139‑1 (MQ=255)
    cATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCTACCATGCCATAGGTATCGAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATCCCTGCGCCAGCACTTGTGTCAGa                                                                                                                         >  1:169783/1‑141 (MQ=255)
      ttGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCTGAAc                                                                                                                       >  2:157681/1‑141 (MQ=255)
          tAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATAGCACACGTTACTAAT‑TCTTGGAACGGTTTCAAATCTGCTTGCATACCAGCGCCAGCACTTTTTTAAGAACCGGc                                                                                                                   >  2:230662/1‑141 (MQ=255)
                     cATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTaaaa                                                                                                        >  1:345449/1‑141 (MQ=255)
                     cATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGACCGAACCGGCAATTGTTaaaa                                                                                                        >  1:169494/1‑141 (MQ=255)
                      aTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCTGCAATTGTTAAAAc                                                                                                       <  1:199433/141‑1 (MQ=255)
                             gACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttcttt                                                                                                <  1:118846/141‑1 (MQ=255)
                              aCCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTa                                                                                               >  1:182970/1‑141 (MQ=255)
                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGc                                                                                   >  2:161448/1‑141 (MQ=255)
                                            tCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTCAAACTTTCTTTAAAGCCATTCAGCtt                                                                                 >  2:4702/1‑141 (MQ=255)
                                                 aGTAACGGGGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTTGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTc                                                                            <  2:102895/141‑1 (MQ=255)
                                                 aGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGACAGCTTGCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTc                                                                            <  1:296695/141‑1 (MQ=255)
                                                  gTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTcc                                                                           >  1:155269/1‑141 (MQ=255)
                                                     aCGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCCGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTGCACTCCTAc                                                                        >  2:262060/1‑141 (MQ=255)
                                                      cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTTCTTGCATACCTGCGCCAGCACTTTTGTCAGAATCGGGAATTGTTAAAACTTTCTTTAAAGACATTGAGATTCACTCGTACa                                                                       >  2:156547/1‑141 (MQ=255)
                                                       gATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCTTACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCActactaaaa                                                                      >  2:68234/1‑138 (MQ=255)
                                                            cAGTTAAAGCGACCATGCCATAAGCAT‑TAATGTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                 >  1:91649/1‑141 (MQ=255)
                                                                           tGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                  >  1:372266/1‑141 (MQ=255)
                                                                             ctatccGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatat                                                <  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)
                                                                                                                                          |                                                                                                                              
CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAAT‑TCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGA  >  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: ‑