#!/bin/bash -Ceuo pipefail
sentieon driver -r Homo_sapiens_assembly38.fasta \
 \
--algo TNfilter \
--tumor_sample HCC1395_tumor --normal_sample reference-NA12878 --max_event_count 5 \
-v positive_somatic_control_1.tnseq.vcf.gz \
 --contamination positive_somatic_control_1.tnseq.contamination_data.tsv  \
--tumor_segments positive_somatic_control_1.tnseq.segments \
--orientation_priors positive_somatic_control_1.tnseq.orientation_data.tsv \
positive_somatic_control_1.tnseq.filtered.vcf.gz

cat <<-END_VERSIONS > versions.yml
"DAQ:CONTROL_VARIANT_CALLING:CONTROL_SOMATIC_VC:BAM_VARIANT_CALLING_SOMATIC_TNSEQ:SENTIEON_TNFILTER":
    sentieon: $(echo $(sentieon driver --version 2>&1) | sed -e "s/sentieon-genomics-//g")
END_VERSIONS