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

Error!Deepmod2 analysis of human methylation #29

Open
zzangcy opened this issue Sep 2, 2024 · 2 comments
Open

Error!Deepmod2 analysis of human methylation #29

zzangcy opened this issue Sep 2, 2024 · 2 comments

Comments

@zzangcy
Copy link

zzangcy commented Sep 2, 2024

Hi,
I encountered an error when I used the test BAM for analysis. It can't be resolved. I'm worried. Please help me.
And our server does not have a GPU,The python version is v3.12,the code and error message are in the file.
Thanks a lot,
Freya

dm2-error.log

@umahsn
Copy link
Collaborator

umahsn commented Sep 2, 2024

Hi,

How did you basecall your FAST5 files? Can you share the exact command? DeepMod2 needs move table data either in BAM file or in the FAST5 file, depending upon how you basecall. It seems like your BAM file does not have move tables. If the move tables are in FAST5 file (only if you basecalled with Guppy v <6.4), use --fast5_move parameter.

If neither FAST5, nor the BAM file has move tables, you would need to rebasecall your reads using --moves_out parameter in Guppy, as shown here: https://github.com/WGLab/DeepMod2/blob/main/docs/Example.md#2-methylation-calling-from-fast5-files-with-guppy-basecalling. If you basecall with Guppy v>6.4, then the move table will be in BAM file and DeepMod2 will detect it and you wouldn't need to use --fast5_move parameter.

Umair

@zzangcy
Copy link
Author

zzangcy commented Sep 3, 2024

Hi,

How did you basecall your FAST5 files? Can you share the exact command? DeepMod2 needs move table data either in BAM file or in the FAST5 file, depending upon how you basecall. It seems like your BAM file does not have move tables. If the move tables are in FAST5 file (only if you basecalled with Guppy v <6.4), use --fast5_move parameter.

If neither FAST5, nor the BAM file has move tables, you would need to rebasecall your reads using --moves_out parameter in Guppy, as shown here: https://github.com/WGLab/DeepMod2/blob/main/docs/Example.md#2-methylation-calling-from-fast5-files-with-guppy-basecalling. If you basecall with Guppy v>6.4, then the move table will be in BAM file and DeepMod2 will detect it and you wouldn't need to use --fast5_move parameter.

Umair

Hi,
Because the FAST5 files comes from a doctor who knows nothing about sequencing, I can't provide specific information about basecall, that's ok,I fully understand your explanation now👍, which is very important to me.Thanks so mush!🌼
Freya

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