There is a function called geneBody_coverage.py in the Python tool RseQC, which determines whether the coverage of Reads is uniform along the transcripts, o if there is a 3' or 5' bias.
Now I want to apply this tool to my research question: I want analyse my coverage within the 5'UTR. However, due to the fact that the 5'UTR of each transcript is different, I have to normalize it first. So in the end I want to have a plot with "Percentile of 5'UTR" on the x-axis and "Read number" on the y-axis.
Does anyone of you know a convient alternative for performing this analysis in R? I already have my bam file in R (imported via "readGAlignments")+ I searched all 5'UTRs with "fiveUTRsByTranscript".
Thanks for any ideas