malonge / RagTag

Tools for fast and flexible genome assembly scaffolding and improvement
MIT License
470 stars 47 forks source link

Scaffolding error: 'RuntimeError: There are no useful alignments. Check output alignment files.` #145

Closed bahiyahazli closed 1 year ago

bahiyahazli commented 1 year ago

Hello.

I encountered some problems when trying to scaffold. This is my script:

ragtag.py scaffold -w -u --debug -t 32 ~/wgs-bahiyah/pre/contig-bwa/ref/KT862811.fasta ~/wgs-bahiyah/pre/scaffold-ragtag/ragtag.correct.fasta -o ~/wgs-bahiyah/pre/scaffold-ragtag/

Below are my python and related dependencies info: Python3 Unimap, minimap2 and nucmer function well mummer==3.23 numpy==1.24.0 intervaltree==3.10.0 networkx==2.8.8 pysam==0.20.0

And here is my output:

Sat Dec 24 12:50:03 2022 --- VERSION: RagTag v2.1.0 Sat Dec 24 12:50:03 2022 --- CMD: ragtag.py scaffold -w -u --debug -t 32 /home/bahiyah/wgs-bahiyah/pre/contig-bwa/ref/KT862811.fasta /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ragtag.correct.fasta -o /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ Sat Dec 24 12:50:03 2022 --- INFO: Mapping the query genome to the reference genome Sat Dec 24 12:50:03 2022 --- INFO: Overwriting pre-existing file: /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ragtag.scaffold.asm.paf Sat Dec 24 12:50:03 2022 --- INFO: Running: minimap2 -x asm5 -t 32 /home/bahiyah/wgs-bahiyah/pre/contig-bwa/ref/KT862811.fasta /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ragtag.correct.fasta > /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ragtag.scaffold.asm.paf 2> /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ragtag.scaffold.asm.paf.log Sat Dec 24 12:50:03 2022 --- INFO: Finished running : minimap2 -x asm5 -t 32 /home/bahiyah/wgs-bahiyah/pre/contig-bwa/ref/KT862811.fasta /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ragtag.correct.fasta > /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ragtag.scaffold.asm.paf 2> /home/bahiyah/wgs-bahiyah/pre/scaffold-ragtag/ragtag.scaffold.asm.paf.log Sat Dec 24 12:50:03 2022 --- INFO: Reading whole genome alignments Sat Dec 24 12:50:03 2022 --- INFO: Filtering and merging alignments Traceback (most recent call last): File "/home/bahiyah/miniconda3/envs/bioinfo/bin/ragtag_scaffold.py", line 4, in import('pkg_resources').run_script('RagTag==2.1.0', 'ragtag_scaffold.py') File "/home/bahiyah/miniconda3/envs/bioinfo/lib/python3.8/site-packages/pkg_resources/init.py", line 672, in run_script self.require(requires)[0].run_script(script_name, ns) File "/home/bahiyah/miniconda3/envs/bioinfo/lib/python3.8/site-packages/pkg_resources/init.py", line 1479, in run_script exec(script_code, namespace, namespace) File "/home/bahiyah/miniconda3/envs/bioinfo/lib/python3.8/site-packages/RagTag-2.1.0-py3.8.egg/EGG-INFO/scripts/ragtag_scaffold.py", line 577, in File "/home/bahiyah/miniconda3/envs/bioinfo/lib/python3.8/site-packages/RagTag-2.1.0-py3.8.egg/EGG-INFO/scripts/ragtag_scaffold.py", line 488, in main RuntimeError: There are no useful alignments. Check output alignment files.

Do you mind checking where the hiccup is? I've re-git clone, re-pip install and even re-python3 setup.py install, few times already. Additionally, I've tried to ./run_run_tests.sh in the /tests/integration_tests/ folder and the test gives out error as well. Probably because of this same issue.

Thank you so much. Will dig more onto what's actually happening.

EmilieSmeets22 commented 1 year ago

Hey @bahiyahazli Did you figure out the issue? I just got the exact same error.

Thanks already

Sameer1028 commented 1 year ago

I am facing the same problem. Can someone please provide their insights?

malonge commented 1 year ago

This is not an error with RagTag. It's just saying that there were no alignments that it could use to inform scaffolding. To debug, please analyze the alignment files themselves. And, as you have done, I do recommend using --debug so that you can check out the intermediate alignment files as well.