LUYAO REN 3 лет назад
Родитель
Сommit
cd27e588ff
3 измененных файлов: 8 добавлений и 8 удалений
  1. +3
    -3
      tasks/benchmark.wdl
  2. +1
    -1
      tasks/filter_vcf_bed.wdl
  3. +4
    -4
      workflow.wdl

+ 3
- 3
tasks/benchmark.wdl Просмотреть файл

task benchmark { task benchmark {
File vcf
File filtered_vcf
File benchmarking_dir File benchmarking_dir
File ref_dir File ref_dir
File? qc_bed File? qc_bed
String sample = basename(vcf,".hc.vcf")
String sample = basename(filtered_vcf,".filtered.vcf")
String fasta String fasta
String docker String docker
String cluster_config String cluster_config


export HGREF=/cromwell_root/tmp/reference_data/GRCh38.d1.vd1.fa 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 bgzip ${filtered_vcf} -c > ${sample}.rtg.vcf.gz
/opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg index -f vcf ${sample}.rtg.vcf.gz /opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg index -f vcf ${sample}.rtg.vcf.gz


if [ ${qc_bed} ];then if [ ${qc_bed} ];then

+ 1
- 1
tasks/filter_vcf_bed.wdl Просмотреть файл

} }
output { output {
File filtered_vcf = "${project}.filtered.vcf" File filtered_vcf = "${project}.filtered.vcf"
File? filtered_bed = "benchmark_region_query_bed.bed"
Array[File] filtered_bed = glob("benchmark_region_query_bed.bed")
} }
} }

+ 4
- 4
workflow.wdl Просмотреть файл



call benchmark.benchmark as benchmark_D5_vcf { call benchmark.benchmark as benchmark_D5_vcf {
input: input:
vcf=filter_vcf_bed_D5_vcf.filtered_vcf,
filtered_vcf=filter_vcf_bed_D5_vcf.filtered_vcf,
benchmarking_dir=benchmarking_dir, benchmarking_dir=benchmarking_dir,
ref_dir=ref_dir, ref_dir=ref_dir,
qc_bed=filter_vcf_bed_D5.filtered_bed, qc_bed=filter_vcf_bed_D5.filtered_bed,


call benchmark.benchmark as benchmark_D6_vcf { call benchmark.benchmark as benchmark_D6_vcf {
input: input:
vcf=filter_vcf_bed_D6.filtered_vcf,
filtered_vcf=filter_vcf_bed_D6.filtered_vcf,
benchmarking_dir=benchmarking_dir, benchmarking_dir=benchmarking_dir,
ref_dir=ref_dir, ref_dir=ref_dir,
qc_bed=filter_vcf_bed_D6.filtered_bed, qc_bed=filter_vcf_bed_D6.filtered_bed,


call benchmark.benchmark as benchmark_F7_vcf { call benchmark.benchmark as benchmark_F7_vcf {
input: input:
vcf=filter_vcf_bed_F7_vcf.filtered_vcf,
filtered_vcf=filter_vcf_bed_F7_vcf.filtered_vcf,
benchmarking_dir=benchmarking_dir, benchmarking_dir=benchmarking_dir,
ref_dir=ref_dir, ref_dir=ref_dir,
qc_bed=filter_vcf_bed_F7.filtered_bed, qc_bed=filter_vcf_bed_F7.filtered_bed,


call benchmark.benchmark as benchmark_M8_vcf { call benchmark.benchmark as benchmark_M8_vcf {
input: input:
vcf=filter_vcf_bed_M8_vcf.filtered_vcf,
filtered_vcf=filter_vcf_bed_M8_vcf.filtered_vcf,
benchmarking_dir=benchmarking_dir, benchmarking_dir=benchmarking_dir,
ref_dir=ref_dir, ref_dir=ref_dir,
qc_bed=filter_vcf_bed_M8.filtered_bed, qc_bed=filter_vcf_bed_M8.filtered_bed,

Загрузка…
Отмена
Сохранить