hi mortiz, you always adress me to that page and i always read that but still i ahve problem.:-(
this is my main script that i used this subroutine
#!/usr/bin/perl -w
use strict;
my %genome1;
my %genome2;
my @files;
my $file;
#input first genome file name.
print "Enter the first genome fasta file=\n";
my $filename=<>;
#analyze the genome.getting domains,number of domains and their combin
+ations.
analyze_genome($filename);
%genome1=n_nn();
for my $key (keys %genome1){
print "$genome1{$key}\n";
}
#delet useless files.
@files = ("domainlist","result","first");
for my $file (@files) {
unlink $file
or warn "Couldn't remove '$file': $!";
};
#input second genome file name.
print"Enter the second genome fasta file=\n";
$filename=<>;
analyze_genome($filename);
%genome2=n_nn();
#delet usless files.
@files = ("domainlist","result","first");
for my $file (@files) {
unlink $file
or warn "Couldn't remove '$file': $!";
};
#calculate domain versatility for each domain in the genomes.
versatility(%genome1,%genome2);
exit;
|