smithlabcode / ribotricer

A tool for accurately detecting actively translating ORFs from Ribo-seq data
http://doi.org/djv4
GNU General Public License v3.0
28 stars 8 forks source link

Ribotricer Output #141

Closed bshim181 closed 11 months ago

bshim181 commented 11 months ago
Screenshot 2023-07-21 at 3 05 00 PM Screenshot 2023-07-21 at 3 05 14 PM

Hello, This is my first time using ribotricer to interpret the ribosome sequencing dataset. I was kind of confused by what the proportion meant in the metagene profile graph. is the proportion of reads that are 29nt and 30nt long / total number of reads that in the bam file that have aligned to the translating ORFs?

Also, I was only able to find metagene profile for reads of these two lengths (29 nt and 30 nt). In the RibORF characterization,

Screenshot 2023-07-21 at 3 08 52 PM

I was able to notice level of periodicity in reads of other lengths as well. I was wondering what might have happened to reads of those lengths.

Screenshot 2023-07-21 at 3 09 47 PM

This is the bam file summary.

saketkc commented 11 months ago

We should clarify the meaning of the percentage in the documentation. But to clarify, ribotricer only plots metagene profile of read lengths that pass the phase score cutoff. In this case you just have two read lengths (29,30) that meet this criterion. Hope this helps. Feel free to reopen if there are any follow up questions.