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

Filter Bad MLST Profiles from PubMLST #328

Open
Nilad opened this issue Aug 18, 2022 · 0 comments
Open

Filter Bad MLST Profiles from PubMLST #328

Nilad opened this issue Aug 18, 2022 · 0 comments

Comments

@Nilad
Copy link

Nilad commented Aug 18, 2022

Hi,

Recently, i notice that few Profile like ST-2724 in "Acinetobacter baumannii#1" schema or ST-2609 in "Haemophilus influenzae" schema contains "N" letter index instead of allele index.

So i got this error when i want to download theses schemas...

ariba pubmlstget "Haemophilus influenzae" mlst_hinfluenza_test
WARNING: spades not found in path. Looked for spades.py
Traceback (most recent call last):
  File "/usr/local/bin/ariba", line 312, in <module>
    args.func(args)
  File "/usr/local/lib/python3.8/dist-packages/ariba/tasks/pubmlstget.py", line 11, in run
    preparer.run()
  File "/usr/local/lib/python3.8/dist-packages/ariba/pubmlst_ref_preparer.py", line 81, in run
    self.profile = mlst_profile.MlstProfile(profile_file, duplicate_warnings=True)
  File "/usr/local/lib/python3.8/dist-packages/ariba/mlst_profile.py", line 15, in __init__
    self._load_input_file()
  File "/usr/local/lib/python3.8/dist-packages/ariba/mlst_profile.py", line 29, in _load_input_file
    type_tuple = tuple(int(row[x]) for x in self.genes_list)
  File "/usr/local/lib/python3.8/dist-packages/ariba/mlst_profile.py", line 29, in <genexpr>
    type_tuple = tuple(int(row[x]) for x in self.genes_list)
ValueError: invalid literal for int() with base 10: 'N'

Maybe, it can possible to check if a no numeric value is present and remove the corresponding ST profile...

Thanks in advance.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant