|
|
|
|
|
|
|
|
# bcftools norm -m -both ${vcf} | bcftools norm -f ${ref_dir}/${fasta} -Ov -o ${basename}.norm.vcf |
|
|
# bcftools norm -m -both ${vcf} | bcftools norm -f ${ref_dir}/${fasta} -Ov -o ${basename}.norm.vcf |
|
|
# Split multiallelic sites |
|
|
# Split multiallelic sites |
|
|
bcftools norm -m -both ${vcf} -o ${basename}.norm.vcf |
|
|
bcftools norm -m -both ${vcf} -o ${basename}.norm.vcf |
|
|
|
|
|
bcftools view -O z -f 'PASS' ${basename}.norm.vcf -o ${basename}.PASS.vcf.gz |
|
|
>>> |
|
|
>>> |
|
|
|
|
|
|
|
|
runtime { |
|
|
runtime { |
|
|
|
|
|
|
|
|
} |
|
|
} |
|
|
|
|
|
|
|
|
output { |
|
|
output { |
|
|
File norm_vcf = "${basename}.norm.vcf" |
|
|
|
|
|
|
|
|
File norm_vcf = "${basename}.PASS.vcf.gz" |
|
|
} |
|
|
} |
|
|
} |
|
|
} |