GoekeLab / m6anet

Detection of m6A from direct RNA-Seq data
https://m6anet.readthedocs.io/
MIT License
104 stars 19 forks source link

mod_ratio availability in data.indiv_proba.csv #107

Closed sunsetyerin closed 1 year ago

sunsetyerin commented 1 year ago

I realized that, unlike data.site_proba.csv, data.indiv_proba.csv doesn't include mod_ratio result. I wonder if mod_ratio can be supported in the data.site_proba.csv result. How can I get m6A methylation level from data.indiv_proba.csv?

chrishendra93 commented 1 year ago

hi @sunsetyerin ,

data.indiv_proba.csv does not contain mod ratio since each row describes read level probability score. So if you need an estimate for how likely each read to be modified, the probability_modified column in data.indiv_proba.csv will tell you that. The mod ratio in data.site_proba.csv is derived by thresholding the read probability score and calculating the number of reads that are modified from each site