SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
featureCounts: a universal read summarization program shi Bioinformatics 149 05-02-2017 09:05 AM
Bug in featureCounts? brinda Bioinformatics 1 01-19-2014 04:22 PM
Possible tophat 2.0.6 bug? panagiotis RNA Sequencing 1 12-20-2012 07:50 AM
Is it a bug in BWA? SD2010Bioinfo Bioinformatics 4 10-11-2012 05:27 PM
is it bug in cufflinks? jay2008 Bioinformatics 0 10-10-2012 02:16 PM

Reply
 
Thread Tools
Old 01-23-2014, 04:55 PM   #1
brinda
Junior Member
 
Location: boston

Join Date: Jan 2014
Posts: 2
Default Bug in featureCounts?

We are using featureCounts for both and single-end and pair-end data. In the command we always have the '-p' flag on, assuming it will be ignored for single-end data. However, we notice that the function of '-M' is affected because of '-p' flag.

When we run single-end data without the -M flag, all the reads with NH > 1 are ignored, as expected. In the test data attached, we get a read count of 4. With -M flag, we get a count of 16, including reads with NH > 1, again, as expected.

However, when we run single-end data without the -M flag but with '-p' flag, the multi-mapping reads (NH > 1) are no longer ignored and we get a read count of 16.

The '-p' should have no effect on filtering multi-mapping reads. I'm attaching a sample sam and gtf file where this behavior may be replicated.
brinda is offline   Reply With Quote
Old 01-23-2014, 06:35 PM   #2
shi
Wei Shi
 
Location: Australia

Join Date: Feb 2010
Posts: 235
Default

As stated in the program usage, the '-p' option should only be used for the summarization of paired-end reads. So you should not use it for the single-end reads.

featureCounts tries to find the pair for each read when '-p' is specified. This causes a problem when single-end reads are provided and '-p' option was incorrectly specified.

In the next release (late Februray hopefully), featureCounts will automatically detect the read type (single-end or paired-end) and it will ignore '-p' option when single-end reads are provided.

Wei
shi is offline   Reply With Quote
Reply

Tags
featurecounts, rna-seq

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 08:48 PM.


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