fiberseq/fibertools-rs

Does m6a prediction require fiber-seq preperation

Closed this issue · 3 comments

Hey Mitchell,

Could you clarify if fibertools is able to predict m6a from normally prepared HiFi reads (with full kinetic data) or if they require the special fiber-seq preparation before sequencing? I ran ft on some normally prepared data and the output bam had all the new m6a tags, but the new preprint seems to suggest only the latter is true. The IGV looked similar to yours here, but I'm not sure if mine are all nonsense.

best,
Alex

Hi Alex. Thanks for using the tool!

If you wish to use it as a chromatin accessibility assay then you do need to use the Fiber-seq protocol before sequencing. However if you want to detect endogenous m6a then what you have is sufficient (note, eukaryotes to my knowledge don't have endogenous DNA m6a.)

What is your use case?

Cheers,
Mitchell

I had several samples of cattle sequenced on a Revio with the full kinetics (but no special preparation) so was thinking I could potentially get an extra layer of data for free out if it. It doesn't sound like that is the case though.

Yeah sorry. The m6a calls you got are almost certainly just false positives.