Search
Question: Counting ambiguously mapped reads per feature
0
gravatar for rubi
16 months ago by
rubi70
rubi70 wrote:

Hi,

I have a bam file in which ~50% of the reads map ambiguously (since the library is from small RNAs). For each feature (e.g., gene) in my GTF file I'd like to add a count of 1 if a read from that bam file maps to that feature uniquely and 1/n if a read maps to n locations (that feature being one of them).

Is there any Bioconductor code that can achieve that?

 

ADD COMMENTlink modified 16 months ago by Gordon Smyth32k • written 16 months ago by rubi70
3
gravatar for Gordon Smyth
16 months ago by
Gordon Smyth32k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth32k wrote:

Rsubread::featureCounts() with countMultiMappingReads=TRUE and fraction=TRUE does exactly that.

ADD COMMENTlink modified 16 months ago • written 16 months ago by Gordon Smyth32k
3
gravatar for Steve Lianoglou
16 months ago by
Genentech
Steve Lianoglou12k wrote:

You might try using tools like salmon, kallisto, or RSEM, which provide a more principled (and community-vetted) approach to solving this exact issue.

ADD COMMENTlink written 16 months ago by Steve Lianoglou12k
2

Here are my Bioc2016 slides on using these quantifiers with Bioconductor software packages:

http://mikelove.github.io/bioc2016

In particular, see slide 4 for an idea how you can do better than assigning 1/n to all of the potential alignment locations.

ADD REPLYlink written 16 months ago by Michael Love15k
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.2.0
Traffic: 270 users visited in the last hour