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

Using GTDB-TK for a databas of partial bacterial genome #426

Closed
mkoohim opened this issue Aug 21, 2022 · 2 comments
Closed

Using GTDB-TK for a databas of partial bacterial genome #426

mkoohim opened this issue Aug 21, 2022 · 2 comments
Labels
stale Awaiting submitter feedback.

Comments

@mkoohim
Copy link

mkoohim commented Aug 21, 2022

Hello,

I have a database biosynthetic gene cluster (BGCs) from different bacteria. BGCs are clusters of adjacent genes which is part of whole genome of the bacteria. I used GTDB-TK to make a phylogeny tree for these BGCs but unfortunately as its not whole genome I couldn't get any results and in identify step there was not any marker to match with these BGCs.
I am wondering if there is any way to use GTDB-TK on these BGCs database?
Ive attached here one of BGC file:

seq1.zip

@aaronmussig
Copy link
Member

Hi, unfortunately the genome will need hits to marker genes since GTDB-Tk needs hits to either the 53, or 120 single-copy marker genes (for archaea and bacteria respectively).

You may be able to use the ani_rep command of GTDB-Tk to try and find the closest representative genome in ANI space, but note you may need a less-stringent alignment fraction value (this can be set with --min_af)

@aaronmussig aaronmussig added the stale Awaiting submitter feedback. label Sep 30, 2022
@pchaumeil
Copy link
Collaborator

Ticket closed due to inactivity.

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

No branches or pull requests

3 participants