Question: Using Rate4Site To Calculate Mutation Rate Per Gene?
gravatar for anin.gregory
7.9 years ago by
United States
anin.gregory90 wrote:

I would like to calculate the mutation rate per gene based on an MSA. In this past post (C: Mutation rates for a specific gene?) someone suggested using Rate4Site to calculate mutation rate per gene. However the program is specifically designed to calculate the conservation score for each amino acid, with a lower value indicating higher conservation.

I read the associated papers Rate4Site: an algorithmic tool for the identification of functional regions in proteins by surface mapping of evolutionary determinants within their homologues and Comparison of Site-Specific Rate-Inference Methods for Protein Sequences: Empirical Bayesian Methods Are Superior, and I don't understand how you can get from the conservation score to a mutation rate. Maybe I am missing something in the papers. Has anyone tried using Rate4Site to calculate mutation rate per amino acid / per gene??

gene mutation msa • 3.6k views
ADD COMMENTlink modified 7.9 years ago by Asaf8.5k • written 7.9 years ago by anin.gregory90
gravatar for Asaf
7.9 years ago by
Asaf8.5k wrote:

The rate of each position (the r parameter which is the only parameter estimated) appears in the output sent to stdout or to a file using the option -y filename. The conservation scores are just z-scores of these rates, normalized for each gene.

ADD COMMENTlink written 7.9 years ago by Asaf8.5k

I'm sorry. I still don't understand how this is the mutation rate. It is saying that the rate is 0.9999. This does not make sense, since mutation rates are usually in 10^-8 orders of magnitude for my system. Am I doing something wrong? I an MSA and use the default NJ tree option. Thanks in advance.

ADD REPLYlink written 7.8 years ago by anin.gregory90
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2339 users visited in the last hour