I am looking to systematically scan genomes for balancing selection. However, I am going to be doing this analysis within a species, with 66 replicates. In doing so, I have been recommended to stay away from programs that typically are used for analyses across species.
Does anyone have a suggestion for dN/dS software to use for between population analysis?
Thanks for your time.
1 answer
Although PAML is probably the gold standard for this kind of analysis, it is kind of awkward to script around it and its output is considered painful to parse (I never personally did so). The ETE Toolkit is really nice and flexible for programmatically performing such tasks, with the added benefit of offering good documentation and pretty graphical output.
Log in to answer this question.
-
How codeml detects positive selection
written by Amay •
-
CODEML: Ideal setting for generating dn/ds values for a gene family within a species.
written by venura
-
Identify fixed differences between population in WGS data (VCF Format)
written by goatsrunfaster •
-
Problem with CODEML parameters
written by aeonswhisper •
-
Calculate dN/dS ratio for protein coding genes using PAML
written by b.shrestha •
-
dN/dS within and between species
written by krc3004 •
-
Multiple Sequence Alignment for very large data sets
written by dthorbur
-
Is there any tool to calculate dn/ds for orthologues?
written by ThulasiS
-
Gene-specific selection with counting method when no synonymous subsitutions
written by ksw •
-
is Hapmap or 1000 genome VCF data is from diseased or healthy one
written by Being Bioinformatician
I'm no expert but I'm wondering how should the computation be different for inter and intra-species analyses and why a typical tool like PAML wouldn't be suitable. You can check this discussion to see if it helps: Best Practices/Softwares To Calculate Ka/Ks Ratio
I think this issue has been resolved. It appears to have been a misunderstanding between a bioinformatician and myself. I believe they were concerned about the ML tree, rather than actually calculating the ratio. Since I have such a big data set, we are going to use RAxML or ExaML to generate the ML tree and either HyPhy or potentially the ETE Toolkit to calculate the dN/dS ratio. Though, a lot of background reading is necessary before choosing.
thanks for following up