Skip to content

RagTag overly sensitive when identifying minimap2 SAM format flag #207

@ASLeonard

Description

@ASLeonard

Adding additional mm2 params of -xasm5 throws an error because RagTag sees a "-a" flag indicating SAM output from minimap. The actual argument is equivalent in minimap2's eyes as -x asm5 which is valid, because there is no "-" and "a" in the same string. Pretty annoying to catch all these cases, but raising in case anyone else hits this.

all_flags = "".join([i for i in self.params_string.split(" ") if i.startswith("-")])
if "a" in all_flags:
raise ValueError("Alignments must not be in SAM format (-a).")

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions