Skip to content

This issue was moved to a discussion.

You can continue the conversation there. Go to discussion →

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

“Disentangling failed: 'Multiple isolated embplant_pt components detected! Broken or contamination?” AND the assemble length range different from reference #308

Closed
GaoXY99 opened this issue Jan 17, 2024 · 2 comments

Comments

@GaoXY99
Copy link

GaoXY99 commented Jan 17, 2024

Hi, I have some problems during chloroplast genome assembling and need your help. I ran several times, upload the following two commands, log, and the results. Thank you very much for your attentiion.

COMMAND1:
get_organelle_from_reads.py -1 seq/PYA1/PYA1_clean_R1.fastq.gz -2 seq/PYA1/PYA1_clean_R2.fastq.gz -o temp/round50_k12_seedetc/PYA1_new -t 24 -F embplant_pt -R 50 -k 21,29,39,45,59,65,79,85,99,105,119,141
parameter1
parameter1_get_org.log.txt

COMMAND2:
get_organelle_from_reads.py -1 seq/PYA1/PYA1_clean_R1.fastq.gz -2 seq/PYA1/PYA1_clean_R2.fastq.gz -o temp/round50_k12_seedetc/PYA1_new -t 24 -F embplant_pt -R 50 -k 21,29,39,45,59,65,79,85,99,105,119,141 -w 85 --max-reads inf --reduce-reads-for-coverage inf
parameter2
parameter2_get_org.log.txt

  1. Disentangling fail reported in these two log file
    Disentangling failed: 'Multiple isolated embplant_pt components detected! Broken or contamination?'
    Disentangling failed: 'No new connections.'
    How to deal with it?

  2. assembled sequence length of quadripartite is obviously different from the reference
    eg: assembled SSC sequence length is 13kb, the SSC of reference is 18kb
    assembled IR (32 kb) is longer than the reference(27kb).

  3. It seemed that "-w 85 --max-reads inf --reduce-reads-for-coverage inf" improve the depth, however, it has not reduced the contig number.

@JianjunJin
Copy link
Collaborator

Thank you for your query. Based on your runs, it appears that the organelle sufficient graph has already been achieved, so adjusting parameters may not be necessary (refer to our FAQ for more details).

The primary challenge is the disentanglement of graphs with extremely high depth contigs. We have a new branch addressing this issue, although it's not yet officially released on conda. You can check it out here - the intermediate_graph branch. Since it's not available on bioconda, please install using Option 2 or 3. Don't hesitate to contact me for any assistance.

@smallfishcui
Copy link

I encountered the same error while using the assembly done by flye. I adjusted the parameter --min-depth 100 --max-depth 3000 and it finally completed the analyses. I changed the parameters for multiple times...it doesn't work if the range is smaller or larger than these.
best,
Cui

Repository owner locked and limited conversation to collaborators Apr 3, 2024
@JianjunJin JianjunJin converted this issue into discussion #324 Apr 3, 2024

This issue was moved to a discussion.

You can continue the conversation there. Go to discussion →

Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants