issues about finalizing with run-asm-pipeline-post-review.sh

[zjuzexian]

Hi all,
I meet a problem when finalizing the output after re-order the chromosome by juicer.
And every time when after filling gaps. (finalize-output.sh)
The program died with the temp.fasta collaped with error like this:


And the .FINAL.fasta is much smaller(900Mb) than the input .final.fasta(1.5G).

And I checked the .final.cprops, it may died after the last scaf in temp.fa(scaf_30_RagTag:::fragment_5:::debris) and then reach these scaffolds:
scaf_30_RagTag:::fragment_6:::overhang_1 730 923994
scaf_30_RagTag:::fragment_6:::overhang_2:::gap 731 6
scaf_30_RagTag:::fragment_7:::debris:::overhang_1:::gap 732 94
scaf_30_RagTag:::fragment_7:::debris:::overhang_2 733 10906
And I checked in the .finale cprops there is no other scaffolds with "overhang", only fragment or debris.

I have checked that for .final.fasta, these sequneces(scaf_30_RagTag:::fragment_6:::overhang_1 toscaf_30_RagTag:::fragment_7:::debris:::overhang_2 ) are all exist.

I wonder why it happen and how can I fix this?

Thanks!

[zjuzexian]

It seems that the FINAL.assembly file before my adjust habe already had these fragment overhang, which doesn't exist in the previous .final.assembly file

[dudcha]

My guess is that you are reviewing the wrong file in Juicebox. JB should load rawchrom.hic and rawchrom.assembly (final.hic and final.assembly in previous versions). These files should not contain gapncontigs and overhangs, that’s one step before.

…

On Apr 8, 2022, at 8:12 AM, zjuzexian ***@***.***> wrote:  It seems that after I changed the chromosome border with juicebox, the .review.assembly file changed some fragment into overhang, which doesn't exist in the previous .assembly file — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you are subscribed to this thread.

[cui-yutong]

I met the same problem. Did you get the solution?
Thanks,

[zjuzexian]

Yes, just as dudcha said, you should check if you are using the right file from juicebox.