+++ /dev/null
-#!/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.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");
-
-
-}
-
-