Feature requests: read input as bam

[grenaud]

For BWA mem (or aln). Would it be possible to read the input as unaligned BAM? If we have a mix of merged reads as single-end reads and unmerged reads as pairs, we need to run the program twice and merge.