search

VCCRI / Sierra

Discover differential transcript usage from polyA-captured single cell RNA-seq data
GNU General Public License v3.0
49 stars 17 forks source link

Question on counting overlapping peaks #72

Open HenriettaHolze opened 1 month ago

HenriettaHolze commented 1 month ago

Hi, I'm looking at a transcript with a shorter an longer version of an exon.
I was wondering how you count reads for those peaks?
Like reads that map partially to the shorter exon but fully to the longer, or reads that map to the longer but the junction shows they belong to the shorter? I supply my own peak list using the UTR and CDS regions from my gtf file at the moment. Thanks!

rj-patrick commented 1 month ago

Hi @HenriettaHolze, The peaks are counted individually - in the end we test for the relative usage so we're not attempting to perform transcript reconstruction. So in your case reads overlapping the short and long exons will be counted for both.