Question: How should i determine outFilterMultimapNmax in rna seq analysis, P.s im working on pseudogenes.
gravatar for sabaghianamir70
13 months ago by
sabaghianamir7010 wrote:

how do i know what the outFilterMultimapNmax must be for my analysis? consider that,i am working on pseudogenes.

rna-seq star pseudogenes • 568 views
ADD COMMENTlink modified 13 months ago by h.mon29k • written 13 months ago by sabaghianamir7010

You forgot to mention that this parameter is used with STAR (I have added the tag) and that your previous thread is here: Gene and pseudogenes

ADD REPLYlink written 13 months ago by Kevin Blighe60k
gravatar for h.mon
13 months ago by
h.mon29k wrote:

First of all, if you are interested in multi-mappers quantification, you do not want to use htseq-count, you want to use RSEM, Salmon or kallisto. Although your current post doesn't say you are using htseq-count, your previous post does.

According to the post Multimapping Reads and Pseudogenes, STAR has fairly stringent parameters and one mismatch between the best match and second-best match is enough to exclude the second-best mapping position, so in general, I would expect for most pseudogenes it is not necessary to alter --outFilterMultimapNmax, unless the duplication is fairly recent and / or gene and pseudogene are identical over a large stretch (longer than your sequencing read length).

Based on your past thread, it seems you are interested in a particular gene / pseudogene pair, as opposed to interested in transcription expression of all pseudogenes in the genome. If that is the case, you have to find out how similar PER3 and PER4 (and potentially PER1 and PER2) are. You can also test different --outFilterMultimapNmax settings and visualize the regions of interest with IGV to see how mapping there is affected.

ADD COMMENTlink written 13 months ago by h.mon29k

im working on rna sequencing project, and i want to find transcripts of pseudogenes with very low expression rate, but im tried that before and i get zero counts,(using Star, FeatureCount not htseq and limma for the DEG), i want to know what options i must set, to get some count from it.(some one told me for that i must use --outFilterMultimapNmax), for example Per3, it cant be 0, its an active gene, and has oscillation through the time and has very low expression rate but not Zero.

ADD REPLYlink written 13 months ago by sabaghianamir7010
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: 1336 users visited in the last hour