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

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.

rna-seq featurecounts • 487 views
ADD COMMENTlink written 8 months ago by dec986110
3
gravatar for genomax
8 months ago by
genomax55k
United States
genomax55k 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 8 months ago • written 8 months ago by genomax55k

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 8 months ago by dec986110
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.3.0
Traffic: 1723 users visited in the last hour