pangenome / odgi

Optimized Dynamic Genome/Graph Implementation: understanding pangenome graphs
https://doi.org/10.1093/bioinformatics/btac308
MIT License
193 stars 39 forks source link

atomicbitvector assertion error in untangle #489

Open esrice opened 1 year ago

esrice commented 1 year ago

When running this command:

odgi untangle -i k_locus.og -r 'bGalGal1b#0#chrZ:11159192-11400491' -m 256 -P > test.bed

I receive the following error:

[odgi::untangle] warning: no step index specified. Building one with a sample rate of 8. This may take additional time. A step index can be provided via -a, --step-index. A step index can be built using odgi st
epindex.
[odgi::algorithms::stepindex] Collecting Steps Progress: 18.37% @ 1.84e+03/s elapsed: 00:00:00:00 remain: 00:[odgi::algorithms::stepindex] Collecting Steps Progress: 100.00% @ 1.06e+03/s elapsed: 00:00:00:00 re
main: 00:00:00:00
[odgi::algorithms::stepindex] Building Progress: 100.00% @ 1.05e+03/s elapsed: 00:00:00:00 remain: 00:00:00:00
[odgi::algorithms::untangle] untangling 294 queries with 1 targets
[odgi::algorithms::untangle] establishing initial cuts for 294 paths
odgi: /opt/conda/conda-bld/odgi_1671316371355/work/deps/atomicbitvector/include/atomic_bitvector.hpp:159: bool atomicbitvector::atomic_bv_t::set(size_t, std::memory_order): Assertion `idx < _size' failed.
Aborted

Running v0.8.2-0-g8715c55. Let me know if there's any other information that would be helpful. Thanks!

ekg commented 1 year ago

I would guess that there is a mismatch between yout BED file and the graph paths. Have you confirmed they match? You may need to apply odgi procbed if the graph is a subset of a larger one but the BED is not. Just a thought.

To be sure we can understand the best thing is to share the test case. Is this possible?

On Wed, Mar 29, 2023, 23:03 Edward S. Rice @.***> wrote:

When running this command:

odgi untangle -i k_locus.og -r 'bGalGal1b#0#chrZ:11159192-11400491' -m 256 -P > test.bed

I receive the following error:

[odgi::untangle] warning: no step index specified. Building one with a sample rate of 8. This may take additional time. A step index can be provided via -a, --step-index. A step index can be built using odgi st epindex. [odgi::algorithms::stepindex] Collecting Steps Progress: 18.37% @ 1.84e+03/s elapsed: 00:00:00:00 remain: 00:[odgi::algorithms::stepindex] Collecting Steps Progress: 100.00% @ 1.06e+03/s elapsed: 00:00:00:00 re main: 00:00:00:00 [odgi::algorithms::stepindex] Building Progress: 100.00% @ 1.05e+03/s elapsed: 00:00:00:00 remain: 00:00:00:00 [odgi::algorithms::untangle] untangling 294 queries with 1 targets [odgi::algorithms::untangle] establishing initial cuts for 294 paths odgi: /opt/conda/conda-bld/odgi_1671316371355/work/deps/atomicbitvector/include/atomic_bitvector.hpp:159: bool atomicbitvector::atomic_bv_t::set(size_t, std::memory_order): Assertion `idx < _size' failed. Aborted

Running v0.8.2-0-g8715c55. Let me know if there's any other information that would be helpful. Thanks!

— Reply to this email directly, view it on GitHub https://github.com/pangenome/odgi/issues/489, or unsubscribe https://github.com/notifications/unsubscribe-auth/AABDQEJDPSLJ7WA724N2QEDW6SPR5ANCNFSM6AAAAAAWMNFDHE . You are receiving this because you are subscribed to this thread.Message ID: @.***>

esrice commented 1 year ago

Thanks for the fast reply. I wasn't using a range from a bed file, but rather the argument to -r is the name of the path from "odgi paths" as recommended by the tutorial "Untangling the pangenome" under the "Linearize the C4 region" heading. The og is indeed a subgraph of a larger graph, but I'm trying to untangle using the whole path of bGalGal1b through the subgraph, as the tutorial does for the hg38 path. How might I go about doing this? Thanks!

On Thu, Mar 30, 2023, 7:36 AM Erik Garrison @.***> wrote:

I would guess that there is a mismatch between yout BED file and the graph paths. Have you confirmed they match? You may need to apply odgi procbed if the graph is a subset of a larger one but the BED is not. Just a thought.

To be sure we can understand the best thing is to share the test case. Is this possible?

On Wed, Mar 29, 2023, 23:03 Edward S. Rice @.***> wrote:

When running this command:

odgi untangle -i k_locus.og -r 'bGalGal1b#0#chrZ:11159192-11400491' -m 256 -P > test.bed

I receive the following error:

[odgi::untangle] warning: no step index specified. Building one with a sample rate of 8. This may take additional time. A step index can be provided via -a, --step-index. A step index can be built using odgi st epindex. [odgi::algorithms::stepindex] Collecting Steps Progress: 18.37% @ 1.84e+03/s elapsed: 00:00:00:00 remain: 00:[odgi::algorithms::stepindex] Collecting Steps Progress: 100.00% @ 1.06e+03/s elapsed: 00:00:00:00 re main: 00:00:00:00 [odgi::algorithms::stepindex] Building Progress: 100.00% @ 1.05e+03/s elapsed: 00:00:00:00 remain: 00:00:00:00 [odgi::algorithms::untangle] untangling 294 queries with 1 targets [odgi::algorithms::untangle] establishing initial cuts for 294 paths odgi: /opt/conda/conda-bld/odgi_1671316371355/work/deps/atomicbitvector/include/atomic_bitvector.hpp:159: bool atomicbitvector::atomic_bv_t::set(size_t, std::memory_order): Assertion `idx < _size' failed. Aborted

Running v0.8.2-0-g8715c55. Let me know if there's any other information that would be helpful. Thanks!

— Reply to this email directly, view it on GitHub https://github.com/pangenome/odgi/issues/489, or unsubscribe < https://github.com/notifications/unsubscribe-auth/AABDQEJDPSLJ7WA724N2QEDW6SPR5ANCNFSM6AAAAAAWMNFDHE

. You are receiving this because you are subscribed to this thread.Message ID: @.***>

— Reply to this email directly, view it on GitHub https://github.com/pangenome/odgi/issues/489#issuecomment-1489726831, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABLL5JJD4BPYJOS25QHYY2DW6ULUNANCNFSM6AAAAAAWMNFDHE . You are receiving this because you authored the thread.Message ID: @.***>