-
Notifications
You must be signed in to change notification settings - Fork 12
Description
Describe the bug
I used Juicebox to edit an assembly and generated the review.assembly file.
Using that file, I tried to generate a new scaffold (fasta, agp, bed) with the following :
python juicebox_assembly_converter.py -a <modified_assembly> -f <original_premodification_fasta>
Unfortunately juicebox_assembly_converter.py complains that the naming is wrong (see below).
Checking for breaks listed in .assembly and making them...
Traceback (most recent call last):
File "../../juicebox_scripts/juicebox_assembly_converter.py", line 875, in
simple_chr_names=simple_chr_names)
File "../../juicebox_scripts/juicebox_assembly_converter.py", line 91, in process
sequences = self._add_breaks(sequences, assembly_map)
File "../../juicebox_scripts/juicebox_assembly_converter.py", line 300, in _add_breaks
raise BadContigNameError("Unbroken contig {} failed to map!!".format(fragment_name))
main.BadContigNameError: Unbroken contig ctg00000001.1 failed to map!!
ctg00000001.1 doesn't exist in my assembly.
Indeed the review file seems to contain wrong contig names: "ctgXXXXXX" instead of the "h1tgXXXXXXX" of the original premodification FASTA file.
Is this a bug of Juicebox or of the converter script?
Thanks for your help.
Laurent Falquet
OS
MacOS Sequoia, Juicebox v.2.15