[Bioperl-l] Getting Alignment Score from Bio::Tools::Run::Alignment::TCoffee

Edward Wijaya ewijaya at singnet.com.sg
Wed Apr 20 05:26:56 EDT 2005


Hi,
For the above TCoffee alignment module

Is there a way to:
1. get the alignment score (any method for this)
2. Pass input sequence as an array (not fasta file)
   e.g. input is "@seq = qw(ATCCGGA AGGGGAA)" etc.

Below is my current working program.

Thanks and hope to hear from you again.
Regards,
Edward WIJAYA

__BEGIN__
#!/usr/bin/perl -w
use strict;
use Bio::AlignIO;
use Data::Dumper;

BEGIN {$ENV{TCOFFEEDIR} = '/home/edward/MyBioTool/T-COFFEE_distribution_Version_2.03/bin/'; }
use Bio::Tools::Run::Alignment::TCoffee;

my $inputfilename = "hm02r.fasta";
my $outfile = "tcoffee.out";
my @params = (
              '-in'    => $inputfilename,
              '-outfile' => $outfile,
          );
my $factory = Bio::Tools::Run::Alignment::TCoffee->new(@params);
my $aln = $factory->align($inputfilename);

__END__


More information about the Bioperl-l mailing list