Question: featureCounts not recognizing paired-end alignment file.. why?
gravatar for dec986
5 months ago by
United States
dec986100 wrote:


I am running featureCounts on publicly available paired-end data. After alignment, I run featureCounts, but there are very few successfully assigned reads (only 25.2% successfully assigned).

I see

|              Paired-end : no

Why is featureCounts not seeing that this data is paired-end? Is this why I am getting so few reads successfully assigned?


rna-seq featurecounts • 326 views
ADD COMMENTlink written 5 months ago by dec986100
gravatar for genomax
5 months ago by
United States
genomax49k wrote:

Are you adding -p option to signify that you want fragments counted (not reads, in case of PE)? Just to be sure the alignment was done using PE reads at the same time, correct?

Also remember that multi-mapping reads are not counted by default so that may explain the low % of assignment.

ADD COMMENTlink modified 5 months ago • written 5 months ago by genomax49k

I was confused by the menu option,

-p                  If specified, fragments (or templates) will be counted
                      instead of reads. This option is only applicable for
                      paired-end reads.

this made it sound like it automatically detects whether or not it is paired end.

ADD REPLYlink written 5 months ago by dec986100
Please log in to add an answer.


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