Skip to content

protein evidence (blast or diamond) - minos_metrics_blastp_tophit fails  #56

@dadrasarmin

Description

@dadrasarmin

Hi!

I tried to feed a bunch of protein fasta files as evidence to the minos and I noticed a part of my pipeline fails always. After investigation, I noticed that minos_metrics_blastp_combine finished successfully but minos_metrics_blastp_tophit was resubmitted 3 times due to Input files updated by another job which was not true. I checked the time stamp and this job was submitted too early (almost at the same time as minos_metrics_blastp_combine). I am using minos on a local machine with a huge amount of resources. Maybe you cannot reproduce the issue with an HPC. After trying many things, this seems to solve the problem for me. I wanted to report it for future investigation by your team.

Changing this

run:
from minos.scripts.get_blast_tophit import get_blast_tophit
get_blast_tophit(input[0], output[0], params.pident_threshold, params.qcov_threshold)

To this:

    shell:
        """
        python -c 'import os; from minos.scripts.get_blast_tophit import get_blast_tophit; get_blast_tophit("{input}", "{output}", {params.pident_threshold}, {params.qcov_threshold})'
        """

Best regards,
Armin

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions