Browse Source

first commit

master
YaqingLiu 5 years ago
parent
commit
be69f5a27d
1 changed files with 4 additions and 1 deletions
  1. +4
    -1
      tasks/bcftools.wdl

+ 4
- 1
tasks/bcftools.wdl View File



bcftools view ${pon_vcf1} -Oz -o /cromwell_root/tmp/bcftools/$(basename $pon_vcf1).gz bcftools view ${pon_vcf1} -Oz -o /cromwell_root/tmp/bcftools/$(basename $pon_vcf1).gz
bcftools index /cromwell_root/tmp/bcftools/$(basename $pon_vcf1).gz


bcftools view ${pon_vcf2} -Oz -o /cromwell_root/tmp/bcftools/$(basename $pon_vcf2).gz bcftools view ${pon_vcf2} -Oz -o /cromwell_root/tmp/bcftools/$(basename $pon_vcf2).gz
bcftools index /cromwell_root/tmp/bcftools/$(basename $pon_vcf2).gz


bcftools view ${pon_vcf3} -Oz -o /cromwell_root/tmp/bcftools/$(basename $pon_vcf3).gz bcftools view ${pon_vcf3} -Oz -o /cromwell_root/tmp/bcftools/$(basename $pon_vcf3).gz
bcftools index /cromwell_root/tmp/bcftools/$(basename $pon_vcf3).gz


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

Loading…
Cancel
Save