Hi. This is my first time using featureCounts and I need your help (🚨). I have single-end reads and used HTSeq-count several times with the -m (mode) option specified to --intersection-nonempty https://htseq.readthedocs.io/en/release_0.11.1/count.html. However, I am now trying to run featureCounts and not 100 percent sure that which option of in the featureCounts is equal to -m intersection-nonempty in HTSeq-count.
Any help will be appreciated.
Thanks.
•
written
by
Apex92 •
1 answer
featureCounts does not implement all the features of htseq-counts. Conceptually the tool works differently and uses different internal representations.
I don't think that the intersection-nonempty mode is implemented.
you should think about why you want that to be set and what it actually means in the larger context.
why would you want that parameter to be set? is it because you want one feature to be fully covered, or completely included? etc perhaps other parameters of featureCounts do an even better job achieving that goal.
in my opinion the example in the htseq count documentation does not fully explain and define what intersection-nonempty does
Log in to answer this question.
-
htseq-count : dosen't show options in bash ln ubuntu
written by Foad •
-
htseq-count without GTF?
written by blur
-
How many of the reads end up being counted by HTSeq count
written by Apex92 •
-
HT-seq-count default parameter to "score"secondary alignments contradicts its documentation that it…
written by Zee_S
-
Best practices for counting genes with nanopore RNAseq data
written by Florian Bernard
-
Which htseq-count mode to use for handling overlaps during RNA-Seq quantitation
written by skhan •
-
Htseq-count Why Do I have more nofeature readswhen i'm changing nonunique option ?
written by bjmnpnd •
-
Is it necessary to name sort single end reads before htseq count
written by natsterbug •
-
Help with a reading count command using HTSeq
written by flavobacteria •
-
Using HTSeq for paired-end files
written by Whoknows
Are you sure you want HTSeq-count? There are better tools for evaluating transcript levels nowadays, if this is what you are indeed doing.
I want to use featureCounts but not sure how to specify the mapping option as I did in HTSeq count (
-m intersection-nonempty)