@@ -89,9 +89,16 @@ process manta {
8989
9090 script:
9191 """
92+ if [ -z "\$ SLURM_CPUS_ON_NODE" ]; then
93+ CPUS=1
94+ else
95+ CPUS=\$ SLURM_CPUS_ON_NODE
96+ fi
97+
98+ echo "CPUS: \$ CPUS"
9299 configManta.py --normalBam $bamfile --referenceFasta $params . ref_fasta --runDir testRun
93100 cd testRun
94- ./runWorkflow.py -m local -j \$ SLURM_CPUS_ON_NODE
101+ ./runWorkflow.py -m local -j \$ CPUS
95102 gunzip -c results/variants/diploidSV.vcf.gz > ../manta.vcf
96103 """
97104}
@@ -145,9 +152,15 @@ process fermikit {
145152
146153 script:
147154 """
148- fermi2.pl unitig -s3g -t\$ SLURM_CPUS_ON_NODE -l150 -p sample sample.fq.gz > sample.mak
155+ if [ -z "\$ SLURM_CPUS_ON_NODE" ]; then
156+ CPUS=1
157+ else
158+ CPUS=\$ SLURM_CPUS_ON_NODE
159+ fi
160+
161+ fermi2.pl unitig -s3g -t\$ CPUS -l150 -p sample sample.fq.gz > sample.mak
149162 make -f sample.mak
150- run-calling -t\$ SLURM_CPUS_ON_NODE $params . ref_fasta sample.mag.gz > calling.sh
163+ run-calling -t\$ CPUS $params . ref_fasta sample.mag.gz > calling.sh
151164 bash calling.sh
152165 vcf-sort -c sample.sv.vcf.gz > fermikit.vcf
153166 bgzip -c fermikit.vcf > fermikit.vcf.gz
@@ -356,6 +369,12 @@ process variant_effect_predictor {
356369 exit 0
357370 fi
358371
372+ if [ -z "\$ SLURM_CPUS_ON_NODE" ]; then
373+ CPUS=1
374+ else
375+ CPUS=\$ SLURM_CPUS_ON_NODE
376+ fi
377+
359378 variant_effect_predictor.pl \
360379 -i "\$ INFILE" \
361380 --format "\$ FORMAT" \
@@ -373,7 +392,7 @@ process variant_effect_predictor {
373392 --total_length \
374393 --canonical \
375394 --ccds \
376- \$ ( test "\$ SLURM_CPUS_ON_NODE " -gt 1 && echo "--fork \$ SLURM_CPUS_ON_NODE " ) \
395+ \$ ( test "\$ CPUS " -gt 1 && echo "--fork \$ CPUS " ) \
377396 --fields Consequence,Codons,Amino_acids,Gene,SYMBOL,Feature,EXON,PolyPhen,SIFT,Protein_position,BIOTYPE \
378397 --assembly "\$ ASSEMBLY" \
379398 --offline
0 commit comments