exchange samtools mpileup with bcftools mpileup [OPTIONS] -f ref.fa in.bam [in2.bam […]] https://samtools.github.io/bcftools/bcftools.html#common_options