[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