diff --git a/alignment/dnaseqalign.sh b/alignment/dnaseqalign.sh index ffa99c473ba1a0588583188813609fde8757dae4..219f7f822c4cfd2f14c4287c73ba241f88342113 100644 --- a/alignment/dnaseqalign.sh +++ b/alignment/dnaseqalign.sh @@ -12,7 +12,7 @@ usage() { exit 1 } OPTIND=1 # Reset OPTIND -while getopts :r:x:y:p:u:h opt +while getopts :r:x:y:p:uh opt do case $opt in r) index_path=$OPTARG;; diff --git a/variants/cnvkit.sh b/variants/cnvkit.sh index d4ba922c136de31ea9d0b8440f558cca0d55e312..40b72185e46627dae8078e28777608c3e83b9453 100644 --- a/variants/cnvkit.sh +++ b/variants/cnvkit.sh @@ -21,6 +21,8 @@ done shift $(($OPTIND -1)) +index_path='/project/shared/bicf_workflow_ref/GRCh38/clinseq_prj/' + # Check for mandatory options if [[ -z $pair_id ]] || [[ -z $sbam ]]; then usage @@ -33,13 +35,13 @@ baseDir="`dirname \"$0\"`" source /etc/profile.d/modules.sh module load cnvkit/0.9.0 -cnvkit.py coverage ${sbam} /project/shared/bicf_workflow_ref/GRCh38/UTSWV2.cnvkit_targets.bed -o ${pair_id}.targetcoverage.cnn -cnvkit.py coverage ${sbam} /project/shared/bicf_workflow_ref/GRCh38/UTSWV2.cnvkit_antitargets.bed -o ${pair_id}.antitargetcoverage.cnn +cnvkit.py coverage ${sbam} ${index_path}/UTSWV2.cnvkit_targets.bed -o ${pair_id}.targetcoverage.cnn +cnvkit.py coverage ${sbam} ${index_path}/UTSWV2.cnvkit_antitargets.bed -o ${pair_id}.antitargetcoverage.cnn if [[ $umi == 'umi' ]] then -cnvkit.py fix ${pair_id}.targetcoverage.cnn ${pair_id}.antitargetcoverage.cnn /project/shared/bicf_workflow_ref/GRCh38/UTSWV2.uminormals.cnn -o ${pair_id}.cnr +cnvkit.py fix ${pair_id}.targetcoverage.cnn ${pair_id}.antitargetcoverage.cnn ${index_path}/UTSWV2.uminormals.cnn -o ${pair_id}.cnr else -cnvkit.py fix ${pair_id}.targetcoverage.cnn ${pair_id}.antitargetcoverage.cnn /project/shared/bicf_workflow_ref/GRCh38/UTSWV2.normals.cnn -o ${pair_id}.cnr +cnvkit.py fix ${pair_id}.targetcoverage.cnn ${pair_id}.antitargetcoverage.cnn ${index_path}/UTSWV2.normals.cnn -o ${pair_id}.cnr fi cnvkit.py segment ${pair_id}.cnr -o ${pair_id}.cns