#!/bin/bash -Ceuo pipefail
samtools \
    stats \
    --threads 4 \
    --reference Homo_sapiens_assembly38.fasta \
     \
    Sig_18_tissue.md.bam \
    > Sig_18_tissue.stats

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