differences between featureCounts and htseq-count
1
0
Entering edit mode
18 months ago
golany89 • 0

Hi, I have RNA-seq data (from ribodepleted samples) and i'm trying to map my reads to features using either htseq or featureCounts. The input for both is a bam file (reads were aligned using STAR) and a GTF file from ENSEMBLE. Both assign reads to exons by default.

FeatureCounts reported that ~33% of my reads were assigned to features, however when I'm trying to check how many reads were assigned using htseq, I saw that only 2% of my reads were assigned (There is no report in htseq so I summerized the second column in the results file, without the last 5 rows that contain the "no-feature" and "ambiguous" etc).

What is causing this difference between the two algorithms? Thank you very much!

RNA-Seq • 1.0k views
ADD COMMENT
2
Entering edit mode

The differences should be very minor between the two tools. If you post the exact commands used, it will be easier to see what causes the differences.

ADD REPLY
2
Entering edit mode

Is it possible you're using the wrong stranded argument in htseq? The vast majority of stranded libraries (e.g. using the Illumina kit) should use --stranded=reverse NOT --stranded=yes

ADD REPLY
0
Entering edit mode

Thank you very much, this solved the problem.

ADD REPLY
1
Entering edit mode
ADD COMMENT

Login before adding your answer.

Traffic: 1478 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