8000 R2DT fails with complex header lines · Issue #64 · r2dt-bio/R2DT · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

R2DT fails with complex header lines #64

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

Closed
blakesweeney opened this issue Jan 12, 2022 · 0 comments · Fixed by #157
Closed

R2DT fails with complex header lines #64

blakesweeney opened this issue Jan 12, 2022 · 0 comments · Fixed by #157

Comments

@blakesweeney
Copy link
Contributor

It can't handle some complex header lines like:

>DB|TEXT|MORE Some nice name

It seems that having several | characters causes esl-sfetch to fail with:

Creating SSI index for /rna/r2dt/temp/sequence.fasta...    done.
Indexed 1 sequences (1 names).
SSI index written to file /rna/r2dt/temp/sequence.fasta.ssi
sh: 1: TEXT: not found
sh: 1: MORE: not found
seq DB not found in SSI index for file /rna/r2dt/temp/sequence.fasta

# R2DT :: visualise RNA secondary structure using templates
# Version 1.2 (August 10, 2021)
# https://github.com/RNAcentral/R2DT
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# R2DT :: visualise RNA secondary structure using templates
# Version 1.2 (August 10, 2021)
# https://github.com/RNAcentral/R2DT
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Visualising sequence DB|TEXT|MORE using the HS_LSU_3D model from ribovision_lsu
Traceback (most recent call last):
  File "/rna/r2dt/r2dt.py", line 583, in <module>
    cli()
  File "/usr/local/lib/python3.4/dist-packages/click/core.py", line 764, in __call__
    return self.main(*args, **kwargs)
  File "/usr/local/lib/python3.4/dist-packages/click/core.py", line 717, in main
    rv = self.invoke(ctx)
  File "/usr/local/lib/python3.4/dist-packages/click/core.py", line 1137, in invoke
    return _process_result(sub_ctx.command.invoke(sub_ctx))
  File "/usr/local/lib/python3.4/dist-packages/click/core.py", line 956, in invoke
    return ctx.invoke(self.callback, **ctx.params)
  File "/usr/local/lib/python3.4/dist-packages/click/core.py", line 555, in invoke
    return callback(*args, **kwargs)
  File "/usr/local/lib/python3.4/dist-packages/click/decorators.py", line 17, in new_func
    return f(get_current_context(), *args, **kwargs)
  File "/rna/r2dt/r2dt.py", line 163, in draw
    force_draw(force_template, fasta_input, output_folder, seq_id)
  File "/rna/r2dt/r2dt.py", line 517, in force_draw
    ribovision.visualise('lsu', fasta_input, output, seq_id, model_id)
  File "/rna/r2dt/utils/ribovision.py", line 51, in visualise
    raise ValueError("Failed esl-sfetch for: %s" % rnacentral_id)
ValueError: Failed esl-sfetch for: DB|TEXT|MORE

A complete fix of the issue is to transform all fasta headers to something that will always work (1, 2, 3) and then run the steps on this transformed file. Before creating the final results for the user all ids should be replaced with the original ones.

@github-project-automation github-project-automation bot moved this to Backlog in R2DT Issues Oct 4, 2024
@AntonPetrov AntonPetrov moved this from Backlog to In progress in R2DT Issues Oct 13, 2024
@github-project-automation github-project-automation bot moved this from In progress to Done in R2DT Issues Dec 22, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
Status: Done
Development

Successfully merging a pull request may close this issue.

1 participant
0