Merck / Merck/deepbgc

ValueError: End location (8780) must be greater than or equal to start location (9428)

Open
#98 0 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Dominant language
Jupyter Notebook
Stars
161
Forks
29
PR merge metrics
No merged PRs in 30d

Description

Hi,

I'm running version 0.1.31 and I'm getting the ValueError

part of error message is below:

Processing record #221974: k141_5772
Updating record molecule type to DNA
Finding genes in record: k141_5772
Reusing already existing HMMER hmmscan result: 111X30.10/tmp/k141_5772.pfam.domtbl.txt
End location (8780) must be greater than or equal to start location (9428)
Traceback (most recent call last):
File "/lustre/home/2101111992/miniconda3/envs/deepbgc/lib/python3.7/site-packages/deepbgc/main.py", line 113, in main
run(argv)
File "/lustre/home/2101111992/miniconda3/envs/deepbgc/lib/python3.7/site-packages/deepbgc/main.py", line 102, in run
args.func.run(**args_dict)
File "/lustre/home/2101111992/miniconda3/envs/deepbgc/lib/python3.7/site-packages/deepbgc/command/pipeline.py", line 177, in run
step.run(record)
File "/lustre/home/2101111992/miniconda3/envs/deepbgc/lib/python3.7/site-packages/deepbgc/pipeline/annotator.py", line 35, in run
pfam_annotator.annotate()
File "/lustre/home/2101111992/miniconda3/envs/deepbgc/lib/python3.7/site-packages/deepbgc/pipeline/pfam.py", line 122, in annotate
location = self._get_pfam_loc(best_hsp.query_start, best_hsp.query_end, protein)
File "/lustre/home/2101111992/miniconda3/envs/deepbgc/lib/python3.7/site-packages/deepbgc/pipeline/pfam.py", line 55, in _get_pfam_loc
return FeatureLocation(start, end, strand=feature.strand)
File "/lustre/home/2101111992/miniconda3/envs/deepbgc/lib/python3.7/site-packages/Bio/SeqFeature.py", line 800, in init
f"End location ({self.end}) must be greater than "
ValueError: End location (8780) must be greater than or equal to start location (9428)

DeepBGC failed with ValueError: End location (8780) must be greater than or equal to start location (9428)
k141_5772.pfam.domtbl.txt
k141_5772.pfam.proteins.txt
k141_5772.prodigal.nucl.txt
k141_5772.prodigal.proteins.txt
Thank you for your help!

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start in deepbgc/pipeline/pfam.py, especially _get_pfam_loc, using the traceback and the attached HMMER and Prodigal output files as inputs. Reproduce the failure with the k141_5772 data and trace how the reported query coordinates reach FeatureLocation. Done means the failing input is handled without the ValueError and the pipeline completes annotation for this case.

Written by the indexing model from the issue text.

Assessment

Tech stack
python
Domain
bioinformatics
Issue type
Bug
Difficulty
3/5
Estimated time
1-2 days
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
35/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.