Skip to content

Assertion `newEdgeNodes.count(std::make_pair(reverse(pos), reverse(before))) == 1 || newEdgeNodes.count(std::make_pair(before, pos)) == 1' failed. #63

@olekto

Description

@olekto

Hi,
I'm running with MBG 1.0.10, but with verkko 1.0 beta 2 since the 1.0 release haven't been updated in Conda yet.

removed 1 tips
try resolve k=1601, replaced 2 nodes with 5 nodes, unitigified 8 nodes to 4 nodes
try resolve k=1602, replaced 6 nodes with 17 nodes, unitigified 26 nodes to 9 nodes
try resolve k=1603, replaced 7 nodes with 23 nodes, unitigified 36 nodes to 16 nodes
try resolve k=1604, replaced 10 nodes with 26 nodes, unitigified 32 nodes to 16 nodes
try resolve k=1605, replaced 2 nodes with 6 nodes, unitigified 8 nodes to 4 nodes
try resolve k=1606, replaced 5 nodes with 12 nodes, unitigified 14 nodes to 7 nodes
try resolve k=1607, replaced 3 nodes with 9 nodes, unitigified 12 nodes to 6 nodes
try resolve k=1608, replaced 3 nodes with 10 nodes, unitigified 15 nodes to 5 nodes
try resolve k=1609, replaced 3 nodes with 10 nodes, unitigified 16 nodes to 6 nodes
try resolve k=1610, replaced 5 nodes with 15 nodes, unitigified 24 nodes to 10 nodes
try resolve k=1611, replaced 3 nodes with 7 nodes, unitigified 10 nodes to 5 nodes
try resolve k=1612, replaced 7 nodes with 17 nodes, unitigified 24 nodes to 12 nodes
try resolve k=1613MBG: src/UnitigResolver.cpp:1816: ResolutionResult resolve(ResolvableUnitigGraph&, const HashList&, size_t, std::vector<PathGroup>&, const phmap::flat_hash_set<long unsigned int>&, size_t, bool): Assertion `newEdgeNodes.count(std::make_pair(reverse(pos), reverse(before))) == 1 || newEdgeNodes.count(std::make_pair(before, pos)) == 1' failed.
./buildGraph.sh: line 39: 26304 Aborted                 (core dumped) /cluster/projects/nn8013k/programs/miniconda3/envs/verkko/bin/MBG $iopt -t 4 -k 1001 -r 15000 -R 4000 -w 100 --kmer-abundance 1 --unitig-abundance 2 --error-masking=collapse-msat --output-sequence-paths ../1-buildGraph/paths.gaf --out ../1-buildGraph/hifi-resolved.gfa

I got the same issue with 1.0.9, but wanted to test with 1.0.10 before reporting.

This is for an Atlantic salmon. Trying with only HiFi reads so far. I could share data if you want (55 GB uncompressed corrected data).

Thank you.

Ole

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