[Bioperl-l] Inconsistency in Alignment Score from Bio::Tools::Run::Alignment::TCoffee versus actual command-line program

Edward WIJAYA ewijaya at singnet.com.sg
Wed Apr 20 08:08:04 EDT 2005


Hi,

Below was my code that does the alignment for two sequences.
But how come it gives different result compare to the actual  
command-line/web result?
Zero versus 100  (note the * marked line)

Could it be a bug in the module?

---BioPerl Result----
PAIRWISE_ALIGNMENT [No Tree]
T-COFFEE, Version_2.03(Wed Jul 11 14:38:06 PDT 2001)
Notredame, Higgins, Heringa, JMB(302)pp205-217,2000
CPU   0 sec
SCORE 0          *
NSEQ  2
LEN   15

seq_0           AAAAAAAAAAAAAAA
seq_1           AAAAAAAAAAAAAAA
                 ***************

=VERSUS=

---Command Line Result---
									   *
CLUSTAL FORMAT for T-COFFEE Version_2.03, CPU=0.00 sec, SCORE=100, Nseq=2,  
Len=15

seq_0           AAAAAAAAAAAAAAA
seq_1           AAAAAAAAAAAAAAA
                 ***************

---Web Result--
http://www.ch.embnet.org/wwwtmp/run_2316.html



__BEGIN__
#!/usr/bin/perl -w
use strict;
use Bio::AlignIO;
BEGIN {$ENV{TCOFFEEDIR} =
        '/home/edward/MyBioTool/T-COFFEE_Ver2.03/bin/'; }
use Bio::Tools::Run::Alignment::TCoffee;

my @params = (
               '-outfile' => 'tcof.out',
               '-maxlen'  => '100'
           );
my $factory = Bio::Tools::Run::Alignment::TCoffee->new(@params);

my @array = qw (AAAAAAAAAAAAAAAA
                  AAAAAAAAAAAAAAAA);

foreach  ( 0..$#array )
{
      push @seqs2, (Bio::Seq->new(-seq => $array[$_],
                                  -display_id => 'seq_'.$_
                                 ));
}

my  $seq_array_ref = \@seqs2;
my  $aln2 = $factory->align($seq_array_ref);  #final line that prints to  
STDOUT

__END__

--
Regards,
Edward WIJAYA
SINGAPORE


More information about the Bioperl-l mailing list