steineggerlab / Metabuli

Metabuli: specific and sensitive metagenomic classification via joint analysis of DNA and amino acid.
GNU General Public License v3.0
118 stars 10 forks source link

Database creation #73

Open VinaykiyaVridhi opened 3 months ago

VinaykiyaVridhi commented 3 months ago

I am trying to create a custom database with nt fastas, but getting a new error each time. Is there a way I could use my Kraken db to run metabuli?

Additionally, I tried running metabuli with one of your databases, in my output file how do I get both the nt column but also the aa column.

Appreciate any help!

Regards!

jaebeom-kim commented 3 months ago

Thank you for reaching out! Let me fix the problem of nt DB creation. Kraken DB is not compatible with Metabuli.

Related to the output question, nt and aa are not printed currently. I think I can amend Metabuli to print nt, but print aa might be more complex.

Thanks!

VinaykiyaVridhi commented 3 months ago

Thank you for your response! In your paper you have mentioned about the specificity in protein level. How can I use metabuli to do a protein level analysis?

jaebeom-kim commented 3 months ago

Metabuli performs protein level search and DNA level search at the same time. So you're already conducting the protein level analysis!

VinaykiyaVridhi commented 2 months ago

Thank you for responding. Where exactly can I find the protein level output, because when I ran the tool, I get only one output file which shows nt related results.

jaebeom-kim commented 2 months ago

Sorry, I misunderstood what you said. Metabuli does not report which protein/gene the query sequence is derived from. It only reports taxon. I am working on Metabuli2 to report proteins/genes. Please look forward to the next version.

VinaykiyaVridhi commented 2 months ago

Thank you for your prompt response.

Unfortunately, I am still unable to view the taxon at the protein level. Could you please advise if there is a specific parameter that needs to be added to include this information in the output?