#!/bin/bash -euo pipefail
filter_hotspot_vcf.py \
    --input-vcf 1136_5WB-T1-TDNA-1_A23WJ5GLT4_9_vs_1136_5WB-N1-BDNA-1_B23KFG5LT4_1.tnseq.hmnz.vcf.gz \
    --hotspot-bed FINAL_hotspot_list_07302025_hg38_sorted_w_anno_mods.bed \
    --output-vcf 1136_5WB-T1-TDNA-1_A23WJ5GLT4_9_vs_1136_5WB-N1-BDNA-1_B23KFG5LT4_1.tnseq.hotspot.hmnz.vcf

python3 -c "import pysam; pysam.tabix_compress('1136_5WB-T1-TDNA-1_A23WJ5GLT4_9_vs_1136_5WB-N1-BDNA-1_B23KFG5LT4_1.tnseq.hotspot.hmnz.vcf', '1136_5WB-T1-TDNA-1_A23WJ5GLT4_9_vs_1136_5WB-N1-BDNA-1_B23KFG5LT4_1.tnseq.hotspot.hmnz.vcf.gz', force=True); pysam.tabix_index('1136_5WB-T1-TDNA-1_A23WJ5GLT4_9_vs_1136_5WB-N1-BDNA-1_B23KFG5LT4_1.tnseq.hotspot.hmnz.vcf.gz', preset='vcf', force=True)"

cat <<-END_VERSIONS > versions.yml
"NFCORE_HARMONIZE:HARMONIZE:FILTER_HOTSPOT_VCF":
    python: $(python --version | sed 's/Python //')
    filter_hotspot_vcf: "1.0.0"
    pysam: $(python3 -c "import pysam; print(pysam.__version__)")
END_VERSIONS