#!/bin/bash -euo pipefail
filter_hotspot_vcf.py \
--input-vcf 1136_5WB-T1-TDNA-1_A23WJ5GLT4_8_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_8_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_8_vs_1136_5WB-N1-BDNA-1_B23KFG5LT4_1.tnseq.hotspot.hmnz.vcf', '1136_5WB-T1-TDNA-1_A23WJ5GLT4_8_vs_1136_5WB-N1-BDNA-1_B23KFG5LT4_1.tnseq.hotspot.hmnz.vcf.gz', force=True); pysam.tabix_index('1136_5WB-T1-TDNA-1_A23WJ5GLT4_8_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