#!/bin/bash -Ceuo pipefail
diff_dv_sentieon.py \
--dv-vcf germline_control_2.dv_rescue.vcf.gz \
--sentieon-vcf germline_control_2.haplotyper.unfiltered.vcf.gz \
--output-vcf germline_control_2.dv_only_sites.vcf.gz
tabix -p vcf germline_control_2.dv_only_sites.vcf.gz
cat <<-END_VERSIONS > versions.yml
"DAQ:CONTROL_VARIANT_CALLING:CONTROL_GERMLINE_VC:BAM_VARIANT_CALLING_SENTIEON_HAPLOTYPER_RF:DIFF_DV_SENTIEON":
python: $(python --version | sed 's/Python //g')
cyvcf2: $(python -c "import cyvcf2; print(cyvcf2.__version__)")
END_VERSIONS