Skip to content

hsmetrics fails on non-human genomes #146

@percyfal

Description

@percyfal

We currently use agilents sureselect per default for calculating hs_metrics. hsmetrics calculations will fail for non-human genomes since the sequence dictionaries differ (i.e. chromosome lists):

Caused by: net.sf.samtools.util.SequenceUtil$SequenceListsDifferException: Sequence dictionaries are not the same size (1, 93)

Suggested solutions:

  1. run hs_metrics manually for relevant samples for now

  2. put target files in directory structure in which genome can be deduced, and if genome_build differs, don't run hs_metrics

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions