[Bioperl-l] Could not go on with Bio::Tools::Run::Alignment::Pal2Nal

galiantus at gmail.com galiantus at gmail.com
Tue Jun 4 05:12:45 UTC 2013


Dear colleagues!

I have a problem with running pal2nal wrapper. Here is entaire script, wich 
I try to use:

#!/usr/bin/perl
BEGIN { $ENV{PAL2NALDIR} = 'C:\KaKs\pal2nal.v14.tar\pal2nal.v14' }
use Bio::Tools::Run::Alignment::Pal2Nal;
use Bio::SeqIO;
use Bio::AlignIO;

my $alignio = Bio::AlignIO->new(-file => 'pep.aln',
                 -format => 'Clustalw');    
my $alnpr = $alignio->next_aln;

#Input DNA sequences
$seqio = Bio::SeqIO->new(-file => 'nuc.fasta',  '-format' => 'fasta');
    
    my $DNA1 = $seqio->next_seq();
    my $DNA2 = $seqio->next_seq();

# Make a Pal2Nal factory
$factory = Bio::Tools::Run::Alignment::Pal2Nal->new();

my $alnDP = $factory->run($alnpr, [$DNA1,$DNA2]);

$out = Bio::AlignIO->new(-file => ">4paml.phy",
                         -format => "phylip" );
$out->write_aln($aln);

The warning message I recieved is:

--------------------- WARNING ---------------------
MSG: Must provide a Bio::Align::AlignI object when calling write_aln
---------------------------------------------------

What is going wrong?

I also attached both data files.

Thanks for your kindly help!

Vladimir Popov
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nuc.fasta
Type: application/octet-stream
Size: 4929 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20130603/b1653553/attachment-0008.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pep.aln
Type: application/octet-stream
Size: 3259 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20130603/b1653553/attachment-0009.obj>


More information about the Bioperl-l mailing list