#!/bin/bash -Ceuo pipefail
sentieon driver \
\
-r Homo_sapiens_assembly38.fasta \
-t 12 \
--interval germline_rescue_targets.bed \
-i HCC1395_BL.recalibrated.bam \
\
--algo Haplotyper --given germline_control_1.dv_only_sites.vcf.gz --emit_mode variant germline_control_1.forcecall.unfiltered.vcf.gz \
cat <<-END_VERSIONS > versions.yml
"DAQ:CONTROL_VARIANT_CALLING:CONTROL_GERMLINE_VC:BAM_VARIANT_CALLING_SENTIEON_HAPLOTYPER_RF:SENTIEON_FORCECALL":
sentieon: $(echo $(sentieon driver --version 2>&1) | sed -e "s/sentieon-genomics-//g")
END_VERSIONS