From 05c69e728e078495e571a3f002fa1bda01654f3c Mon Sep 17 00:00:00 2001 From: rajar Date: Wed, 2 Nov 2011 21:41:32 +0000 Subject: [PATCH] deleted svn path=/; revision=209 --- interpro/RUN_TMHMM.pl | 66 ------------------------------------------- 1 file changed, 66 deletions(-) delete mode 100644 interpro/RUN_TMHMM.pl diff --git a/interpro/RUN_TMHMM.pl b/interpro/RUN_TMHMM.pl deleted file mode 100644 index fa9212e..0000000 --- a/interpro/RUN_TMHMM.pl +++ /dev/null @@ -1,66 +0,0 @@ -#!/usr/bin/perl - use strict; - use warnings; - -#################################################################################### -# create the variables # -# keys is $gene_header, values is $seq # -# Just replace the word species in this file with the species that you are working.# -# species word appears 4 times in this file. # -#################################################################################### - -my %seq_hash; -my $gene_header; -my $seq; -my $file; -my $key; -my $value; -my $status; - -$file = "/lemma/Rajani/iprscan/species.fa"; -open(in_file,$file); - - while() - { - my $line = $_; - chomp $line; - - - if($line =~ /^\>/) - { - if(defined($gene_header)) - { - $seq_hash{$gene_header} = $seq; - } - $gene_header = $line; - $seq = ""; - } - else - { - $seq = "$seq"."$line"; - } - $seq_hash{$gene_header} = $seq; - - } -# To display the keys and values on the screen -# foreach $key(%seq_hash) -# { -# print "Gene is $key and the sequence is $seq_hash{$gene_header} \n"; -# } - - - - while (($key, $value) = each(%seq_hash)) - { - # print $key."\n".$value."\n"; - - open(FILE1,">species.seq") || die("Cannot Open File"); - print FILE1 "$key\n"; - print FILE1 "$value"; - - $status = system("/lemma/SignalP/TMHMM2.0c/bin/tmhmm tmhmm species.seq >> species_tmhmm.txt"); - - - } - - -- 2.34.1