marbl / merqury

k-mer based assembly evaluation
Other
272 stars 19 forks source link

I don't understand how to get the right K mer size #120

Closed Axze-rgb closed 3 months ago

Axze-rgb commented 3 months ago

Hello, I know it seems obvious and easy, but if I use your example with the human genome:

best_k.sh 6.2GB 0.001 
genome: 6.2GB
tolerable collision rate: 0.001
6.2983

It's not 21

What am I missing?

Thanks a lot

arangrhie commented 3 months ago

My best_k.sh isn't smart to parse out the . and GB. Please use 6200000000 for 6.2 GB.

Axze-rgb commented 3 months ago

Ahhh, thank you :)