Ver código fonte

remove -t

master
LUYAO REN 4 anos atrás
pai
commit
b4df32d3ce
3 arquivos alterados com 10 adições e 58 exclusões
  1. +4
    -21
      tasks/bcf2vcf.wdl
  2. +4
    -27
      tasks/delly.wdl
  3. +2
    -10
      workflow.wdl

+ 4
- 21
tasks/bcf2vcf.wdl Ver arquivo

task bcf2vcf { task bcf2vcf {
File ins_bcf
File ins_bcf_index
File del_bcf
File del_bcf_index
File dup_bcf
File dup_bcf_index
File inv_bcf
File inv_bcf_index
File bnd_bcf
File bnd_bcf_index

File bcf
File bcf_index
String sample_name String sample_name
String docker String docker
String disk_size String disk_size


set -o pipefail set -o pipefail
set -e set -e
/opt/hall-lab/bcftools-1.9/bin/bcftools view ${ins_bcf} > ${sample_name}.ins.vcf
/opt/hall-lab/bcftools-1.9/bin/bcftools view ${del_bcf} > ${sample_name}.del.vcf
/opt/hall-lab/bcftools-1.9/bin/bcftools view ${dup_bcf} > ${sample_name}.dup.vcf
/opt/hall-lab/bcftools-1.9/bin/bcftools view ${inv_bcf} > ${sample_name}.inv.vcf
/opt/hall-lab/bcftools-1.9/bin/bcftools view ${bnd_bcf} > ${sample_name}.bnd.vcf
/opt/hall-lab/bcftools-1.9/bin/bcftools view ${bcf} > ${sample_name}.vcf
>>> >>>


runtime { runtime {
} }


output { output {
File ins_vcf = "${sample_name}.ins.vcf"
File del_vcf = "${sample_name}.del.vcf"
File dup_vcf = "${sample_name}.dup.vcf"
File inv_vcf = "${sample_name}.inv.vcf"
File bnd_vcf = "${sample_name}.bnd.vcf"
File ins_vcf = "${sample_name}.vcf"
} }
} }



+ 4
- 27
tasks/delly.wdl Ver arquivo

nt=$(nproc) nt=$(nproc)
export OMP_NUM_THREADS=$nt export OMP_NUM_THREADS=$nt


/opt/delly/src/delly call -t DEL -g ${ref_dir}/${fasta} -o ${sample_name}.delly.del.somatic.bcf -x /opt/delly/excludeTemplates/human.hg38.excl.tsv ${tumor_bam} ${normal_bam}
/opt/delly/src/delly call -g ${ref_dir}/${fasta} -o ${sample_name}.delly.somatic.bcf -x /opt/delly/excludeTemplates/human.hg38.excl.tsv ${tumor_bam} ${normal_bam}


/opt/delly/src/delly filter -t DEL -f somatic -o ${sample_name}.delly.del.somatic.filtered.bcf -s samples.tsv ${sample_name}.delly.del.somatic.bcf
/opt/delly/src/delly filter -f somatic -o ${sample_name}.delly.somatic.filtered.bcf -s samples.tsv ${sample_name}.delly.somatic.bcf


/opt/delly/src/delly call -t DUP -g ${ref_dir}/${fasta} -o ${sample_name}.delly.dup.somatic.bcf -x /opt/delly/excludeTemplates/human.hg38.excl.tsv ${tumor_bam} ${normal_bam}

/opt/delly/src/delly filter -t DUP -f somatic -o ${sample_name}.delly.dup.somatic.filtered.bcf -s samples.tsv ${sample_name}.delly.dup.somatic.bcf

/opt/delly/src/delly call -t INS -g ${ref_dir}/${fasta} -o ${sample_name}.delly.ins.somatic.bcf -x /opt/delly/excludeTemplates/human.hg38.excl.tsv ${tumor_bam} ${normal_bam}

/opt/delly/src/delly filter -t INS -f somatic -o ${sample_name}.delly.ins.somatic.filtered.bcf -s samples.tsv ${sample_name}.delly.ins.somatic.bcf

/opt/delly/src/delly call -t INV -g ${ref_dir}/${fasta} -o ${sample_name}.delly.inv.somatic.bcf -x /opt/delly/excludeTemplates/human.hg38.excl.tsv ${tumor_bam} ${normal_bam}

/opt/delly/src/delly filter -t INV -f somatic -o ${sample_name}.delly.inv.somatic.filtered.bcf -s samples.tsv ${sample_name}.delly.inv.somatic.bcf

/opt/delly/src/delly call -t BND -g ${ref_dir}/${fasta} -o ${sample_name}.delly.bnd.somatic.bcf -x /opt/delly/excludeTemplates/human.hg38.excl.tsv ${tumor_bam} ${normal_bam}

/opt/delly/src/delly filter -t BND -f somatic -o ${sample_name}.delly.bnd.somatic.filtered.bcf -s samples.tsv ${sample_name}.delly.bnd.somatic.bcf


>>> >>>


} }


output { output {
File del_bcf = "${sample_name}.delly.del.somatic.filtered.bcf"
File del_bcf_index = "${sample_name}.delly.del.somatic.filtered.bcf.csi"
File dup_bcf = "${sample_name}.delly.dup.somatic.filtered.bcf"
File dup_bcf_index = "${sample_name}.delly.dup.somatic.filtered.bcf.csi"
File ins_bcf = "${sample_name}.delly.ins.somatic.filtered.bcf"
File ins_bcf_index = "${sample_name}.delly.ins.somatic.filtered.bcf.csi"
File inv_bcf = "${sample_name}.delly.inv.somatic.filtered.bcf"
File inv_bcf_index = "${sample_name}.delly.inv.somatic.filtered.bcf.csi"
File bnd_bcf = "${sample_name}.delly.bnd.somatic.filtered.bcf"
File bnd_bcf_index = "${sample_name}.delly.bnd.somatic.filtered.bcf.csi"
File bcf = "${sample_name}.delly.del.somatic.filtered.bcf"
File bcf_index = "${sample_name}.delly.del.somatic.filtered.bcf.csi"
} }
} }



+ 2
- 10
workflow.wdl Ver arquivo



call bcf2vcf.bcf2vcf as bcf2vcf { call bcf2vcf.bcf2vcf as bcf2vcf {
input: input:
ins_bcf=delly.ins_bcf,
ins_bcf_index=delly.ins_bcf_index,
del_bcf=delly.del_bcf,
del_bcf_index=delly.del_bcf_index,
dup_bcf=delly.dup_bcf,
dup_bcf_index=delly.dup_bcf_index,
inv_bcf=delly.inv_bcf,
inv_bcf_index=delly.inv_bcf_index,
bnd_bcf=delly.bnd_bcf,
bnd_bcf_index=delly.bnd_bcf_index,
bcf=delly.bcf,
bcf_index=delly.bcf_index,
docker=BCFdocker, docker=BCFdocker,
disk_size=disk_size, disk_size=disk_size,
cluster_config=SMALLcluster_config, cluster_config=SMALLcluster_config,

Carregando…
Cancelar
Salvar