[Bioperl-l] Bio::Tools::Run::Phylo::Phyml, tree_string

Tristan Lefebure tristan.lefebure at gmail.com
Tue Jul 26 14:47:43 UTC 2011


Hi there,
I am not quite sure I understand why tree_string() from Bio::Tools::Run::Phylo::Phyml returns 
a string that looks like that (I removed the end of the tree):

Tree is BIONJ(((((((('92':0.0114354726,'12':0.0472591023)0.0000000000:0.0000005859,...

Why do we have this 'Tree is BIONJ' thing?

A quick look at the code in the _run() function gives :

	{
	open(my $FH_TREE, "<", $tree_file)
	    || $self->throw("Phyml call ($command) did not give an output: $?");
	local $/;
	$self->{_tree} .= <$FH_TREE>;
    }

Why appending something to $self->{_tree}? What about?
	$self->{_tree} = <$FH_TREE>;

I was about to fill a bug report, but then I saw that in Phyml.t:

    is substr($factory->tree_string, 0, 9), 'BIONJ(SIN', 'tree_string()';

Well, I am lost. Any help much appreciated...

--
Tristan



More information about the Bioperl-l mailing list