From: elserj <elserj@localhost>
Date: Thu, 23 Sep 2010 22:02:41 +0000 (+0000)
Subject: Fix gene header scope error and remove > symbol from gene header
X-Git-Url: http://gitweb.planteome.org/?a=commitdiff_plain;h=b4aed119b0fe9ce7482bd585a47a4d87a0834a86;p=old-jaiswallab-svn%2F.git

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 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 = "";