forked from albertwcheng/RNASeqMappingScripts3
-
Notifications
You must be signed in to change notification settings - Fork 0
/
tophatMapSE.sh
47 lines (33 loc) · 1.12 KB
/
tophatMapSE.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
38
39
40
41
42
43
44
45
#!/bin/bash
source fileUtils.sh
source ./initvars.sh
source $tophatshvar
if [ ! -e $tophatOutputDir ];then
mkdir $tophatOutputDir
fi
#requestEmptyDirWithWarning $tophatOutputDir
saveIFS=$IFS
IFS=`echo -en ","`
declare -a samples=($Samples)
IFS=$saveIFS
###nodenum=1
for sample in ${samples[*]}; do
echo "initiate tophat for sample $sample"
vname=${sample}
lfilelist=${!vname}
#lfilelist=`echo $lfilelist | tr "," "|"`
echo left file list $lfilelist
###nodehostname="episode-0$nodenum"
#qsubcommand="qsub -v SAMPLENAME=$sample,LFILELIST=\"$lfilelist\",RFILELIST=\"$rfilelist\",SCRIPTDIR=$scriptDir -q $jobQueue -e $queueJobStdWritePath/$sample.tophatqueue.err -m a -M $queueJobEmail -o $queueJobStdWritePath/$sample.tophatqueue.out $scriptDir/tophatMapJob.sh" ###-l nodes=$nodehostname
#echo $qsubcommand
#eval $qsubcommand
mkdir $tophatOutputDir/${sample}
bsubcommand="bsub $scriptDir/tophatMapJobArgListSE.sh ${sample} $lfilelist $scriptDir"
echo $bsubcommand
eval $bsubcommand
###nodenum=`expr $nodenum + 1`
#declare -a filelist=(${!vname})
#for fil in ${filelist[*]}; do
# echo "file $fil"
#done
done