#!/bin/bash -Ceuo pipefail
samtools \
    stats \
    --threads 1 \
    --reference Homo_sapiens_assembly38.fasta \
    1029_11H-N1-BDNA-1_B23KGCJLT4_1.recalibrated.bam \
    > 1029_11H-N1-BDNA-1_B23KGCJLT4_1.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