alpaca call -t 1 --fdr 0.05 'HG00100' < all.bcf > calls.2.bcf