add formats.sam.pairs() to plot insert size distribution given BAM output