ok.Here is another thought fragment: I don't like ....let me think how this is called.... ah yes .. i don't like Mr Clever guys. this is what i absolutely hate in the scientific field.

i just need to know if it is any way to extract the information of $treeout to a file.thats all.

ok.Here is another thought fragment: I don't like ....let me think how this is called.... ah yes .. i don't like Mr Clever guys. this is what i absolutely hate in the scientific field. i just need to know if it is any way to extract the information of $treeout to a file.thats all. ps check my exp

:)

A: yada yada

B: What are you talking about, show code?

A: yada yada left

B: So what are you talking about, show code?

A: I just need help smartass

B: I DONT KNOW WHAT YOURE ASKING, SHOW CODE?

#!/usr/bin/perl

$ENV{CLUSTALDIR} = 'C:\\\\\\\\\\\\\\\\\\;
use Bio::Tools::Run::Alignment::Clustalw;

my $file = <>; # Get file name (fasta file) from command prompt.

#  Build a clustalw alignment factory
my $factory = Bio::Tools::Run::Alignment::Clustalw->new(-matrix => 'BL
+OSUM',-outfile=>'out', -output =>'clustalw'); 
my $ktuple = 2;
$factory->ktuple($ktuple);

my $inseq = Bio::SeqIO->new(
                    -file=> "<$file",
                    -format=> 'fasta'
                   );
my $seq;
my @seq_array;
while ($seq = $inseq->next_seq) {
    push(@seq_array, $seq);
}
# Now we do the actual alignment.
my $seq_array_ref = \@seq_array;
my $aln = $factory->align($seq_array_ref);

#constarcting the Tree
my $alnio = Bio::AlignIO->new(-file => 'out', -format=>'clustalw');

my $dfactory = Bio::Tree::DistanceFactory->new(-method => 'NJ');

my $stats = Bio::Align::ProteinStatistics->new;

my $treeout = Bio::TreeIO->new(-format => 'newick');

while( my $aln = $alnio->next_aln ) 
{
 my $mat = $stats->distance(-method => 'Kimura',
                            -align  => $aln);
 my $tree = $dfactory->make_tree($mat);
 $treeout->write_tree($tree);
}
my $outfile = 'C:\Users\Ioulia\Desktop\outree.tre';
my $treeout = Bio::TreeIO->new(-format => 'newick', -file => ">$outfil
+e");
exit;
[download]