Question: Calculate dN/dS ratio for protein coding genes using PAML
gravatar for b.shrestha
23 months ago by
b.shrestha0 wrote:


I am trying to calculate dN/dS ratio for protein coding genes. I have sequence alignment for a protein coding gene for about 30 species that includes an outgroup. I do not have much experience on using PAML package. I am familiar with tweaking parameters in codeml.ctl file for different models - site, branch, and branch site. However, I am not sure what kind of analysis would be better to get dN and dS values for each species in the alignment. I tired model = 1 and got branch specific dN/dS values but I am not sure branch specific omega values corresponds to omega value for each species. Or would it be better to do pairwise comparison using runmode = -2 and look at dN/dS value compare to outgroup?

It would be great to get some help with interpreting results.

Thank you!

ADD COMMENTlink written 23 months ago by b.shrestha0
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: 1365 users visited in the last hour