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% 52.4 / ‑6.5 22T14K (ACA→AAA) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (13/7);  total (13/9)
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.

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAA  >  CP000730/638477‑638749
                                                                                                                                         |                                                                                                                                       
cATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttgt                                                                                                                                      <  2:23365/141‑1 (MQ=255)
            aTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCaa                                                                                                                          >  1:252425/1‑141 (MQ=255)
             tGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAGCGGCAAt                                                                                                                         >  1:230550/1‑141 (MQ=255)
                     cATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTaaaa                                                                                                                 <  1:91496/141‑1 (MQ=255)
                               ccATGTATCTTTATAAATAGTAACGGTGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTaa                                                                                                       >  1:209183/1‑141 (MQ=255)
                                 aTGTATCTTTATCCATAGTAACGATGGCAGTCAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCAAACCTGCGCCAGCACTAATGTCAAAACCCGCAATTGCTAAAACTTTCTTTAAAg                                                                                                     <  2:237874/141‑1 (MQ=255)
                                                 aGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTTTTGGAACGTTATCAAATCTGCTTGCATACCTTCTCCAGCACTTTTGTCAGAACATGCATTTGATAAAACATTCTTTAAAGCCATTCACCTTCACTc                                                                                     >  2:15889/1‑141 (MQ=255)
                                                  gTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTcc                                                                                    >  2:1342/1‑141 (MQ=255)
                                                            cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                          >  1:112230/1‑141 (MQ=255)
                                                            cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                          >  1:201034/1‑141 (MQ=255)
                                                             aGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCCCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTATGGCGATTGACATTCACAACTACataataat                                                                         <  2:118291/141‑1 (MQ=255)
                                                                             ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGGACCGGGAATTGTTAAAAATTTCATATCATTCGAAGAGCAAAACTACTAAGTAGTAATATTGTATTCatcatat                                                         >  1:133886/1‑141 (MQ=255)
                                                                             ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTCGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatat                                                         >  2:130649/1‑141 (MQ=255)
                                                                              cATACCTATCTAATTCGTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATATTAATATTGTATTCatcatatc                                                        <  2:48494/141‑1 (MQ=255)
                                                                                  cGTATCTAATTCTTGCAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttt                                                    <  2:217013/141‑1 (MQ=255)
                                                                                      tCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa                                                >  1:13079/1‑141 (MQ=255)
                                                                                                                    tGCATACCTGCGCCAGCACCTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAatat                  <  1:103533/141‑1 (MQ=255)
                                                                                                                          ccTGCGCCACCCCTGTTGCCAGACCCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGTTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatg            <  1:120124/141‑1 (MQ=255)
                                                                                                                           cttcgcCAGCACTTTTGTCAGAACCGTCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAATCTAATTGAAAAATATTAAGCATTCAATATTtgatga           <  1:214306/138‑1 (MQ=255)
                                                                                                                                   gCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaa   >  1:240586/1‑141 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTTTTGaaa  >  2:60237/1‑141 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTTaaa  >  2:22876/1‑141 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa  <  1:130649/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAA  >  CP000730/638477‑638749

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