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

A. baumannii Pasteur scheme not in BLAST database #20

Closed
VGalata opened this issue Jun 6, 2016 · 1 comment
Closed

A. baumannii Pasteur scheme not in BLAST database #20

VGalata opened this issue Jun 6, 2016 · 1 comment
Assignees
Labels

Comments

@VGalata
Copy link

VGalata commented Jun 6, 2016

Dear Torsten,

It seems that the Pasteur A. baumannii scheme is not included in the BLAST database. Using the scheme abaumannii_2 returns no results:

bin/mlst --scheme abaumannii_2 <input FASTA> --debug

Retuned (after skipping matches for the Oxford A. baumannii scheme):

$VAR1 = [
          'abaumannii_2',
          '-',
          '-/-/-/-/-/-/-'
        ];

Also the inspection of the BLAST database returned no matches for the genes from the Pasteur scheme.

blastdbcmd -db db/blast/mlst.fa -entry all | grep 'Pas_.*'

Used version: mlst 2.5

@tseemann tseemann self-assigned this Jun 7, 2016
@tseemann tseemann added the bug label Jun 7, 2016
@tseemann
Copy link
Owner

tseemann commented Jun 7, 2016

@VGalata you are right - for some reason my script https://github.com/tseemann/mlst/blob/master/scripts/mlst-make_blast_db
does not include schemes with _2 in them.

I suspect I had a good reason at the time but I am not sure why now!

I have fixed the problem. Thank you for the report. Please use the latest github version (git clone).

@tseemann tseemann closed this as completed Jun 7, 2016
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants