Clustal errors
  • I have been trying to get clustal to run with bioperl for what seems like ages. A couple of the ways I have tried are posted here
    <!-- m --><a class="postlink" href=""></a><!-- m -->

    Please take a quick look if you have gotten clustal to run.

    OR if you could post an example of a working clustal call...
  • It is difficult to test the code without the input files. The third error is a compilation error. You need a semi-colon after this line:
    use Bio::Align::AlignI;
  • Having added the semicolon , I now just get a warning that the call crashed:
    MSG: Clustalw call ( align -infile=/tmp/bVAPEIWTAf/osmBXTWzVC -output=gcg -ktuple=1 -outfile=/tmp/bVAPEIWTAf/3IeGnfpb5V >/dev/null 2>/dev/null) crashed: 32512

    I just have no idea why.

    try.fasta is standard fasta:

  • for the first one, the exception complains about the 'output' tag. Remove it,
    @params = ('ktuple' =&gt; 2, 'matrix' =&gt; 'BLOSUM', 'outfile' =&gt; 'clustal_out&#46;fa');
    then uncomment the line of code:
      #$inputfilename = '/home/users/ostman/genomes/drosophila/clustal/try&#46;fasta';
    and it works.

    Running the third script I got a different exception message...
    MSG: Got a sequence with no letters in - cannot guess alphabet []
    STACK Bio::PrimarySeq::_guess_alphabet /usr/lib/perl5/site_perl/5.8.0/Bio/
    STACK Bio::PrimarySeq::seq /usr/lib/perl5/site_perl/5.8.0/Bio/
    STACK Bio::PrimarySeq::new /usr/lib/perl5/site_perl/5.8.0/Bio/
    STACK Bio::LocatableSeq::new /usr/lib/perl5/site_perl/5.8.0/Bio/
    STACK Bio::AlignIO::msf::next_aln /usr/lib/perl5/site_perl/5.8.0/Bio/AlignIO/
    STACK Bio::Tools::Run::Alignment::Clustalw::_run /usr/lib/perl5/site_perl/5.8.0/Bio/Tools/Run/Alignment/
    STACK Bio::Tools::Run::Alignment::Clustalw::align /usr/lib/perl5/site_perl/5.8.0/Bio/Tools/Run/Alignment/
    STACK main::vary_params
    STACK toplevel

    In the vary_params subroutine uncomment the line of code:
    #        $string = $aln-&gt;consensus_string(); # Get consensus of alignment
    and it works.