From 492eac9afe6b54f2c6144efbd8793d30ad7223f2 Mon Sep 17 00:00:00 2001 From: rajar Date: Wed, 2 Nov 2011 20:53:06 +0000 Subject: [PATCH] File removed svn path=/; revision=203 --- interpro/RUN_TMHMM.pl | 52 ------------------------------------------- 1 file changed, 52 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 6aaf9a2..0000000 --- a/interpro/RUN_TMHMM.pl +++ /dev/null @@ -1,52 +0,0 @@ -#!/usr/bin/perl -#use strict; -# create the variables -# keys is $gene_header, values is $seq -my %seq_hash; -my $gene_header; -my $saeq; -$file = "/lemma/Rajani/iprscan/Toly_new.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; - - } - -#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,">Toly.seq") || die("Cannot Open File"); - print FILE1 "$key\n"; - print FILE1 "$value"; - -$status7 = system("/lemma/SignalP/TMHMM2.0c/bin/tmhmm tmhmm Toly.seq >> Toly_tmhmmdata.txt"); - - -} - - -- 2.34.1