Entering edit mode
8.6 years ago
startup_biostar
▴
20
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