Description
Describe the bug
NextDenovo stuck in the raw_align phase (both db_stat and db_split have been completed) and the nextdenovo.sh.e log shows that Minimap2 continues to output mapped 1 sequences, but there is no error message.Now that you can eliminate the problem with system resources (enough to use) and no error messages are reported, the nextdenovo process just disappeared on its own, I wonder why
Error message
here is my output information:
(base) [handanlu@cn3 01.nextdenovo]$ cat nextdenovo_nohup.out
nohup: ignoring input
[61142 INFO] 2025-03-17 21:34:56 NextDenovo start...
[61142 INFO] 2025-03-17 21:34:56 version:2.5.2 logfile:pid61142.log.info
[61142 INFO] 2025-03-17 21:34:56 skip mkdir: /data/nextdenovo_out
[61142 WARNING] 2025-03-17 21:34:56 backup /data/nextdenovo_out/01.raw_align to /data/nextdenovo_out/01.raw_align.backup.v0
[61142 INFO] 2025-03-17 21:34:56 mkdir: /data/nextdenovo_out/01.raw_align
[61142 WARNING] 2025-03-17 21:34:56 backup /data/nextdenovo_out/02.cns_align to /data/nextdenovo_out/02.cns_align.backup.v0
[61142 INFO] 2025-03-17 21:34:56 mkdir: /data/nextdenovo_out/02.cns_align
[61142 WARNING] 2025-03-17 21:34:56 backup /data/nextdenovo_out/03.ctg_graph to /data/nextdenovo_out/03.ctg_graph.backup.v0
[61142 INFO] 2025-03-17 21:34:56 mkdir: /data/nextdenovo_out/03.ctg_graph
[61142 INFO] 2025-03-17 21:35:01 Total jobs: 1
[61142 INFO] 2025-03-17 21:35:01 Submitted jobID:[61172] jobCmd:[/data/nextdenovo_out/01.raw_align/01.db_stat.sh.work/db_stat1/nextDenovo.sh] in the local_cycle.
[61142 INFO] 2025-03-17 21:42:09 db_stat done
[61142 INFO] 2025-03-17 21:42:09 updated options:
rerun: 3
task: all
deltmp: 1
read_type: ont
job_type: local
input_type: raw
parallel_jobs: 1
read_cutoff: 20k
seed_depth: 45.0
pa_correction: 1
seed_cutfiles: 3
genome_size: 470m
seed_cutoff: 88771
job_prefix: nextDenovo
blocksize: 41247962893
ctg_cns_options: -p 10
nextgraph_options: -a 1
input_fofn: /data/ont.fofn
workdir: /data/nextdenovo_out
minimap2_options_map: -x map-ont
sort_options: -m 100g -t 20 -k 40 -k 40
raw_aligndir: /data/nextdenovo_out/01.raw_align
cns_aligndir: /data/nextdenovo_out/02.cns_align
ctg_graphdir: /data/nextdenovo_out/03.ctg_graph
minimap2_options_raw: -t 20 -x ava-ont --min-occ-floor=500 -K 1000
correction_options: -p 10 -max_lq_length 20000 -r ont -min_len_seed 100000
minimap2_options_cns: -t 20 -x ava-ont -k 21 -w 21 --minlen 20000 --maxhan1 50000
[61142 INFO] 2025-03-17 21:42:09 summary of input data:
file: /data/nextdenovo_out/01.raw_align/input.reads.stat
[Read length stat]
Types Count (#) Length (bp)
N10 14094 198757
N20 34153 155761
N30 58765 130588
N40 87541 113180
N50 120401 100000
N60 157360 89357
N70 198567 80384
N80 244375 72223
N90 295458 64608
Types Count (#) Bases (bp) Depth (X)
Raw 352830 34898987631 74.25
Filtered 0 0 0.00
Clean 352830 34898987631 74.25
*Suggested seed_cutoff (genome size: 470.00Mb, expected seed depth: 45, real seed depth: 45.00): 88771 bp
[61142 INFO] 2025-03-17 21:42:14 Total jobs: 1
[61142 INFO] 2025-03-17 21:42:14 Submitted jobID:[62141] jobCmd:[/data/nextdenovo_out/01.raw_align/02.db_split.sh.work/db_split1/nextDenovo.sh] in the local_cycle.
[61142 INFO] 2025-03-17 21:50:57 db_split done
[61142 INFO] 2025-03-17 21:50:57 Total jobs: 9
[61142 INFO] 2025-03-17 21:50:57 Submitted jobID:[67488] jobCmd:[/data/nextdenovo_out/01.raw_align/03.raw_align.sh.work/raw_align1/nextDenovo.sh] in the local_cycle.
This is the contents of my cfg file:
[General]
job_type = local # local, slurm, sge, pbs, lsf
job_prefix = nextDenovo
task = all # all, correct, assemble
rewrite = no # yes/no
deltmp = yes
parallel_jobs = 1 # 4个并行任务(共使用约80核)
input_type = raw # raw, corrected
read_type = ont # clr, ont, hifi
input_fofn = ont.fofn
workdir = nextdenovo_out
[correct_option]
read_cutoff = 20k
genome_size = 470m # estimated genome size
pa_correction = 1
sort_options = -m 100g -t 20 -k 40
minimap2_options_raw = -t 20 -x ava-ont --min-occ-floor=500 -K 1000 # 超多线程+内存优化 # P/parallel_jobs
correction_options = -p 10 -max_lq_length 20000 -r ont -min_len_seed 100000 # P/pa_correction
[assemble_option]
minimap2_options_cns = -t 20 -x ava-ont -k 21 -w 21 --minlen 20000 --maxhan1 50000 # -t P/parallel_jobs
nextgraph_options = -a 1
Genome characteristics
My genome size (470M) and ONT data characteristics (70x coverage,N50=98k,Q=19)