search

fiberseq / fibertools-rs

Tools for fiberseq data written in rust.
https://fiberseq.github.io/fibertools/fibertools.html
42 stars 5 forks source link

how to state the genome 6mA marks #58

Closed cell-c closed 4 days ago

cell-c commented 1 month ago

Thanks for providing fibertools-rs. It is an excellent tool for detecting 6mA marks using PacBio HiFi sequencing.

I don't find this function in fibertools-rs about the statistics of the genome 6mA marks. How is this function of the software found? Thanks.

mrvollger commented 1 month ago

Can you elaborate on what you are looking for?

mrvollger commented 4 days ago

We have added a command called ft qc that should resolve this.

cell-c commented 4 days ago

Hi, I am sorry for missing your response. Accturaly, I predicted the 6mA and 5mc marks using this tool. I want to obtain the methylation statistics in the reference genome. Currently I map the methylation bam result to the reference genome using pbmm2 and analyse the methylation information in the reference genome based on the coverage of methylation sites. Do you have any suggestions about the statistics? Thank you.

StephanieBohaczuk commented 3 days ago

Maybe you're looking for this? https://github.com/fiberseq/FIRE