|
|
@@ -10,12 +10,13 @@ task bcftools { |
|
|
|
set -e |
|
|
|
mkdir -p /cromwell_root/tmp/bcftools |
|
|
|
|
|
|
|
for i in ${sep=" " vcf} |
|
|
|
for i in $vcf |
|
|
|
do |
|
|
|
bcftools view $i -Oz -o /cromwell_root/tmp/bcftools/$(basename $i.gz) |
|
|
|
bcftools index /cromwell_root/tmp/bcftools/$(basename $i.gz) -f |
|
|
|
done |
|
|
|
|
|
|
|
ls /cromwell_root/tmp/bcftools > a |
|
|
|
cd /cromwell_root/tmp/bcftools |
|
|
|
bcftools merge -m all -f PASS,. --force-samples *.gz | bcftools plugin fill-AN-AC | bcftools filter -i 'SUM(AC)>1' > panel_of_normal.vcf |
|
|
|
>>> |
|
|
@@ -29,5 +30,6 @@ task bcftools { |
|
|
|
|
|
|
|
output { |
|
|
|
File panel_of_normal_vcf = "panel_of_normal.vcf" |
|
|
|
File a = "a" |
|
|
|
} |
|
|
|
} |