SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
featureCounts segmentation fault batel Bioinformatics 24 12-30-2017 11:24 PM
featureCounts: a universal read summarization program shi Bioinformatics 149 05-02-2017 10:05 AM
featureCounts error heso Bioinformatics 8 11-14-2014 01:14 PM
Bug in featureCounts? brinda Bioinformatics 1 01-23-2014 07:35 PM
Bug in featureCounts? brinda Bioinformatics 1 01-19-2014 05:22 PM

Reply
 
Thread Tools
Old 03-16-2015, 03:51 AM   #1
SylvainL
Senior Member
 
Location: Geneva

Join Date: Feb 2012
Posts: 175
Default featureCounts

Hi all,

I am using (for the first time) featureCounts to get gene level counts in a RNAseq experiment... It is mouse RNAseq and polyA RNA were extracted...

A quick look at the results shows up to 15% of reads which are unassigned because of no overlap with any feature.

Do you normally have the same stats?

s.
SylvainL is offline   Reply With Quote
Old 03-16-2015, 04:20 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

I just checked one dataset (mouse liver), where ~5% had no feature, ~3% were ambiguously mapped (a non-stranded protocol), and ~19% were not uniquely mapped. ~75% were counted if we exclude unmapped reads. This is similar in a mouse hippocampus dataset that I just checked. Obviously if you include unmapped reads then these percentages would go down a bit.

So, that unassigned percentage seems a bit high, but perhaps there's a good biological reason for that (e.g., in tumors or sperm).
dpryan is offline   Reply With Quote
Old 03-16-2015, 04:31 AM   #3
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,814
Default

Post #29. Response from Wei Shi (author of featureCounts): http://seqanswers.com/forums/showthr...t=49506&page=2
GenoMax is offline   Reply With Quote
Old 03-16-2015, 04:32 AM   #4
SylvainL
Senior Member
 
Location: Geneva

Join Date: Feb 2012
Posts: 175
Default

Hi dpryan,

I already removed the non uniquely mapped so yes, it's a bit lower if we consider all the mapped reads... and yes, it's sperm...

So I should not worry :9

Thanks a lot
SylvainL is offline   Reply With Quote
Old 03-16-2015, 04:35 AM   #5
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

Yeah, sperm will inherently be different, since you have depletion of mRNAs. My sperm datasets aren't polyA enriched, they're ribo-depleted, so I don't have any numbers for comparison.
dpryan is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:43 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO