Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Setting up input for epigenetic modification detection #324

Open
GianlucaMattei opened this issue Oct 3, 2022 · 1 comment
Open

Setting up input for epigenetic modification detection #324

GianlucaMattei opened this issue Oct 3, 2022 · 1 comment

Comments

@GianlucaMattei
Copy link

Hello,
I would like to run 5hmC and 5mC detection.

in this case should I use the dna_r9.4.1_450bps_modbases_5mc_hac.cfg or the dna_r9.4.1_450bps_modbases_dam-dcm-cpg_hac.cfg to detect all the methylated Cs (not only the in the GpGs context) ?

What is the difference between the two configs and is there any description about them anywhere?

About --ref-mods-all-motifs: is "N 0" correct (to retrieve also other types of modifications) or I should just use "C 0" ? (I was looking at this https://github.com/jkbonfield/hts-specs/blob/methylation/SAMtags.tex#L477)

@marcus1487
Copy link
Collaborator

The dna_r9.4.1_450bps_modbases_dam-dcm-cpg_hac.cfg model was an initial proof of concept model and is not recommended. The dna_r9.4.1_450bps_modbases_5mc_hac.cfg model is an improvement over this and is the best model for R9 data. Though the best results for all-context modified base detection can be achieved with kit14 and the Rerio all-context Remora model. We will continue to develop updates for kit14 models, but are unlikely to update R9 models at this time.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants