Sliding Window Ka/Ks (Dn/Ds) Analysis Via Python?

0

Entering edit mode

10.3 years ago

a1ultima ▴ 840

I am working in Python and performing Ka/Ks (dN/dS) ratio analysis.

I would like to ask for any suggestions for how to go about performing "within gene" or "sliding window" Ka/Ks analysis, such that a Ka/Ks ratio can be given for positions along the gene sequence.

This MATLAB tutorial really sums up exactly what I need, but I am obliged by my boss to use Python.

Anticipated output example

enter image description here

Extra info

I am using Python2.7 and often use Biopython
My users are likely to run the programs from Windows machines

evolution protein • 7.4k views

ADD COMMENT • link 10.3 years ago by a1ultima ▴ 840

1

Entering edit mode

Doing sliding window dn/ds analyses is not a good idea: http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0003746

ADD REPLY • link 10.3 years ago by karl.schmid ▴ 20

0

Entering edit mode

thank you very much for this gem @karl.schmid, but indeed we plan to use sliding windows more for visual understanding as the paper suggests

ADD REPLY • link 10.3 years ago by a1ultima ▴ 840

0

Entering edit mode

Hi, 5 years later I have the exact same problem...did you figure it out? Thanks in advance

ADD REPLY • link 5.1 years ago by mariamcurbelo2061 • 0

Login before adding your answer.