From b4aed119b0fe9ce7482bd585a47a4d87a0834a86 Mon Sep 17 00:00:00 2001
From: elserj <elserj@localhost>
Date: Thu, 23 Sep 2010 22:02:41 +0000
Subject: [PATCH] Fix gene header scope error and remove > symbol from gene
 header

svn path=/; revision=51
---
 interactome_scripts/load_species_db.pl | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/interactome_scripts/load_species_db.pl b/interactome_scripts/load_species_db.pl
index e5aad70..0747cce 100755
--- 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 = "";
-- 
2.34.1