YaqingLiu 4 anni fa
parent
commit
d43d762386
3 ha cambiato i file con 3 aggiunte e 36 eliminazioni
  1. +2
    -7
      tasks/bcftools.wdl
  2. +0
    -28
      tasks/bcftools1.wdl
  3. +1
    -1
      workflow.wdl

+ 2
- 7
tasks/bcftools.wdl Vedi File

set -e set -e
mkdir -p /cromwell_root/tmp/bcftools mkdir -p /cromwell_root/tmp/bcftools


echo $vcf > b
for i in $vcf
for i in ${sep=" " vcf}
do do
bcftools view $i -Oz -o /cromwell_root/tmp/bcftools/$(basename $i.gz) bcftools view $i -Oz -o /cromwell_root/tmp/bcftools/$(basename $i.gz)
bcftools index /cromwell_root/tmp/bcftools/$(basename $i.gz) -f bcftools index /cromwell_root/tmp/bcftools/$(basename $i.gz) -f
done done


ls /cromwell_root/tmp/bcftools > a
cd /cromwell_root/tmp/bcftools
bcftools merge -m all -f PASS,. --force-samples *.vcf.gz | bcftools plugin fill-AN-AC | bcftools filter -i 'SUM(AC)>1' > panel_of_normal.vcf
bcftools merge -m all -f PASS,. --force-samples /cromwell_root/tmp/bcftools/*.vcf.gz | bcftools plugin fill-AN-AC | bcftools filter -i 'SUM(AC)>1' > panel_of_normal.vcf
>>> >>>


runtime { runtime {
output { output {
File panel_of_normal_vcf = "panel_of_normal.vcf" File panel_of_normal_vcf = "panel_of_normal.vcf"
File a = "a"
File b = "b"
} }
} }

+ 0
- 28
tasks/bcftools1.wdl Vedi File

task bcftools {

Array[File] vcf
String docker
String cluster_config
String disk_size

command <<<
set -o pipefail
set -e
mkdir -p /cromwell_root/tmp/bcftools
for i in ${sep=" " vcf}
do
echo $i > b
done
>>>

runtime {
docker: docker
cluster: cluster_config
systemDisk: "cloud_ssd 40"
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
}
output {
File b = "b"
}
}

+ 1
- 1
workflow.wdl Vedi File

import "./tasks/bcftools1.wdl" as bcftools
import "./tasks/bcftools.wdl" as bcftools


workflow {{ project_name }} { workflow {{ project_name }} {

Loading…
Annulla
Salva