Hello everyone, usually, when I have an alignment of nucleotide, in order to get a global dNdS value for the whole alignment, I use sqinr and the function kaks which calculates a pairwise dNdS between each sequences, then I just calculate the mean of all the dNdS values and here is my global dNdS ratio.

But here it is just a punctual value.

And I wondered if there was techniques (program) that allows to calculate a dNdS and that gives a pvalue associated in order ton conclude that the dNdS in significantly different from 0 for exemple.

Thank your for you help.