diff --git a/variants/svcalling.sh b/variants/svcalling.sh index 0c3b97d85cd3c35dc2a7c739e54aef0d8124755e..b737b98abb0210a828a0ed67ab765ce23b03fa0c 100755 --- a/variants/svcalling.sh +++ b/variants/svcalling.sh @@ -110,7 +110,11 @@ then then svaba run -p $NPROC -G ${reffa} -t ${sbam} -n ${normal} -a ${pair_id} else - svaba run -p $NPROC -G ${reffa} -t ${sbam} -a ${pair_id} + bamlist='' + for i in *.bam; do + bamlist="$bamlist -t ${i}" + done + svaba run -p $NPROC -G ${reffa} -a ${pair_id} $bamlist fi #Create SV FILE vcf-concat ${pair_id}.svaba.unfiltered*sv.vcf | perl -pe 's/\.consensus|\.bam//g' | vcf-sort| bgzip > ${pair_id}.svaba.unfiltered.sv.vcf.gz