[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