Counting reads for gene transcripts with overlapping loci
0
5
Entering edit mode
8.8 years ago
Ruet ▴ 70

I am working with the human genome and have been allowing single mapping during alignment with Tophat. When I count reads that map to a particular loci, if a read falls on two gene transcripts, I do not count the ambiguous read. The transcripts I am interested in are snRNAs and snoRNAs but is seems like the transcripts of interest are falling under longer transcripts. Therefore, I count reads from the smaller snRNA transcripts as well as some reads from the longer transcript. I was wondering how to counts reads for overlapping transcripts especially when the shorter gene transcripts falls entirely under a longer gene transcript?

RNA-Seq gene R • 3.3k views
ADD COMMENT
0
Entering edit mode

Is your RNA-seq protocol stranded? Are you using this information in your read counting?

ADD REPLY
2
Entering edit mode

Hi Mark,

Yes, my RNA-seq protocol is stranded and I am using the information for read counting in my analysis. I am using the "union" model in the figure below, and I am not counting ambiguous reads.

image

Thanks!

ADD REPLY

Login before adding your answer.

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