1 |
gpertea |
23 |
#!/bin/tcsh -f |
2 |
|
|
# |
3 |
|
|
# 1 2 3 4 5 6 7 |
4 |
|
|
#my ($file, $numpass, $slice_num, $last, $skipped, $total, $user)=@ARGV; |
5 |
|
|
set fq=$1 |
6 |
|
|
set numpass=$2 |
7 |
|
|
set slnum=$3 |
8 |
|
|
set islast=$4 |
9 |
|
|
set slskipped=$5 |
10 |
|
|
set sltotal=$6 |
11 |
|
|
set fesi=$7 |
12 |
|
|
|
13 |
|
|
set esf=$fesi:t |
14 |
|
|
set espath=$fesi:h |
15 |
|
|
set fb=$esf:r |
16 |
|
|
set fb=$fb:r |
17 |
|
|
set srvlog=srv_${fb}_$slnum.log |
18 |
|
|
set exout=$fq.exout |
19 |
|
|
|
20 |
|
|
#goto TEST_LOOP |
21 |
|
|
|
22 |
|
|
set sv=`ps auxwww | grep exonerate-server | grep 'port 3804'` |
23 |
|
|
|
24 |
|
|
if ( $#sv > 2 ) then |
25 |
|
|
echo "Existing server on this machine on port 3804, aborting.." |
26 |
|
|
exit 241 |
27 |
|
|
endif |
28 |
|
|
|
29 |
|
|
# echo "Launching server.." |
30 |
|
|
set cpath = `pwd` |
31 |
|
|
cd $espath |
32 |
|
|
exonerate-server --maxconnections 1 --verbosity 1 --proteinwordlen 4 --port 3804 $esf >& $cpath/$srvlog & |
33 |
|
|
set bpid = $! |
34 |
|
|
#echo "Background process ID= $bpid" |
35 |
|
|
|
36 |
|
|
cd $cpath |
37 |
|
|
|
38 |
|
|
while ( `grep listening $srvlog`"x" == "x" ) |
39 |
|
|
sleep 2 |
40 |
|
|
end |
41 |
|
|
sleep 2 |
42 |
|
|
echo "Running: pexo_slice.pl -o $exout $fq" |
43 |
|
|
pexo_slice.pl -o $exout $fq |
44 |
|
|
/bin/rm -f $fq |
45 |
|
|
fltexonerate -p 75 -c 70 < $exout > $exout.p75c70.gff3 |
46 |
|
|
echo '<Done.' >> /dev/stderr |
47 |
|
|
sleep 2 |
48 |
|
|
|
49 |
|
|
kill $bpid |
50 |
|
|
|
51 |
|
|
exit 0 |