Compute SNP pileup at reference positions in one or more input bam files. Output is ready for the R package facets