|
|
@@ -14,7 +14,10 @@ task benchmark { |
|
|
|
set -o pipefail |
|
|
|
set -e |
|
|
|
nt=$(nproc) |
|
|
|
export HGREF=/cromwell_inputs/*/reference_data/GRCh38.d1.vd1.fa |
|
|
|
mkdir -p /cromwell_root/tmp |
|
|
|
cp -r ${ref_dir} /cromwell_root/tmp/ |
|
|
|
|
|
|
|
export HGREF=/cromwell_root/tmp/reference_data/GRCh38.d1.vd1.fa |
|
|
|
|
|
|
|
/opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg bgzip ${vcf} -c > ${sample}.rtg.vcf.gz |
|
|
|
/opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg index -f vcf ${sample}.rtg.vcf.gz |
|
|
@@ -40,6 +43,8 @@ task benchmark { |
|
|
|
} |
|
|
|
|
|
|
|
output { |
|
|
|
File rtg_vcf = "${sample}.rtg.vcf.gz" |
|
|
|
File rtg_vcf_index = "${sample}.rtg.vcf.gz.tbi" |
|
|
|
File gzip_vcf = "${sample}.vcf.gz" |
|
|
|
File gzip_vcf_index = "${sample}.vcf.gz.tbi" |
|
|
|
File roc_all_csv = "${sample}.roc.all.csv.gz" |