This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bray-Curtis dissimilarity for comparing genomes, how to normalize?

Hi everyone,

I am comparing whole genomes of a closely related species using Nucmer. I want to calculate the %of similarity based on the bases aligned from each genome using Bray-Curtis dissimilarity (I had also tried Sorensen coefficient). The problem is that some genomes have very different genome sizes (let's say ranging from 40 Mb to 170 Mb). In those cases, I get similarity values above 1 (or 100%) which is impossible. I have tried some normalizations as those recommended by Somerfield 2008, and Yoshioka (2008) but nothing worked.

Some suggestions? alternatives?

genome mummer bray-curtis normalization sorensen
• 1,631 views • link
written by Buffo 246

I suggest you try FastANI. It compares sequences without actual alignment by calculating k-mer similarity, which in most cases is related to sequence identity. The upside for your purposes is that sequences of different lengths can be used.

• 0 views • link
written by Mensur Dlakic 3.0K

FastANI is faster than others programs such as Nucmer but it has a few disadvantages such as:

.- Is not recommended for fragmented or short genomes (N50 be ≥10 Kbp).
.- No ANI output is reported for a genome pair if ANI value is much below 80%.

But in my case would be a good alternative. Thanks.

• 0 views • link
written by Buffo 246

0 answers

No answers yet.

Log in to answer this question.

More posts like this