if true, use MD tag together with CIGAR to recover reference bases
position from which to start
position before which to stop
That is, the range of positions is half-open interval
[max(start_from, first mapped read start position),
min(end_at, last mapped end position among all reads))
if true, skip sites with zero coverage
Creates a pileup range from a range of reads. Note that all reads must be aligned to the same reference.
See PileupColumn documentation for description of range elements. Note also that you can't use std.array.array() function on pileup because it won't make deep copies of underlying data structure. (One might argue that in this case it would be better to use opApply, but typically one would use std.algorithm.map on pileup columns to obtain some numeric characteristics.)