|
|
@@ -10,6 +10,7 @@ task bcftools { |
|
|
|
set -e |
|
|
|
mkdir -p /cromwell_root/tmp/bcftools |
|
|
|
|
|
|
|
echo $vcf > b |
|
|
|
for i in $vcf |
|
|
|
do |
|
|
|
bcftools view $i -Oz -o /cromwell_root/tmp/bcftools/$(basename $i.gz) |
|
|
@@ -18,7 +19,7 @@ task bcftools { |
|
|
|
|
|
|
|
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 |
|
|
|
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 |
|
|
|
>>> |
|
|
|
|
|
|
|
runtime { |
|
|
@@ -31,5 +32,6 @@ task bcftools { |
|
|
|
output { |
|
|
|
File panel_of_normal_vcf = "panel_of_normal.vcf" |
|
|
|
File a = "a" |
|
|
|
File b = "b" |
|
|
|
} |
|
|
|
} |