-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathBatch_SCLAP
executable file
·43 lines (39 loc) · 1.1 KB
/
Batch_SCLAP
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
#!/usr/bin/perl
######## Sub Chloroplast Loclaization using ADABOOST predictor (SCLAP)#####################
# Developed by SARAVANAN VIJAYAKUMAR
# Centre for Bioinformatics, Pondicherry University
# brsaran@gmail.com
# Date: 15-03-2012
#
use List::Util qw[max];
use CGI qw(:standard);
print "Enter Protein File:";
$input = <STDIN>;
open(FIL1,"$input");
@Seq_hold1 = <FIL1>;
foreach(@Seq_hold1){
if(substr(@Seq_hold1[$j],0,1) eq ">"){
$i=$i+1;
@Seq_submit1[$i]= @Seq_hold1[$j]."\n";
@FASTA_Header[$i]= substr($Seq_hold1[$j],1,11);
$j++;
}
else{
@Rawseq1[$i].= @Seq_hold1[$j];
$Rawseq1[$i]=~s/\s//g;
@Seq_submit1[$i].= @Seq_hold1[$j];
$j++;
$k++;
}
}
for($m=1;$m<=@Rawseq1;$m++){
$x = int(rand(10002)).".fasta";
open(FL,">$x");
print FL $Rawseq1[$m];
close(FL);
$FileOh = $FASTA_Header[$m];
$FileOh=~s/\W//g;
`./SCLAP $x $FileOh`;
unlink($x);
print "$FASTA_Header[$m] Completed !\n"
}