#!/bin/bash -Ceuo pipefail
sentieon driver \
    -t 12 \
    -r Homo_sapiens_assembly38.fasta \
     \
    -i HCC1395_BL.recalibrated.bam -i HCC1395_tumor.recalibrated.bam \
    --interval xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_chr21_padded100bp.bed \
    --algo TNhaplotyper2 \
    --tumor_sample HCC1395_tumor --normal_sample HCC1395_BL \
    --call_germline_sites \
    --germline_vcf af-only-gnomad.hg38.vcf.gz \
    --pon pon_tnseq_42_curated_v4.vcf.gz \
    HCC1395_tumor_vs_HCC1395_BL.tnseq.vcf.gz \
    --algo OrientationBias --tumor_sample HCC1395_tumor HCC1395_tumor_vs_HCC1395_BL.tnseq.orientation_data.tsv \
    --algo ContaminationModel --tumor_sample HCC1395_tumor --normal_sample HCC1395_BL --vcf af-only-gnomad.hg38.vcf.gz --tumor_segments HCC1395_tumor_vs_HCC1395_BL.tnseq.segments HCC1395_tumor_vs_HCC1395_BL.tnseq.contamination_data.tsv

cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:BAM_VARIANT_CALLING_SOMATIC_ALL:BAM_VARIANT_CALLING_SOMATIC_TNSEQ:SENTIEON_TNHAPLOTYPER2":
    sentieon: $(echo $(sentieon driver --version 2>&1) | sed -e "s/sentieon-genomics-//g")
END_VERSIONS