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

User-defined blast directory? #6

Open
LovellHAGSC opened this issue Dec 21, 2020 · 2 comments
Open

User-defined blast directory? #6

LovellHAGSC opened this issue Dec 21, 2020 · 2 comments

Comments

@LovellHAGSC
Copy link

Dear Dr. Derelle,
Thank you for your work with Broccoli.
Quick question, is there support for a user-supplied blast directory, like orthofinder -b? I'd like to port my orthology pipelines to Broccoli and compare with orthofinder, but to do this, I'd need to feed it the identical blast results.
Thanks,
John

@rderelle
Copy link
Owner

Dear John,

Short answer: I'm afraid such option isn't available.

Firstly, Broccoli is not compatible with Blast, which has a different cigar output format, and can only be ran with Diamond.
Secondly, Broccoli deletes all Diamond output files to save disk-space, and there is no implemented option to restart it with already available Diamond output.

If using Diamond is ok with your task, you would need to modify the script corresponding to step 2 in order to not perform a de-novo Diamond search and instead to use the files present in dir_step2.
I could help you with that.

best,
Romain

@LovellHAGSC
Copy link
Author

Thank you.
Yeah, diamond is fine. I'll shoot you an email and we can discuss so that this thread doesn't get crazy.
Cheers,
John

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