Skip to content

Graph distance logic for skipping repetitive anchors in Giraffe doesn't account for the anchor length itself #4701

@adamnovak

Description

@adamnovak

Here we're supposed to be looking at the length of the space between seeds in the graph if we skip a seed.

if (next_skippable->is_skippable() && next_skippable_it+1 != chain.end() &&
total_graph_distance+graph_distance < this->max_skipped_bases) {

But we're just summing up the graph distance from the previous anchor to this one, and from this one to the next anchor. We aren't accounting for the length of the anchor itself either for the comparison here, or for the update to total_graph_distance later. So we're systematically underestimating the net distance.

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