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
- I am using Python2.7 and often use Biopython
- My users are likely to run the programs from Windows machines