Skip to content
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

TSEBRA failed with keep_ab_initio.cfg (but worked with default.cfg) #28

Open
Huiting120 opened this issue Feb 7, 2023 · 1 comment
Open

Comments

@Huiting120
Copy link

Huiting120 commented Feb 7, 2023

Hi,
I was running TSEBRA with both the default.cfg and keep_ab_initio.cfg. The first run with the default.cfg was successful, however, the run the keep_ab_initio.cfg failed with the following error:

$ ../TSEBRA/bin/tsebra.py -g braker1_augustus.hints.gtf,braker2_augustus.hints.gtf -e braker1_hintsfile.gff,braker2_hintsfile.gff -c keep_ab_initio.cfg -o braker.combined.gtf
### READING GENE PREDICTION: [braker1_augustus.hints.gtf]
### READING GENE PREDICTION: [braker2_augustus.hints.gtf]
### READING EXTRINSIC EVIDENCE: [braker1_hintsfile.gff]
### READING EXTRINSIC EVIDENCE: [braker2_hintsfile.gff]
### BUILD OVERLAP GRAPH
### ADD FEATURES TO TRANSCRIPTS
### SELECT TRANSCRIPTS
Traceback (most recent call last):
  File "/home/tech/Documents/TSEBRA/keep_ab_initio/../TSEBRA/bin/tsebra.py", line 228, in <module>
    main()
  File "/home/tech/Documents/TSEBRA/keep_ab_initio/../TSEBRA/bin/tsebra.py", line 105, in main
    combined_prediction = graph.get_decided_graph()
  File "/home/tech/Documents/TSEBRA/TSEBRA/bin/overlap_graph.py", line 433, in get_decided_graph
    self.decide_graph()
  File "/home/tech/Documents/TSEBRA/TSEBRA/bin/overlap_graph.py", line 408, in decide_graph
    self.edges[key].node_to_remove = self.decide_edge(self.edges[key])
  File "/home/tech/Documents/TSEBRA/TSEBRA/bin/overlap_graph.py", line 346, in decide_edge
    if diff > self.para[f'e_{i+1}']:
KeyError: 'e_5'

The default.cfg has e_5 and e_6, but the keep_ab_initio.cfg does not. But I've ran other annotations without e_5 and e_6 pre-set in the configuration file before and this error did not occur. I wonder what caused the error this time. Any insights would be helpful.

Thanks in advance.

@Huiting120
Copy link
Author

I see the scripts and default.cfg was updated last week. Does the new script require e_5 and e_6? If so, can you provide an updated keep_ab_initio.cfg please?
Thanks

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant