Skip to content

Include instrument name in output metadata TSV #80

@etal

Description

@etal

PacBio instruments include their name in the header of the unaligned BAM they produce. For each read group (there may be multiple), denoted by the @RG tag , the PL key has the sequencing platform (e.g. PACBIO) and the PM key has the platform model (e.g. REVIO or SEQUELII).

Without requiring additional user input, these values could be extracted from the input BAM and emitted as an additional column in the output metadata TSV.

Fragile example code:

rg = pysam.AlignmentFile(path, "rb").header["RG"][0]
instrument_name = f"{rg['PL']} {rg['PM']}"

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions