[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