#!/bin/bash

echo $@

BAM=$1
shift
REF=$1
shift
OUTVCF=$1

CMD="samtools mpileup -t DP,DV,DP4,DPR,SP -Auf $REF $BAM | bcftools call -m -f GQ > $OUTVCF"

echo "CMD: "$CMD
eval $CMD
