will the bam files mapped from Ensembl index be ok to do featurecount with mm10
0
0
Entering edit mode
3.4 years ago
Kai_Qi ▴ 130

Hi:

Today, I suddenly realized that I have been using Emsembl index to do mapping while when do the feature counting I turn to FeatureCounts in Rsubread using the in-built mm10 because it does have GRCm38. Obviously, the annotations are different, but will that affect the accuracy of results?

Thank you very much,

RNA-Seq rna-seq alignment sequence • 808 views
ADD COMMENT
1
Entering edit mode

If chromosome names don't match then no, which is going to be the case here. Use Ensembl GTF to do the counting.

ADD REPLY
0
Entering edit mode

Thank you so much for clarifying. I notice the differences on the successfully assigned alignment:

Successfully assigned alignments : 31662019 (33.6%)

Previously I used the in-built annotation it is 60%-70%. Is this normal?

Thank you so much.

ADD REPLY
0
Entering edit mode

I find the reason. It is depends on what feature I am using to count. Sorry for being naive. I would like to keep the notes here.

ADD REPLY

Login before adding your answer.

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