Script to run Targetp application

svn path=/; revision=202

diff --git a/interpro/RUN_TARGETP.pl b/interpro/RUN_TARGETP.pl
new file mode 100644 (file)
index 0000000..2c83d0d
--- /dev/null
+++ b/interpro/RUN_TARGETP.pl
@@ -0,0 +1,51 @@
+#!/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(<in_file>) 
+{
+        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.fa") || die("Cannot Open File");
+              print FILE1 "$key\n";
+              print FILE1 "$value";
+
+$astatus6 = system("/lemma/SignalP/targetp-1.1/targetp -P /lemma/Rajani/iprscan/Toly.fa >> targetp_4_Toly.out");
+
+}
+
+