Question: Rsamtool or Rbamtools function to view the number of reads between a range
0
gravatar for startup_biostar
3.2 years ago by
startup_biostar10 wrote:

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

rsamtools rna-seq rbamtools R • 939 views
ADD COMMENTlink modified 3.2 years ago by Devon Ryan91k • written 3.2 years ago by startup_biostar10
2
gravatar for Devon Ryan
3.2 years ago by
Devon Ryan91k
Freiburg, Germany
Devon Ryan91k wrote:

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

ADD COMMENTlink written 3.2 years ago by Devon Ryan91k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 933 users visited in the last hour