ValueError: End location (8780) must be greater than or equal to start location (9428)
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
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- 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