浏览代码

reference dataset4.0

master
LUYAO REN 4 年前
父节点
当前提交
43744e524d
共有 3 个文件被更改,包括 8 次插入6 次删除
  1. +1
    -1
      inputs
  2. +5
    -5
      tasks/benchmark.wdl
  3. +2
    -0
      workflow.wdl

+ 1
- 1
inputs 查看文件

{ {
"{{ project_name }}.benchmarking_dir": "oss://pgx-result/renluyao/manuscript/benchmark_calls_v3.0/",
"{{ project_name }}.benchmarking_dir": "oss://pgx-result/renluyao/manuscript_v3.0/reference_dataset_v4.0/",
"{{ project_name }}.fasta": "GRCh38.d1.vd1.fa", "{{ project_name }}.fasta": "GRCh38.d1.vd1.fa",
"{{ project_name }}.disk_size": "500", "{{ project_name }}.disk_size": "500",
"{{ project_name }}.vcf_idx": "{{ vcf_idx }}", "{{ project_name }}.vcf_idx": "{{ vcf_idx }}",

+ 5
- 5
tasks/benchmark.wdl 查看文件

File vcf_idx File vcf_idx
File benchmarking_dir File benchmarking_dir
File ref_dir File ref_dir
String sample = basename(vcf,".rtg.vcf.gz")
String sample
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


if [[ ${sample} =~ "LCL5" ]];then if [[ ${sample} =~ "LCL5" ]];then
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL5.afterfilterdiffbed.vcf.gz ${vcf} -f ${benchmarking_dir}/LCL5.high.confidence.bed.gz --threads $nt -o ${sample} -r ${ref_dir}/${fasta}
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL5.voted.mendelian.vcf.gz ${vcf} -f ${benchmarking_dir}/LCL5.highconfidence.bed --threads $nt -o ${sample} -r ${ref_dir}/${fasta}
elif [[ ${sample} =~ "LCL6" ]]; then elif [[ ${sample} =~ "LCL6" ]]; then
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL6.afterfilterdiffbed.vcf.gz ${vcf} -f ${benchmarking_dir}/LCL6.high.confidence.bed.gz --threads $nt -o ${sample} -r ${ref_dir}/${fasta}
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL6.voted.mendelian.vcf.gz ${vcf} -f ${benchmarking_dir}/LCL6.highconfidence.bed --threads $nt -o ${sample} -r ${ref_dir}/${fasta}
elif [[ ${sample} =~ "LCL7" ]]; then elif [[ ${sample} =~ "LCL7" ]]; then
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL7.afterfilterdiffbed.vcf.gz ${vcf} -f ${benchmarking_dir}/LCL7.high.confidence.bed.gz --threads $nt -o ${sample} -r ${ref_dir}/${fasta}
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL7.voted.mendelian.vcf.gz ${vcf} -f ${benchmarking_dir}/LCL7.highconfidence.bed --threads $nt -o ${sample} -r ${ref_dir}/${fasta}
elif [[ ${sample} =~ "LCL8" ]]; then elif [[ ${sample} =~ "LCL8" ]]; then
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL8.afterfilterdiffbed.vcf.gz ${vcf} -f ${benchmarking_dir}/LCL8.high.confidence.bed.gz --threads $nt -o ${sample} -r ${ref_dir}/${fasta}
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL8.voted.mendelian.vcf.gz ${vcf} -f ${benchmarking_dir}/LCL8.highconfidence.bed --threads $nt -o ${sample} -r ${ref_dir}/${fasta}
else else
echo "only for quartet samples" echo "only for quartet samples"
fi fi

+ 2
- 0
workflow.wdl 查看文件

File vcf_idx File vcf_idx
File benchmarking_dir File benchmarking_dir
File ref_dir File ref_dir
String sample
String fasta String fasta
String docker String docker
String cluster_config String cluster_config
vcf_idx=vcf_idx, vcf_idx=vcf_idx,
benchmarking_dir=benchmarking_dir, benchmarking_dir=benchmarking_dir,
ref_dir=ref_dir, ref_dir=ref_dir,
sample=sample,
fasta=fasta, fasta=fasta,
docker=docker, docker=docker,
cluster_config=cluster_config, cluster_config=cluster_config,

正在加载...
取消
保存