Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
  • P PathoFact
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
    • Locked Files
  • Issues 29
    • Issues 29
    • List
    • Boards
    • Service Desk
    • Milestones
    • Requirements
  • Merge requests 0
    • Merge requests 0
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
    • Test Cases
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages and registries
    • Packages and registries
    • Package Registry
    • Container Registry
    • Infrastructure Registry
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Code review
    • Insights
    • Issue
    • Repository
  • Wiki
    • Wiki
  • External wiki
    • External wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • Laura Denies
  • PathoFact
  • Issues
  • #72
Closed
Open
Issue created May 19, 2021 by Michal Zeman@valzip

Test workflow and common run fails at RGI step at biopython

Hi! I was testing your software and I wasn't able to pass the test workflow (same problem is in normal workflow).

Data were single WGS genome in separate contigs.

Is it problem of my dependencies or some typo in the code?

Tested system:

Linux Mint 19.3 Tricia based on Ubuntu 18.04 bionic

conda 4.10.1

python 3.8.3.final.0

although when I have conda environment active

python version is 3.6.4

and biopython 1.78

Error in rule run_RGI:
    jobid: 46
    output: /home/kuskus/ANALYSIS/pathofact/PathoFact_results/PathoFact_intermediate/PathoFact_intermediate/AMR/RGI_results/Cl3778/group_1.RGI.txt
    log: /home/kuskus/ANALYSIS/pathofact/PathoFact_results/PathoFact_intermediate/logs/Cl3778/group_1.RGI.log (check log file(s) for error message)
    conda-env: /home/kuskus/Software/PathoFact/.snakemake/conda/efe64355
    shell:
        rgi main --input_sequence /home/kuskus/ANALYSIS/pathofact/PathoFact_results/PathoFact_intermediate/splitted/Cl3778/group_1.fasta --output_file /home/kuskus/ANALYSIS/pathofact/PathoFact_results/PathoFact_intermediate/PathoFact_intermediate/AMR/RGI_results/Cl3778/group_1.RGI --input_type protein --local --clean  &> /home/kuskus/ANALYSIS/pathofact/PathoFact_results/PathoFact_intermediate/logs/Cl3778/group_1.RGI.log
        (exited with non-zero exit code)

Log file contains several lines:

Traceback (most recent call last):
  File "/home/kuskus/Software/PathoFact/.snakemake/conda/efe64355/bin/rgi", line 2, in <module>
    from app.MainBase import MainBase
  File "/home/kuskus/Software/PathoFact/.snakemake/conda/efe64355/lib/python3.6/site-packages/app/MainBase.py", line 1, in <module>
    from app.settings import *
  File "/home/kuskus/Software/PathoFact/.snakemake/conda/efe64355/lib/python3.6/site-packages/app/settings.py", line 7, in <module>
    from Bio.Alphabet import generic_dna
  File "/home/kuskus/.local/lib/python3.6/site-packages/Bio/Alphabet/__init__.py", line 21, in <module>
    "Bio.Alphabet has been removed from Biopython. In many cases, the alphabet can simply be ignored and removed from scripts. In a few cases, you may need to specify the ``molecule_type`` as an annotation on a SeqRecord for your script to work correctly. Please see https://biopython.org/wiki/Alphabet for more information."
ImportError: Bio.Alphabet has been removed from Biopython. In many cases, the alphabet can simply be ignored and removed from scripts. In a few cases, you may need to specify the ``molecule_type`` as an annotation on a SeqRecord for your script to work correctly. Please see https://biopython.org/wiki/Alphabet for more information.
Assignee
Assign to
Time tracking