K-Estimator 6.1

k esimator

A program to estimate the number of synonymous (Ks) and nonsynonymous substitutions (Ka) per site and the confidence intervals by Monte Carlo simulations.

Download new version 6.1v

New features in v.6.1:

1. Synonymous/Nonsynonymous Weight. In codons with two or three changes, pathways with synonymous changes now can have a different probability than those pathways with nonsynonymous changes.

2. The results of simulations to explore Confidence Intervals for Divergence estimates, also used to test Ka/Ks>1, are now more detailed.

Divergence Estimates: For noncoding regions, the program can estimate the overall (K) number of nucleotide substitutions per site using several multiple-hits at a site correcting methods: Jukes and Cantor (1-parameter; Jukes and Cantor, 1969), Kimura (2-p; Kimura, 1980), Tajima and Nei (Tajima and Nei, 1984), and Tajima (1-p, 2-p and 4-p; Tajima, 1993).

When coding regions are under analysis, K-Estimator applies the method described in Comeron (1995) to estimate Ks and Ka. This method is a modification of method of Li's (1993) and Pamilo and Bianchi's (1993) (LPB) method that better quantifies the actual number of transitions and transversions and reduces stochastic errors (see Comeron, 1995, for details and comparison to previous methods). Three genetic codes can be applied: Universal, Vertebrate mitochondrial, or Drosophila mitochondrial. Furthermore, three different options can be applied to restrict the codons that are under analysis.

Confidence Intervals: K-Estimator obtains the Confidence Intervals (C.I.) of divergence estimates (K for noncoding regions, and Ks and Ka for coding regions) by Monte Carlo simulations. Computer simulations take into account the following parameters: 1) Divergence Value; K or Ks and Ka, 2) number of nucleotides or codons, 3) the transition : transversion (alfa:beta) substitution ratio, and 4) the G+C content for noncoding regions, or the amino acid composition and G+C content at the third position of codons for coding regions. The program can also calculate the exact probability of obtaining any particular divergence value (K for noncoding regions, and Ks, Ka, and Ka/Ks for coding regions).

Citation:

Comeron JM. 1995. A method for estimating the numbers of synonymous and nonsynonymous substitutions per site. J Mol Evol. 41(6):1152-9. PubMed Google Scholar
 

Comeron JM. 1999. K-Estimator: Calculation of the number of nucleotide substitutions per site and the confidence intervals. Bioinformatics. 15(9):763-4.  PubMed Google Scholar