#!/bin/bash -Ceuo pipefail
echo "Sig_18_Blood" > sample_name.list

bcftools \
    mpileup \
    --fasta-ref Homo_sapiens_assembly38.fasta \
     \
    Sig_18_Blood.recalibrated.bam \
    -T SNP_GRCh38_hg38_wChr.bed \
     \
    | bcftools call --output-type v --no-version --ploidy 1 -c \
    | bcftools reheader --samples sample_name.list \
    | bcftools view --output-file Sig_18_Blood.ngscheckmate.vcf.gz --output-type z --no-version



tabix -p vcf -f Sig_18_Blood.ngscheckmate.vcf.gz

bcftools stats Sig_18_Blood.ngscheckmate.vcf.gz > Sig_18_Blood.ngscheckmate.bcftools_stats.txt

cat <<-END_VERSIONS > versions.yml
"NFCORE_SAREK:SAREK:CRAM_SAMPLE_INTEGRITY:BAM_NGSCHECKMATE:BCFTOOLS_MPILEUP":
    bcftools: $(bcftools --version 2>&1 | head -n1 | sed 's/^.*bcftools //; s/ .*$//')
END_VERSIONS