#!/bin/bash -Ceuo pipefail
# Matplotlib cache dir must be writable; default location may not be
export MPLCONFIGDIR=/tmp/matplotlib

all_qc \
    -bam Sig_18_tissue.md.bam \
    -reference Homo_sapiens_assembly38.fasta \
    -beds xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed xgen-exome-hyb-panel-v2_AND_altera_v3_probes_short_hg38_AND_cdx_spike_custom30off4x_hg38_AND_foresight_clarity_footprint.bed altera_v3_targets_postQC_hg38.bed altera_v3_targets_coding_postQC_hg38.bed altera_v3_targets_noncoding_postQC_hg38.bed xgen-exome-hyb-panel-v2-targets-hg38_minus_altera_v3_targets_postQC_hg38.bed FINAL_hotspot_list_07302025_hg38_sorted.bed foresight_clarity_targets_hg38.bed \
    -dbsnp_bed dbsnp138_hg38_common.bed \
    -dbsnp_vcf Homo_sapiens_assembly38.dbsnp138.vcf.gz \
    -contamination_bed dbsnp_hg38_contam_sites.bed \
    -output_folder Sig_18_tissue \
    -threads 31 \


merge_qc -folders Sig_18_tissue/* -output Sig_18_tissue/Sig_18_tissue_full_basic_stats.tsv

head -n1 Sig_18_tissue/Sig_18_tissue_full_basic_stats.tsv > Sig_18_tissue/Sig_18_tissue_basic_stats.tsv
grep -e "$(basename xgen-exome-hyb-panel-v2-targets-hg38_AND_altera_v3_targets_postQC_hg38_AND_foresight_clarity.bed)" Sig_18_tissue/Sig_18_tissue_full_basic_stats.tsv >> Sig_18_tissue/Sig_18_tissue_basic_stats.tsv

cat <<-END_VERSIONS > versions.yml
"DAQ:BAM_POSTDEDUP_QC:SEQTOOL":
    seqtool: 2.0
END_VERSIONS