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