[Bioperl-l] ClustalW Score?

Sendu Bala bix at sendu.me.uk
Thu May 31 08:34:38 UTC 2007


Jason Stajich wrote:
> Do something like:
> 
> my $fh;
> open($fh, "$commandstring |");
> my $score;
> while(<$fh>) {
>    $score = $1 if ($_ =~ /Score:(\d+)/);
> }
> close($fh);
> 
> ... then at the bottom after the alignment is created do:
> 
> $aln->score($score);
> 
> 
> There may be some more debugging b/c if you invoke the quiet => 1  
> parameter there may be an automatic ">& /dev/null" appended to the  
> end of the parameter string that you'll need to figure out how to  
> override.

Is there any particular reason for not having something along these 
lines committed to the module? Shall I go ahead and implement?



More information about the Bioperl-l mailing list