Is there an R package that can calculate alignment statistics from bam/sam files? By alignment statistics I mean number of reads, number of reads uniquely mapped, number of reads mapped to multiple locations.
I found rsubread package can count number of reads that can be mapped. However, it doesn't differentiate between unique alignment and multiple alignment.
Thanks in advance!