[BUG] IF crashes with --gbk option in a specific caase
Describe the bug When using the --gbk option, IF crashes if a feature (e.g. cassette, attC site) spans over both borders of a replicon.
To Reproduce Steps to reproduce the behavior: With the attached files
- integron_finder --local-max --mute --outdir ./IFResults --gembase --gembase-path ./Gembase --gbk ./Gembase/Replicons/KLPN001.0321.00476.P003.fst
- See error
Expected behavior
Finish without error and produce output files
Screenshots

Please complete the following information):
OS:
- [ ] Linux
- [ ] Windows
- [X] Mac
Integron_Finder Version:
integron_finder version 2.0.2 Using:
-
Python 3.9.13 | packaged by conda-forge | (main, May 27 2022, 17:01:00) [Clang 13.0.1 ]
-
numpy 1.19.4
-
pandas 1.1.5
-
matplolib 3.3.3
-
biopython 1.78
-
Prodigal V2.6.3: February, 2016
-
INFERNAL 1.1.4 (Dec 2020)
-
HMMER 3.3.2 (Nov 2020); http://hmmer.org/
Additional context Not a big issue, as it works properly if we remove the --gbk option
a feature (e.g. cassette, attC site) spans over both borders of a replicon
You mean, that the feature exists in the circular replicon and is cut in 2 in the linear sequence ?
Yes ! Here's an extraction of the .integrons file in my specific case. You can see that the last protein detected begins at position 190392 and ends at position 299
