powered by
Gets pileup from bam over positions.
bamToPileup(bam, fasta, positions, index, BQoffset = 33)