what is Equivalent R function for the command line below?
samtools view yourBAMfile.bam chr1:567876-568100 | wc -l
How should I approach.
I got a
GRanges object with seqnames, start and end. I would like to count the number of reads between each range in the
GRange object in a
BAM file in R