gravatar for Ready2Rapture

3 hours ago by

Hi all,

So I'm running alternative splicing (STAR > RSEM) on bulk RNAseq with mate read lengths supposed to be 151.

I don't have much experience dealing with isoform analysis & read lengths this long and noticed that while the majority read lengths are 151 (for example 1 million) almost half as many are 150 (for example 550,000).

My intuition is to still set sjdbOver hang as 150 since I read length distributions vary around the targeted length, but wanted to double check.

Source link