9 |
|
|
10 |
if ((! $rc) || (! -d $dataD)) { print STDERR $usage; exit } |
if ((! $rc) || (! -d $dataD)) { print STDERR $usage; exit } |
11 |
|
|
12 |
my %anno = map { chop; ($_ => 1) } `cat $dataD/anno.seed`; |
my %anno; |
13 |
|
if (-s "$dataD/anno.seed") |
14 |
|
{ |
15 |
|
%anno = map { chop; ($_ => 1) } `cat $dataD/anno.seed`; |
16 |
|
} |
17 |
|
else |
18 |
|
{ |
19 |
|
print STDERR "NOTE: you have not included any ASEED genomes in a file $dataD/anno.seed\n"; |
20 |
|
} |
21 |
|
|
22 |
my$peg_to_func = &load_funcs($dataD); |
my$peg_to_func = &load_funcs($dataD); |
23 |
my($peg_to_seq,$seq_to_pegs) = &load_seqs($dataD); |
my($peg_to_seq,$seq_to_pegs) = &load_seqs($dataD); |