search

darelab2014 / Dinopore

DInoPORE: Direct detection of INOsines in native RNA with nanoPORE sequencing
16 stars 3 forks source link

confused with event_length * 3012 #11

Open baishengjun opened 7 months ago

baishengjun commented 7 months ago

Hi @darelab2014, Firstly thanks for developing such an awesome tool. I have been using nanopolish eventalign for most of my work with Nanopore data. Recently I read Dinopore source code from s2.Combine_raw_nanopolish.R script, I was confused that why event_length multiped by 3012? In my understand, the event_length is dwelling time of 5-mer stay in pore. Can you give me some explained? 

  dat=dat[dat$strand=="t" & dat$reference_kmer != "NNNNN" & dat$event_stdv<50,]
  dat[,"count"]=round(3012*dat$event_length)

Thanks a lot, Bai