瀏覽代碼

first commit

master
YaqingLiu 5 年之前
父節點
當前提交
c9738d65d8
共有 1 個文件被更改,包括 0 次插入8 次删除
  1. +0
    -8
      tasks/bcftools.wdl

+ 0
- 8
tasks/bcftools.wdl 查看文件

if [ ${PONmode} ]; then if [ ${PONmode} ]; then
mkdir -p /cromwell_root/tmp/bcftools mkdir -p /cromwell_root/tmp/bcftools
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 view ${pon_vcf2} -Oz -o /cromwell_root/tmp/bcftools/$(basename ${pon_vcf1}.gz) bcftools view ${pon_vcf2} -Oz -o /cromwell_root/tmp/bcftools/$(basename ${pon_vcf1}.gz)
bcftools view ${pon_vcf3} -Oz -o /cromwell_root/tmp/bcftools/$(basename ${pon_vcf1}.gz) bcftools view ${pon_vcf3} -Oz -o /cromwell_root/tmp/bcftools/$(basename ${pon_vcf1}.gz)
cd /cromwell_root/tmp/bcftools cd /cromwell_root/tmp/bcftools
#bcftools index *.vcf.gz
#bcftools index FUSCCTNBC330_pon.vcf.gz
#bcftools index FUSCCTNBC467_pon.vcf.gz

for i in ${sep=" " *.vcf.gz} for i in ${sep=" " *.vcf.gz}
do do
bcftools index $i bcftools index $i
done done

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 *.vcf.gz | bcftools plugin fill-AN-AC | bcftools filter -i 'SUM(AC)>1' > panel_of_normal.vcf
else else
touch panel_of_normal.vcf touch panel_of_normal.vcf
fi fi

>>> >>>


runtime { runtime {

Loading…
取消
儲存