#!/usr/bin/perl -w use strict; use Bio::AlignIO; BEGIN {$ENV{TCOFFEEDIR} = '/home/edward/MyBioTool/T-COFFEE_Ver1.37/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( GGGTGTTATTCAAGCAAAAAAA TTTGGAAGTCAATATTTTGTCG CCTTTTATCTGTTTTGACAGTC ACTGAAAAGCTTAGGAAATGGT TATTTGCAGTGATGTAATCAGC ); 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__