diff --git a/variants/germline_vc.sh b/variants/germline_vc.sh index 3ce46806dfd6f476d9be23ae3b0d5c658c2de02c..cd1a1c71c716cf458f57ce7332674f7df6f2798e 100755 --- a/variants/germline_vc.sh +++ b/variants/germline_vc.sh @@ -100,7 +100,7 @@ then for i in *.bam; do bamlist="$bamlist --bam ${PWD}/${i}" done - cut -f 1 $fbsplit | parallel --delay 1 --jobs 0 "freebayes -f ${index_path}/genome.fa --min-mapping-quality 0 --min-base-quality 20 --min-coverage 10 --min-alternate-fraction 0.01 -C 3 --use-best-n-alleles 3 -r {} ${bamlist} > fb.{}.vcf" + cut -f 1 $fbsplit | parallel --memfree 2G --delay 1 --jobs 0 "freebayes -f ${index_path}/genome.fa --min-mapping-quality 0 --min-base-quality 20 --min-coverage 10 --min-alternate-fraction 0.01 -C 3 --use-best-n-alleles 3 -r {} ${bamlist} > fb.{}.vcf" vcf-concat fb.*.vcf | vcf-sort | vcf-annotate -n --fill-type | bcftools norm -c s -f ${reffa} -w 10 -O z -o ${pair_id}.fb.vcf.gz - elif [[ $algo == 'platypus' ]] then @@ -143,7 +143,15 @@ then opt='' if [[ -n $tbed ]] then - opt="--callRegions ${tbed}.gz" + if [[ -f "${tbed}.gz" ]] + then + opt="--callRegions ${tbed}.gz" + else + cp $tbed panel.bed + bgzip panel.bed + tabix panel.bed.gz + opt="--callRegions panel.bed.gz" + fi fi if [[ $rna == 1 ]] then diff --git a/variants/somatic_vc.sh b/variants/somatic_vc.sh index dcfce30f04eb6a48ff2898d036486fe012ac6bcb..94d4ece15db2520d511057201ea691207a06222f 100755 --- a/variants/somatic_vc.sh +++ b/variants/somatic_vc.sh @@ -103,14 +103,14 @@ then opt='' if [[ -n $tbed ]] then - if [[ -f $tbed ]] + if [[ -f "${tbed}.gz" ]] then opt="--callRegions ${tbed}.gz" else - cp $tbed targetpanel.bed - bgzip targetpanel.bed - tabix targetpanel.bed.gz - opt="--callRegions targetpanel.bed.gz" + cp $tbed panel.bed + bgzip panel.bed + tabix panel.bed.gz + opt="--callRegions panel.bed.gz" fi fi mkdir manta