mahulchak / quickmerge

A simple and fast metassembler and assembly gap filler designed for long molecule based assemblies.
GNU General Public License v3.0
200 stars 31 forks source link

Output problems #66

Open Mondlii opened 3 years ago

Mondlii commented 3 years ago

Hi, I was using quickmerge on assemblies from Flye and Shasta. The run finished without any problems and I have an issue with the output files I received. The run gave me two fasta files, "merged_out.fasta" and "merged_prefix.fasta." This was the same for the anchor summary text files as well. I ran quast and busco on both the fasta outputs and it turns out the "merged_prefix.fasta" has a busco completion of 0% and an N50 of 180kb (which is lower than the one I had before merging) and on the other hand, the "merged_out.fasta" output has a completion score of 93% and an N50 of 500kb. The issue with the anchor files is that the one titled anchor_summary_prefix.txt is empty and only has headers, where as the "anchor_summary_out.txt" has a lot of information. Please advise on why this is and which files do you think I should keep moving forward or do I need a rerun? Thank you.

mahulchak commented 3 years ago

Could you please share your command?

On Wed, Jun 23, 2021, 04:47 Mondlii @.***> wrote:

Hi, I was using quickmerge on assemblies from Flye and Shasta. The run finished without any problems and I have an issue with the output files I received. The run gave me two fasta files, "merged_out.fasta" and "merged_prefix.fasta." This was the same for the anchor summary text files as well. I ran quast and busco on both the fasta outputs and it turns out the "merged_prefix.fasta" has a busco completion of 0% and an N50 of 180kb (which is lower than the one I had before merging) and on the other hand, the "merged_out.fasta" output has a completion score of 93% and an N50 of 500kb. The issue with the anchor files is that the one titled anchor_summary_prefix.txt is empty and only has headers, where as the "anchor_summary_out.txt" has a lot of information. Please advise on why this is and which files do you think I should keep moving forward or do I need a rerun? Thank you.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/mahulchak/quickmerge/issues/66, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABZQH2FTWWYXRKSQKYP73SDTUHCR7ANCNFSM47FT5RNQ .

Mondlii commented 3 years ago

merge_wrapper.py -l 1000 /home/Assess_Drafts/assembly.fasta /home/Assess_Drafts/Assembly.fasta

nucmer --mum -p out /home/Assess_Drafts/Assembly.fasta /home/Assess_Drafts/assembly.fasta
delta-filter -r -q -i 100 -l 1000 out.delta > out.rq.delta
quickmerge -d out.rq.delta -r Assess_Drafts/assembly.fasta -q Assess_Drafts/Assembly.fasta -hco 5 -c 2 -l 260000 -ml 5000 -p quick_merge