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% 59.4 / NA 21T14K (ACA→AAA) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/8);  total (13/8)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAA  >  CP000730/638493‑638749
                                                                                                                         |                                                                                                                                       
tgtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATAACTGCGCCAACACTTTTGTCAGAACCGGCAATTg                                                                                                                       >  1:180215/1‑140 (MQ=255)
    tcttCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTaa                                                                                                                   <  1:28362/137‑1 (MQ=255)
                   gTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGc                                                                                                    >  1:183395/1‑140 (MQ=255)
                        tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACTTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTg                                                                                               <  1:16396/140‑1 (MQ=255)
                                        aTGGCAGTTAAAGCGACCATGGCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTATTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACat                                                                               >  1:150245/1‑140 (MQ=255)
                                            cAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataata                                                                           >  2:183530/1‑140 (MQ=255)
                                            cAGTTAAAGCGAACATGACATACGTATCTAATTCATCGAACGTATAGAAATCAGCGAATATACCGGCTACAGCACTTTTTTCTGATCCGGCTATTTTTAAAACTTTCGTTAAAGCCATTGAGCTTCACTCCTACataata                                                                           <  2:287370/140‑1 (MQ=255)
                                             aGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataataa                                                                          >  1:259210/1‑140 (MQ=255)
                                                   aGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGt                                                                    >  1:221416/1‑140 (MQ=255)
                                                    gCGACAATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTCAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTa                                                                   <  1:185966/140‑1 (MQ=255)
                                                             ccATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCCGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                          >  1:143608/1‑140 (MQ=255)
                                                                       cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTaa                                                >  1:146539/1‑140 (MQ=255)
                                                                                ggAACGATTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCGTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                       >  2:160413/1‑140 (MQ=255)
                                                                                   aCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaaaa                                    >  2:275161/1‑140 (MQ=255)
                                                                                      tctgcAAATCTACTTGCATACCGGCTCCTGCACTTTTGTCATAACCTTCTATTGTTAAAAGCTTCTTTAAATCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTTATTGTAAtata                                 <  2:254397/136‑1 (MQ=255)
                                                                                               cTGCTTGCATACCTGCGCAAGCACGTTTGTCACAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTAACAAATATTAAGAAtt                        >  2:33193/1‑140 (MQ=255)
                                                                                                tGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTc                       <  2:36953/140‑1 (MQ=255)
                                                                                                      cATACCTGCGCCAGCACTTTTGTCAGAGCCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATAtt                 <  1:164214/140‑1 (MQ=255)
                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttg        >  2:99226/1‑140 (MQ=255)
                                                                                                                ccAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCACTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt       <  1:177824/140‑1 (MQ=255)
                                                                                                                     aCTTTTGTCAGAACCAGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa  >  2:200328/1‑140 (MQ=255)
                                                                                                                         |                                                                                                                                       
TGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAA  >  CP000730/638493‑638749

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: