miniasm: asm.c:267: ma_ug_seq: Assertion `sub[id].e - sub[id].s <= ks->seq.l' failed. about miniasm HOT 5 CLOSED

You may have reads with same names or same part splited by space. Check it like this:
zcat input.fa.gz | perl -ne 'print "$1\n" if />(\S+)/' | sort | uniq -d

from miniasm.

Thank you, I used rename.sh in=input.fasta out=output.fasta prefix=test and it fixed the names problem

from miniasm.

For others who might encounter this issue too, I thought I'd just add here that rename.sh is a script in the BBtools package: https://jgi.doe.gov/data-and-tools/bbtools/bb-tools-user-guide/

from miniasm.

So what is the likely issue when one encounters this error but you don't have duplicated IDs?

@RyanBio @mictadlo @roblanf

[M::main] ===> Step 5: generating unitigs <===
miniasm: asm.c:267: ma_ug_seq: Assertion `sub[id].e - sub[id].s <= ks->seq.l' failed.
/bin/bash: line 1: 13225 Aborted                 (core dumped) miniasm -f filtered/DC_concat.fq.gz overlaps/DC_concat.paf.gz > miniasm/DC_concat.gfa

from miniasm.

@mw55309 I don't know. I have only seen this assertion failure when there are duplicated read names.

from miniasm.