forked from davidaray/bioinfo_tools
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathgeneric.sge.sh
39 lines (29 loc) · 1020 Bytes
/
generic.sge.sh
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
#!/bin/bash
#$ -V
#$ -cwd
#$ -S /bin/bash
#$ -N SGE
#$ -o $JOB_NAME.o$JOB_ID
#$ -e $JOB_NAME.e$JOB_ID
#$ -q R2D2
#$ -pe sm 1
#$ -P communitycluster
#Preparation:
#Create a working directory
mkdir <path to working directory>
#Migrate the the working directory
cd <path to working directory>
#Create a soft link to the genome you plan to analyze in the working directory make sure extension of the file is .fa.
ln -s <path to genome file> .
#Create abbreviation for genome name
GENOME_NAME=$(basename <genome file> .fa)
/lustre/work/daray/software/faToTwoBit $GENOME_NAME".fa" $GENOME_NAME".2bit"
#run the generate RM cluster script and set to run against RepBase
/lustre/work/daray/software/generateSGEClusterRun_Chewie.pl \
-nolow \
-twoBit $GENOME".2bit" \
-batch_count 10 \
-species #use -lib option if you have a custom library, otherwise decide the appropriate species designation
#change the permissions and submit the qsub script
chmod u+x qsub.sh doLift.sh
#./qsub.sh