Question: How to specify a branch in the tree automatically for branch-site model codeml?
gravatar for sckinta
4.6 years ago by
United States
sckinta570 wrote:

Hi I am trying to use codeml to evaluate positive selection using multiple nucleotide sequences alignment (phylip format) and maximum likelihood tree (newick format tree from RAxML). 

I used branch site model by specifying model = 2 and NSsites = 2 and test by using null model (fix omega) and alternative model (estimate omega).

To make certain branch as foreground branch, I need to add #1 on that branch in tree files. I know which taxa I am interested in. If I want to label the branch representing smallest monophyletic group which contains the taxa i am interested in (see example below), how to do it automatically? Since I have hundreds of tree to run codeml, an automatic way to label trees is necessary. Also checking tree file that contains hundreds of taxa to find the branch is painful.

(    (A37:0.0392824443,A36:0.2301499674):0.0070550356,    (A35:0.0345817897,        (A34:0.0028844901,A33:0.0086649064):0.0413749361    ):0.0479157566)

For example, I am in interested in A35. The smallest monophyletic group contains A35: is (A35:0.0345817897,        (A34:0.0028844901,A33:0.0086649064):0.0413749361    ). so I will add #1 behind that right parenthness. so the finally labeling should be

(    (A37:0.0392824443,A36:0.2301499674):0.0070550356,    (A35:0.0345817897,        (A34:0.0028844901,A33:0.0086649064)#1:0.0413749361    ):0.0479157566).

But how can do it automatically?


codeml tree • 2.6k views
ADD COMMENTlink written 4.6 years ago by sckinta570

the following Codeml wrapper could interest you: (Easy branch selection and inline visualization of results)

ADD REPLYlink modified 6 months ago by RamRS27k • written 4.6 years ago by jhc2.8k

also available as an interactive command line tool:

ADD REPLYlink modified 5 months ago by RamRS27k • written 4.5 years ago by jhc2.8k
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: 875 users visited in the last hour