Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Duplicate fasta headers cause assertion error #1

Open
lmfaber opened this issue Nov 4, 2019 · 0 comments
Open

Duplicate fasta headers cause assertion error #1

lmfaber opened this issue Nov 4, 2019 · 0 comments
Labels
bug Something isn't working

Comments

@lmfaber
Copy link
Owner

lmfaber commented Nov 4, 2019

Causes this error when creating the read graph from the equivalence class:

2019-11-04 17:14:49 [INFO]: Perform clustering with HDBSCAN. (min_cluster_size: 2)
2019-11-04 17:14:51 [INFO]: Salmon mapping...
2019-11-04 17:15:11 [INFO]: Build full graph...
Traceback (most recent call last):
  File "/mnt/prostlocal/lasse/miniconda3/envs/test_karma/bin/karma.py", line 208, in <module>
    full_graph = ReadGraph.from_equivalence_classes(eqc_file, sequences)
  File "/data/prostlocal/lasse/miniconda3/envs/test_karma/bin/read_graph.py", line 78, in from_equivalence_classes
    assert no_of_contigs == len(max_number_of_mapped_reads.keys())
AssertionError

It would be nice to implement a check for this when reading the fasta file.

@lmfaber lmfaber added the bug Something isn't working label Nov 4, 2019
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant