瀏覽代碼

test regions

master
YaqingLiu 4 年之前
父節點
當前提交
d0db16f39a
共有 3 個檔案被更改,包括 9 行新增3 行删除
  1. +3
    -2
      inputs
  2. +4
    -1
      tasks/Haplotyper.wdl
  3. +2
    -0
      workflow.wdl

+ 3
- 2
inputs 查看文件

"{{ project_name }}.dbsnp": "dbsnp_146.hg38.vcf", "{{ project_name }}.dbsnp": "dbsnp_146.hg38.vcf",
"{{ project_name }}.recaled_bam": "{{ recaled_bam }}", "{{ project_name }}.recaled_bam": "{{ recaled_bam }}",
"{{ project_name }}.SENTIEON_INSTALL_DIR": "/opt/sentieon-genomics", "{{ project_name }}.SENTIEON_INSTALL_DIR": "/opt/sentieon-genomics",
"{{ project_name }}.dbmills_dir": "oss://pgx-reference-data/GRCh38.d1.vd1/",
"{{ project_name }}.db_mills": "Mills_and_1000G_gold_standard.indels.hg38.vcf", "{{ project_name }}.db_mills": "Mills_and_1000G_gold_standard.indels.hg38.vcf",
"{{ project_name }}.cluster_config": "OnDemand bcs.a2.3xlarge img-ubuntu-vpc", "{{ project_name }}.cluster_config": "OnDemand bcs.a2.3xlarge img-ubuntu-vpc",
"{{ project_name }}.docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/sentieon-genomics:v2020.10.07", "{{ project_name }}.docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/sentieon-genomics:v2020.10.07",
"{{ project_name }}.dbsnp_dir": "oss://pgx-reference-data/GRCh38.d1.vd1/", "{{ project_name }}.dbsnp_dir": "oss://pgx-reference-data/GRCh38.d1.vd1/",
"{{ project_name }}.sample": "{{ sample_id }}", "{{ project_name }}.sample": "{{ sample_id }}",
"{{ project_name }}.disk_size": "500", "{{ project_name }}.disk_size": "500",
"{{ project_name }}.recaled_bam_index": "{{ recaled_bam_index }}"
"{{ project_name }}.recaled_bam_index": "{{ recaled_bam_index }}",
"{{ project_name }}.regions": "oss://pgx-reference-data/bed/cbcga/S07604514_Covered.bed",
"{{ project_name }}.interval_padding": "0"
} }



+ 4
- 1
tasks/Haplotyper.wdl 查看文件

String disk_size String disk_size


File? regions File? regions
Int? interval_padding

String INTERVAL = if ${regions} then "--interval ${regions} --interval_padding ${interval_padding}" else ""


command <<< command <<<
set -o pipefail set -o pipefail
export SENTIEON_LICENSE=192.168.0.55:8990 export SENTIEON_LICENSE=192.168.0.55:8990
nt=$(nproc) nt=$(nproc)
${SENTIEON_INSTALL_DIR}/bin/sentieon driver -r ${ref_dir}/${fasta} -t $nt \ ${SENTIEON_INSTALL_DIR}/bin/sentieon driver -r ${ref_dir}/${fasta} -t $nt \
-i ${recaled_bam} --interval ${regions} \
-i ${recaled_bam} ${INTERVAL} \
--algo Haplotyper -d ${dbsnp_dir}/${dbsnp} ${sample}_hc.vcf --algo Haplotyper -d ${dbsnp_dir}/${dbsnp} ${sample}_hc.vcf
>>> >>>

+ 2
- 0
workflow.wdl 查看文件

String cluster_config String cluster_config


File? regions File? regions
Int? interval_padding


call Haplotyper.Haplotyper as Haplotyper { call Haplotyper.Haplotyper as Haplotyper {
input: input:
fasta=fasta, fasta=fasta,
ref_dir=ref_dir, ref_dir=ref_dir,
regions=regions, regions=regions,
interval_padding=interval_padding,
recaled_bam=recaled_bam, recaled_bam=recaled_bam,
recaled_bam_index=recaled_bam_index, recaled_bam_index=recaled_bam_index,
dbsnp=dbsnp, dbsnp=dbsnp,

Loading…
取消
儲存