-
Notifications
You must be signed in to change notification settings - Fork 32
Closed
Milestone
Description
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
Labels
No labels