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

Endless running Tombo resuiggle, reached to 100% in resuiggle but doesn't go to next step, generate output files #437

Open
kwonej0617 opened this issue Jul 13, 2023 · 1 comment

Comments

@kwonej0617
Copy link

kwonej0617 commented Jul 13, 2023

Hello! @marcus1487

I have run Tombo successfully with my dataset earlier. However, I used the same data and run the same pipeline again (I might add some options in resquiggle like corrected-group RawGenomeCorrected_000 --basecall-group Basecall_1D_000 --fit-global-scale --include-event-stdev --ignore-read-locks ), but it didn't work this time. In my error log file, it looks like the resquiggle process was successfully running and reached 100%, but it didn't go to the next step, which are tombs detet_modification step and text_output browser_files. No stats files, wig, bedgraph files generated. Could you please help me with how to solve this problem? Thank you so much!

#preprocess
#singularity exec /share/pkg/containers/tombo/tombo_1.5.1.sif tombo preprocess annotate_raw_with_fastqs --fast5-basedir ${wt_fast5} --fastq-filenames ${wt_fastq} --overwrite --processes 40

#resquiggling
singularity exec /share/pkg/containers/tombo/tombo_1.5.1.sif tombo resquiggle ${wt_fast5} ${ref} --rna --corrected-group RawGenomeCorrected_000 --basecall-group Basecall_1D_000 --overwrite --processes 40 --fit-global-scale --include-event-stdev --ignore-read-locks

#detect modifications
singularity exec /share/pkg/containers/tombo/tombo_1.5.1.sif tombo detect_modifications de_novo --fast5-basedirs ${wt_fast5} --statistics-file-basename wt --corrected-group RawGenomeCorrected_000 --processes 40

#output statistical results
singularity exec /share/pkg/containers/tombo/tombo_1.5.1.sif tombo text_output browser_files --fast5-basedirs ${wt_fast5} --statistics-filename wt.tombo.stats --browser-file-basename wt_rrach --genome-fasta ${ref} --motif-descriptions RRACH:3:m6A --file-types coverage dampened_fraction fraction --corrected-group RawGenomeCorrected_00

[12:18:53] Re-squiggling reads (raw signal to genomic sequence alignment).
^M 0%| | 0/1294636 [00:00<?, ?it/s]^M 0%| | 4/1294636 [00:01<90:04:40, 3.99it/s]^M 0%| | 16/1294636 [00:01<33:47:03, 10.64it/s]^M 0%| | 32/1294636 [00:02<22:31:29, 15.97it/s]^M 0%| | 59/1294636 [00:03<18:19:25, 19.62it/s]^M 0%| | 89/1294636 [00:04<16:11:43, 22.20it/s]^M 0%| | 131/1294636 [00:05<15:07:39, 23.77it/s]^M 0%| | 176/1294636 [00:07<14:19:49, 25.09it/s]^M 0%| | 232/1294636 [00:09<13:58:35, 25.73it/s]^M 0%| | 290/1294636 [00:11<13:39:54, 26.31it/s]^M 0%| | 358/1294636 [00:13<13:35:
...
[11:20:16<04:30, 31.51it/s]^M 99%|█████████▉| 1286300/1294636 [11:20:25<04:24, 31.51it/s]^M 99%|█████████▉| 1286400/1294636 [11:20:40<04:21,
31.50it/s]^M 99%|█████████▉| 1286500/1294636 [11:20:42<04:18, 31.50it/s]^M 99%|█████████▉| 1286600/1294636 [11:20:48<04:15, 31.50it/s]^M 99%|█████████▉| 1286700/1294636 [
11:20:49<04:11, 31.50it/s]^M 99%|█████████▉| 1286800/1294636 [11:20:56<04:08, 31.50it/s]^M 99%|█████████▉| 1286900/1294636 [11:20:57<04:05, 31.50it/s]^M 99%|█████████▉| 12
87100/1294636 [11:20:58<03:59, 31.50it/s]^M 99%|█████████▉| 1287300/1294636 [11:21:01<03:52, 31.50it/s]^M 99%|█████████▉| 1287500/1294636 [11:21:10<03:46, 31.50it/s]^M 99%
|█████████▉| 1287800/1294636 [11:21:12<03:36, 31.51it/s]^M 99%|█████████▉| 1288000/1294636 [11:21:30<03:30, 31.50it/s]^M 99%|█████████▉| 1288100/1294636 [11:21:31<03:27, 3
1.50it/s]^M100%|█████████▉| 1288200/1294636 [11:21:32<03:24, 31.50it/s]^M100%|█████████▉| 1288400/1294636 [11:21:36<03:17, 31.50it/s]^M100%|█████████▉| 1288600/1294636 [11
:21:40<03:11, 31.51it/s]^M100%|█████████▉| 1288800/1294636 [11:21:46<03:05, 31.51it/s]^M100%|█████████▉| 1289000/1294636 [11:21:48<02:58, 31.51it/s]^M100%|█████████▉| 1289
200/1294636 [11:21:55<02:52, 31.51it/s]^M100%|█████████▉| 1289400/1294636 [11:22:03<02:46, 31.51it/s]^M100%|█████████▉| 1289600/1294636 [11:22:08<02:39, 31.51it/s]^M100%|█
████████▉| 1289800/1294636 [11:22:14<02:33, 31.51it/s]^M100%|█████████▉| 1290000/1294636 [11:22:25<02:27, 31.51it/s]^M100%|█████████▉| 1290100/1294636 [11:22:36<02:24, 31.
50it/s]^M100%|█████████▉| 1290200/1294636 [11:22:39<02:20, 31.50it/s]^M100%|█████████▉| 1290400/1294636 [11:22:46<02:14, 31.50it/s]^M100%|█████████▉| 1290600/1294636 [11:2
2:47<02:08, 31.50it/s]^M100%|█████████▉| 1290800/1294636 [11:22:50<02:01, 31.51it/s]^M100%|█████████▉| 1291000/1294636 [11:22:56<01:55, 31.51it/s]^M100%|█████████▉| 129120
0/1294636 [11:23:04<01:49, 31.50it/s]^M100%|█████████▉| 1291400/1294636 [11:23:10<01:42, 31.50it/s]^M100%|█████████▉| 1291600/1294636 [11:23:13<01:36, 31.51it/s]^M100%|███
██████▉| 1291800/1294636 [11:23:24<01:30, 31.50it/s]^M100%|█████████▉| 1291900/1294636 [11:23:28<01:26, 31.50it/s]^M100%|█████████▉| 1292000/1294636 [11:23:30<01:23, 31.50
it/s]^M100%|█████████▉| 1292100/1294636 [11:23:39<01:20, 31.50it/s]^M100%|█████████▉| 1292200/1294636 [11:23:41<01:17, 31.50it/s]^M100%|█████████▉| 1292300/1294636 [11:23:
42<01:14, 31.50it/s]^M100%|█████████▉| 1292400/1294636 [11:23:46<01:10, 31.50it/s]^M100%|█████████▉| 1292500/1294636 [11:23:47<01:07, 31.50it/s]^M100%|█████████▉| 1292623/
1294636 [11:23:47<01:03, 31.51it/s]^M100%|█████████▉| 1293854/1294636 [11:23:48<00:24, 31.54it/s]^M100%|█████████▉| 1294325/1294636 [11:24:00<00:09, 31.54it/s]

@kwonej0617 kwonej0617 changed the title Endless running Tombo resuiggle, reached to 100% in resuiggle but doesn't go to next step, tombo_modifications Endless running Tombo resuiggle, reached to 100% in resuiggle but doesn't go to next step, generate output files Jul 13, 2023
@haowBio
Copy link

haowBio commented Sep 18, 2023

I encountered the same problem when running with docker. After installing it with conda, it ran normally.

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