Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
how to find location of nodes (Read 1841 times)
YaBB Newbies

I love YaBB 1G - SP1!

Posts: 11
how to find location of nodes
Nov 8th, 2012 at 1:10am
I am using HYPHYMP command line to search for lineages with episodic selection.
1) do I need to find recombination breakpoints before the analysis?
2) In the results file there is also information about the nodes. How do I find where the nodes are located . Is there any tree file some where with indication of the nodes?
3) the ps output is hard to read. labels overlap and not the hole tree is in the figure. (see attached file)

Thank you
Back to top
IP Logged
YaBB Administrator

Datamonkeys are forever...

Posts: 1658
Gender: male
Re: how to find location of nodes
Reply #1 - Nov 9th, 2012 at 3:23pm
Hi hanan,

1). Branch-site REL assumes that the alignment you feed it has no recombinants. You could screen the alignment to ensure that that's the case.
2). Internal node names are assigned automatically -- other than the tree ps file (which can sometimes be hard to read, I'll see if I can improve the output), BS-REL does not output node names.
There are several ways to get around this.

(a). You can name the nodes you care about in the tree string, e.g.


In this case the ancestor of a and b will not be given a default name, but rather the name you specified (ab)

(b). If you paste the tree string as input to the following HyPhy batch file, it will print out a tree string with the internal node names shown.

fscanf (stdin, "String", treeString);
Topology T = treeString;
fprintf (stdout, "\n", Format (T,1,1), "\n");


Back to top

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
IP Logged