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

[BUG] IF crashes on a multiple contigs file #103

Closed
1 of 3 tasks
eloilit opened this issue Sep 29, 2022 · 1 comment
Closed
1 of 3 tasks

[BUG] IF crashes on a multiple contigs file #103

eloilit opened this issue Sep 29, 2022 · 1 comment
Assignees
Labels

Comments

@eloilit
Copy link
Contributor

eloilit commented Sep 29, 2022

Describe the bug
IF crashes on a multiple contigs file, using the gembase format, in the middle of the contigs being processed.

To Reproduce
Steps to reproduce the behavior:

  1. integron_finder --local-max --mute --outdir ./IFResults --gembase --gembase-path ./VibrionaceaDBMarch2022_v2 --gbk ./VibrionaceaDBMarch2022_v2/Replicons/VIBR.0322.11443.fna
  2. See error

I located the files to reproduce the error in /pasteur/appa/scratch/elittner/for_issue_IF_crashes
Expected behavior

Normal completion

Screenshots

Traceback (most recent call last):
File "/Users/eloi/Documents/Softwares/anaconda3/envs/IFv2.0.1/bin/integron_finder", line 10, in
sys.exit(main())
File "/Users/eloi/Documents/Softwares/anaconda3/envs/IFv2.0.1/lib/python3.9/site-packages/integron_finder/scripts/finder.py", line 648, in main
integron_res, summary = find_integron_in_one_replicon(replicon, config)
File "/Users/eloi/Documents/Softwares/anaconda3/envs/IFv2.0.1/lib/python3.9/site-packages/integron_finder/scripts/finder.py", line 385, in find_integron_in_one_replicon
integrons = find_integron(replicon, protein_db, integron_max, intI_file, phageI_file, config)
File "/Users/eloi/Documents/Softwares/anaconda3/envs/IFv2.0.1/lib/python3.9/site-packages/integron_finder/integron.py", line 150, in find_integron
attc_left = np.array([i_attc.pos_beg.values[0] for i_attc in attc_ac])
File "/Users/eloi/Documents/Softwares/anaconda3/envs/IFv2.0.1/lib/python3.9/site-packages/integron_finder/integron.py", line 150, in
attc_left = np.array([i_attc.pos_beg.values[0] for i_attc in attc_ac])
IndexError: index 0 is out of bounds for axis 0 with size 0

Please complete the following information):

OS:

  • Linux
  • Windows
  • 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
Only one example failing over a big dataset

@eloilit eloilit added the bug label Sep 29, 2022
bneron added a commit that referenced this issue Oct 13, 2022
when attc site overlap intI and there is only one attc

by the way rename some variables to improve find_integron readibility
and add functional test to cover thie case
@eloilit
Copy link
Contributor Author

eloilit commented Oct 14, 2022

VibrionaceaDBMarch2022_v2.zip

command used: integron_finder --local-max --mute --outdir ./IFResults --gembase --gembase-path ./VibrionaceaDBMarch2022_v2 --gbk ./VibrionaceaDBMarch2022_v2/Replicons/VIBR.0322.11443.0157.fna

@jeanrjc Here are the files to reproduce the error

@bneron bneron closed this as completed Oct 13, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants