You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
When I run abeona v0.40.0 on my data, then I get the following error message:
INFO:cortexpy:Log level is INFO
INFO:cortexpy.traverse:Loading graph: g29549.traverse.ctx
INFO:cortexpy.traverse:Loaded 3346 kmers
INFO:cortexpy.traverse:Making graph consistent
INFO:cortexpy.traverse:Loading links file g29549.ctp.gz
INFO:cortexpy.traverse:Exiting after element 100000
INFO:cortexpy.traverse:Writing seq records to -
INFO:cortexpy.graph.serializer.unitig:Collapsing graph with 3346 kmers
INFO:cortexpy.graph.serializer.unitig:Collapsing complete. Collapsed graph contains 87 unitigs
INFO:cortexpy.graph.interactor:Found 3 incoming tip nodes
INFO:cortexpy.graph.interactor:Found 5 outgoing tip nodes
INFO:cortexpy.graph.interactor:Incoming node 0; 5 outgoing nodes; Path number 0
Traceback (most recent call last):
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/bin/cortexpy", line 11, in <m
sys.exit(main())
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
return getattr(module, method_string)(args.args)
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
for record in seq_record_generator:
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
for idx, val in enumerate(iterator):
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
for rec in seq_record_generator:
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
for pidx, path in enumerate(paths):
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
stack.append(iter(G[child]))
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
successors = list(parent_walker.successors())
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
j_unitigs = list(self.link_successors())
File "/data/winni/proj/olof/abeona_paper_analysis/real_data/.snakemake/conda/e505603e/lib/python3.6/site-packages/c
junctions: {self.link_walker.junctions}"""
ValueError: Links do not appear to match unitigs. Have these links been constructed on a different cortex graph?
current unitig: {'repr': 'GGAGG', 'unitig': 'AGGGAGGAGGAGGGAGAAGGAGGGAGGGAAGAGGAGGGAGAAGGAGGGAGG', 'c
successors: [{'repr': 'G', 'unitig': 'GGAGGAGGGAGAAGGAGGGAGGGAAGAGGAGGGAGAAGGAGGGAGGG', 'coverage': (
available bases: {'G', 'A'}
junctions: defaultdict(<class 'list'>, {'A': ['AGGAG', 'AGGAG'], 'G': ['G', 'GGAGGGGGGGAGGG', 'GGAG',
The text was updated successfully, but these errors were encountered:
I am still digging into this issue, but it appears to be a problem in how cortexpy traverses links in this graph. The subgraph is low complexity and contains many cycles. For now, I suggest using the --no-links argument to turn off link traversal.
When I run abeona v0.40.0 on my data, then I get the following error message:
The text was updated successfully, but these errors were encountered: