Skip to content

Add fallback REF entries for loci missing these in VCF #96

@dnil

Description

@dnil

Is your feature request related to a problem? Please describe.
Generally it is a Good Idea that the STR caller identifies the number of reference copies detected given a reference genome, locus catalog and caller version. But we could certainly have a fallback when these are missing, e.g. TRGT
PacificBiosciences/trgt#71.

Describe the solution you'd like
Parse a default REF entry from locus catalogs and use it for output iff there is no INFO.REF field set for the variant.

Describe alternatives you've considered
Stranger could in principle calculate this from the reference region and allele structure definition given, but for more complex cases would need a reference genome sequence, and this baseline guesstimate could deviate from the implementation in the caller.

Metadata

Metadata

Assignees

No one assigned

    Labels

    enhancementNew feature or request

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions