Hi @Wesserg - After reading your recent publication in Computational Biology and Chemistry, I'm interested in applying SRSFseq to RNA-seq data where I'm expecting clear (but relatively small) differences between the shape profiles over genes. For example, the attached IGV screengrab shows two examples of control vs. treatment where there is a clear reduced coverage in the 3' UTR at the ACTB gene. In the lower example, there is also a less prominent reduction in coverage towards the end of the largest exon. I have 4 replicates for all conditions and these patterns are consistent. Overall, I'm expecting the majority of such changes to be in 3' UTRs.

Do you think SFSRseq would be an appropriate method to identify such differences, and if so what's the best approach to do so? From what I gather, SRSFseq is not available as a stable R package. Do you have any plans to add this to Bioconductor for instance? In the meantime, I'll work from the code available in this repository

fig1.pdf