Fix gene header scope error and remove > symbol from gene header

svn path=/; revision=51

diff --git a/interactome_scripts/load_species_db.pl b/interactome_scripts/load_species_db.pl
index e5aad70b956ae730320747a23d0c665b3a3a5b3d..0747ccec443a5377292941ab5449d84548cb41e7 100755 (executable)
--- a/interactome_scripts/load_species_db.pl
+++ b/interactome_scripts/load_species_db.pl
@@ -69,11 +69,11 @@ foreach my $file (@files) {
        my %seq_hash; # keys is $gene, values are $seq and $gene_header
        my $gene;
        my $seq;
+       my $gene_header;
        
        open(in_file,$file);
        
        while(<in_file>) {
-               my $gene_header;
                my $line = $_;
                chomp $line;
                if($line =~ /^\>/) {
@@ -82,6 +82,7 @@ foreach my $file (@files) {
                                $seq_hash{$gene}->{'gene_header'} = $gene_header;
                        }
                        $gene_header = $line;
+                       $gene_header =~ s/^>//;
                        $gene = find_gene($gene_header,$species);
                        $gene =~ s/^\>//;
                        $seq = "";