#!/bin/bash -Ceuo pipefail
samtools \
    stats \
    --threads 2 \
    --reference Homo_sapiens_assembly38.fasta \
    HCC1395_BL.recalibrated.bam \
    > HCC1395_BL.stats

cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:CRAM_SAMPLE_INTEGRITY:CRAM_QC_RECAL:SAMTOOLS_STATS":
    samtools: $(echo $(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*$//')
END_VERSIONS