Rsamtool or Rbamtools function to view the number of reads between a range
1
0
Entering edit mode
8.6 years ago

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

RNA-Seq Rsamtools rbamtools R • 1.8k views
ADD COMMENT
2
Entering edit mode
8.6 years ago

rbamtools directly provides a count function (bamCount), so just go with that.

ADD COMMENT

Login before adding your answer.

Traffic: 2184 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6