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

minia generates non-deterministic assemblies? #25

Closed
morispi opened this issue Jun 21, 2021 · 2 comments
Closed

minia generates non-deterministic assemblies? #25

morispi opened this issue Jun 21, 2021 · 2 comments

Comments

@morispi
Copy link

morispi commented Jun 21, 2021

Hello,

I'm attempting to use minia to perform local assembly of Linked-Reads data, and realign the obtained contigs to determine breakpoints as a part of a SV calling tool.

However, when running multiple instances of minia, with the same parameters, and on the same input file, it seems to output (slightly) different assemblies.

While all assemblies are very similar, some of them sometimes end up messing up the breakpoints computation (although I can definitely take some of the blame on that problem).

I've tried playing around with the parameters a bit, but couldn't come up with something that would yield deterministic assemblies. So I'm wondering, is that an expected behaviour, or maybe am I overlooking something here?

Thanks.

Best,
Pierre

@rchikhi
Copy link
Member

rchikhi commented Jun 25, 2021

hi Pierre, I recall some simplifications operations are done in parallel, the best way to get deterministic results is therefore to use a single thread, hopefully your instance is small enough.. In fact, see section "Reproducibility" in manual.tex.

@morispi
Copy link
Author

morispi commented Jun 30, 2021

Hello Rayan,

Indeed, it seems like I missed the manual.tex information. :) Works well when using -nb-cores 1 -nb-glue-partitions 200 as recommended. My instances are small enough, so this seems perfect.

Thanks a lot!
Pierre

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