[Bioperl-l] Help for using Clustalw.pm

Sun, Jian jsun at utdallas.edu
Mon Jul 26 17:27:56 EDT 2004


Dear all:
     I tried to align multiple sequence through Clustalw.pm, and I used the source code attached below:
*****************************************************************************************************
#!C:\Perl\bin\perl.exe
use lib "C:\Perl\lib";
 
use Bio::Perl;
use Bio::Tools::Run::Alignment::Clustalw;
use strict;
use warnings;
 
$ENV{CLUSTALDIR} = 'C:/Program Files/Apache Group/Apache2/bin/';
my @params = ('ktuple' => 2, 'matrix' => 'BLOSUM');
  my $factory = Bio::Tools::Run::Alignment::Clustalw->new(@params);
  my $ktuple = 3;
  $factory->ktuple($ktuple);  # change the parameter before executing
  my $str = Bio::SeqIO->new(-file=> 'Clustseq.fa', '-format' => 'Fasta');
 my @seq_array =();
 while ( my $seq = $str->next_seq() ) {push (@seq_array, $seq) ;}
  my $seq_array_ref = \@seq_array;
      # where @seq_array is an array of Bio::Seq objects
  my $aln = $factory->align($seq_array_ref);
**************************************************************************
 
and when I run the pl file, I get this error message:
 
///////////////////////////////////////////////////////////////////////////////////////////
Clustalw program not found as clustalw or not executable.
......  
'clustalw' is not recognized as an internal or external command, 
operable program or batch file.
---------------------EXCEPTION-----------------------------------------
MSG: Clustalw call <clustal align -infile=C:\...\> crashed:256
 
STACK Bio::Tools::Run::Alignment::Clustalw::_run c:/Perl/Site/Lib/Bio/Toold/Run/Alignment/Clustalw.pm:581
STACK Bio::Tools::Run::Alignment::Clustalw::_run c:/Perl/Site/Lib/Bio/Toold/Run/Alignment/Clustalw.pm:507
STACK toplevel test726pl.pl:33
-----------------------------------------------------------------------------------------
///////////////////////////////////////////////////////////////////////////////////////
 
Since I already set the CLUSTALDIR variable, I don't know why the Clustalw is still not execuatable. Did I did the right setting?  Does anyone have this kind of experience? 
 
Thanks in advance.
Jane
 



More information about the Bioperl-l mailing list