Nucleotide Substitution Calculators
0 calculators tagged with “Nucleotide Substitution”
All Calculators
No calculators found for this topic.
Types of Nucleotide Substitutions
- Transitions (Ti): Purine ↔ Purine (A↔G) or Pyrimidine ↔ Pyrimidine (C↔T)
- Transversions (Tv): Purine ↔ Pyrimidine (A↔C, A↔T, G↔C, G↔T)
- Synonymous (silent): Changes codon but not amino acid
- Nonsynonymous: Changes the encoded amino acid
Substitution Rate Models
Jukes-Cantor (JC69)
Equal base frequencies, equal rates for all substitution types. Simplest model. Corrected distance: d = −(3/4)ln(1 − 4p/3)
Kimura 2-Parameter (K80)
Transitions and transversions at different rates. κ = Ti/Tv rate ratio. More biologically realistic than JC.
GTR (General Time Reversible)
Six independent substitution rate parameters + unequal base frequencies. Most parameter-rich common model — best fits most empirical datasets.
dN/dS (ω) as a Selection Indicator
dS = synonymous substitution rate (neutral proxy); dN = nonsynonymous rate.
ω = dN/dS
- ω < 1: Purifying selection — nonsynonymous mutations are deleterious
- ω = 1: Neutral evolution
- ω > 1: Positive (adaptive) selection — amino acid change is favored
Glossary
Frequently Asked Questions
A nucleotide substitution model describes the rates at which nucleotides change from one to another over evolutionary time. Simple models (JC69) assume equal substitution rates; complex models (K80, HKY85, GTR) account for unequal base frequencies, transition-transversion bias, and rate variation among sites. Models correct for multiple substitutions at the same site and are used to accurately estimate evolutionary distances and build phylogenetic trees.
A synonymous (silent) substitution changes a codon but not the encoded amino acid — due to genetic code degeneracy (e.g., CGA and CGU both encode Arginine). A nonsynonymous substitution changes the codon to one encoding a different amino acid. Synonymous rate (dS) estimates neutral evolutionary rate; nonsynonymous rate (dN) reflects substitutions changing protein sequence. The dN/dS ratio reveals whether selection acts on a gene.
ω = dN/dS. ω < 1: purifying selection removes deleterious amino acid changes — most protein-coding genes show ω 0.1–0.3. ω = 1: neutral evolution, no selective constraint on protein sequence. ω > 1: positive (adaptive) selection — amino acid change is favored faster than neutral expectation. ω > 1 is found in immune genes, pathogen-host interaction proteins, and rapidly evolving paralogs.
Model selection tests which substitution model best fits sequence data using AIC, BIC, or likelihood ratio tests. Programs like ModelFinder evaluate many models (JC69 through GTR+Γ+I). The best-fitting model is used for phylogenetic reconstruction to avoid model misspecification bias — using too simple a model can produce artifactual groupings (e.g., long-branch attraction in parsimony analyses).