From: rajar <rajar@localhost>
Date: Wed, 2 Nov 2011 20:14:25 +0000 (+0000)
Subject: Script to run Targetp application
X-Git-Url: http://gitweb.planteome.org/?a=commitdiff_plain;h=a10b218c9b28e35aca3240e5d04daccf1fcb2a3c;p=old-jaiswallab-svn%2F.git

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
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");
+
+}
+
+