Question: Calculating Kn/Ks Per Site
3
gravatar for John
8.1 years ago by
John30
John30 wrote:

Hi all,

I've used PAML on an extensive multiple sequence alignment, running pairwise comparisons. I would also like to calculate strength of selection on individual sites. Is it possible with PAML or any other quickly executable programs?

Cheers,

Michael

evolution protein selection paml • 3.5k views
ADD COMMENTlink modified 8.1 years ago by fransua390 • written 8.1 years ago by John30
2
gravatar for Philippe
8.1 years ago by
Philippe1.9k
Barcelona, Spain.
Philippe1.9k wrote:

Hi Michael,

this is possible with PAML. You have to choose between different models: branch-specific model, site-specific model or branch-site model. They are all detailed in the PAML documentation. You can also check this online course for a more user friendly help. The PAML FAQ can also be of interest for some specific questions.

By the way, the common name for the ratio is either Ka/Ks or dN/dS, but not Kn/Ks.

ADD COMMENTlink written 8.1 years ago by Philippe1.9k
1
gravatar for fransua
8.1 years ago by
fransua390
Spain
fransua390 wrote:

Hi, I am not sure you can find something about selective pressure by site doing pure pairwise analysis (runmode = -2), this is usually used for global distances. On the other hand doing site analysis you can get this info by site, but it is generally a bad idea to use this result when analyzing less than 4 sequences (if you try it, have a look to associated errors of the estimations).

have a look to this post: http://gsf.gc.ucdavis.edu/viewtopic.php?f=1&t=1099

an other program that can do this for you is SLR http://www.ebi.ac.uk/goldman/SLR/ but I am afraid that you will find same restrictions...

I think the best solution, if you can, is to group your sequences by 4.

ADD COMMENTlink written 8.1 years ago by fransua390
Please log in to add an answer.

Help
Access

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