Hi! myfriends,
I'm doing read count of my RNA-seq data using featurecounts, I found a pretty low assigned rate which ranges 20~30%, like below.
FeatureCounts result
I noticed that the total is not accord with my FastQC and STAR' report:
So, I used samtools flagstat to check:
I'm in doubt that whether featureCounts and samtools recounted multi-mapping reads? and if so, I wish your kindly suggestions, Thank you!
1 answer
Those numbers don't look totally out of line to me. You have some reads aligning uniquely to a single feature, which featureCounts will count, some reads aligning uniquely, but not to features, and some reads aligning to multiple places.
I've run featureCounts on bams with secondary reads, and bams where I've filtered them away, and I get the same number of features (though the number of Unassigned_MultiMapping did change)
Log in to answer this question.
-
A question on MEME results
written by Hughie •
-
software that can draw this figure
written by yueli7
-
Need help in selecting the correct output file from miRDeep2?
written by bioinforesearchquestions
-
Anyone knows how to create a figure like this?
written by Anny •
-
featureCounts results discordance with library information
written by Hughie •
-
KaryoplotR, display 2 different chromosome scales on the same plotKaryotype
written by Bastien Hervé
-
Need help in selecting the correct output file from miRDeep2?
written by bioinforesearchquestions
-
Tutorial: CorLevelPlot - Visualise correlation results, e.g., clinical parameter correlations
written by Kevin Blighe
-
E-value results from hmmsearch are not accurate. [HMMER]
written by danzanzu •
-
featurecounts of Rsubread strand specific orientation options
written by ofonov •
I've previously seen examples were a lot of RNA-seq reads map to other parts of the genome (introns, intergenic) - could you test whether that is the case?
If " reads map to other parts of the genome (introns, intergenic)" should I remove the sample from downstream analysis?
featureCountsdoes NOT count multi-mapped reads by default (which you don't want to do anyway). You can force it to do so using-Moption. Check to see if you have rRNA contamination in your data if you see a large amount of multi-mapping.