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

SHARED private genomes are not recognized by BLAST sometimes #2422

Open
Sveta-user opened this issue Sep 1, 2021 · 3 comments
Open

SHARED private genomes are not recognized by BLAST sometimes #2422

Sveta-user opened this issue Sep 1, 2021 · 3 comments
Assignees

Comments

@Sveta-user
Copy link

Greetings,
Sharing private genomes with other user(s) is a very useful feature in PATRIC. And it works smoothly by the most part, except with the Blast tool. Blast breaks (sometimes) if at least one genome in a group is a "shared private". Strangely this problem comes and goes: on some days shared genomes are recognized by BLAST, but other times they break Blast (like today)
The regular private genomes work reliably with blast all the time

In the example below I can view shared genome 1597.515, but not blast against it (not today at least ;))

Error creating BLAST FNA database for 1597.515 (owned by [email protected]): No contig sequences found in genome

image

image

@olsonanl
Copy link

olsonanl commented Sep 1, 2021

This is a different issue than what we were seeing before; in this genome the contig sequences are tagged with sequence types of plasmid and chromosome, where the blast service is looking for sequence types of contig for building the database.

Have contig blasts been working for you with other genomes? I'll need to look into updating the code to properly pull the contigs.

@olsonanl olsonanl self-assigned this Sep 1, 2021
@Sveta-user
Copy link
Author

Interesting... Didn't realize that genome ownership or sharing was not the only issue with blast. Now that you pointed it out, I looked at many genome that work fine with BLAST - and nearly all of them have "contig" as sequence type. But there is one exception: 391904.5 has "chromosome" , but works with blast nonetheless -
image

While 1597.515 that does NOT work with blast has a "plasmid" and 2 "chromosome" labeled contigs. Can it be that only the presence of a "plasmid" contig breaks blast..?

image

Thank you for looking into this issue, Bob!

@Sveta-user
Copy link
Author

image

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

2 participants