-
Notifications
You must be signed in to change notification settings - Fork 1
/
tophatMapR5.sh
executable file
·46 lines (34 loc) · 1.17 KB
/
tophatMapR5.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
#!/bin/bash
source fileUtils.sh
source ./initvars.sh
source $tophatshvar
#requestEmptyDirWithWarning $tophatOutputDir
saveIFS=$IFS
IFS=`echo -en ","`
declare -a samples=($Samples)
IFS=$saveIFS
###nodenum=1
for sample in R5; do
echo "initiate tophat for sample $sample"
vname=${sample}_1
lfilelist=${!vname}
#lfilelist=`echo $lfilelist | tr "," "|"`
vname=${sample}_2
rfilelist=${!vname}
#rfilelist=`echo $rfilelist | tr "," "|"`
echo left file list $lfilelist
echo right file list $rfilelist
###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/tophatMapJobArgList.sh ${sample} $lfilelist $rfilelist $scriptDir"
echo $bsubcommand
eval $bsubcommand
###nodenum=`expr $nodenum + 1`
#declare -a filelist=(${!vname})
#for fil in ${filelist[*]}; do
# echo "file $fil"
#done
done