featureCounts not recognizing paired-end alignment file.. why?
1
0
Entering edit mode
6.3 years ago
dec986 ▴ 370

hello,

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?

thanks.

featurecounts RNA-Seq • 3.0k views
ADD COMMENT
3
Entering edit mode
6.3 years ago
GenoMax 141k

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 COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 2601 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6