How To Compute Dn/Ds Ratio And Get Descriptive Statistics On Thousands Of Sequences ?
2
3
Entering edit mode
12.5 years ago
Bioch'Ti ★ 1.1k

Hi Guys,

I'm looking for the best way to (1) compute dN/dS ratios on thousands couples of sequences and (2) get descriptive statistics on these thousands of sequences (length, GC3, etc...) at the same time?

What are the software available to perform that job?

Cheers

selection statistics next-gen snp • 9.0k views
ADD COMMENT
1
Entering edit mode

To clarify: http://en.wikipedia.org/wiki/Ka/Ks_ratio You are referring to the ratio of Non-synonymous vs. Synonymous substitution, correct?

ADD REPLY
6
Entering edit mode
12.5 years ago
hadasa ★ 1.0k

Have a look at HyPhy(Hypothesis testing using Phylogenies). http://www.datam0nk3y.org/hyphy/doku.php It has a specilized module for Positive and negative selection detection and which can be relevant to what you want to accomplish for a large dataset of sequences It might help you out.

ADD COMMENT
4
Entering edit mode
12.5 years ago
Dave Lunt ★ 2.0k

You could use BioPerl to do this. Calculating dN:dS can be done in a BioPerl script by running PAML. See this page http://www.bioperl.org/wiki/HOWTO:PAML and slide 76 of this presentation http://jason.open-bio.org/Bioperl_Tutorials/ISMB2007/. To count the GC there is a script here http://www.bioperl.org/wiki/Bioperl_scripts#SeqStats and to modify it by codon position should be quite possible.

ADD COMMENT
0
Entering edit mode

Or indeed, any of the Bio* projects: they all are well suited to building simple sequence analysis pipelines for many sequences.

ADD REPLY

Login before adding your answer.

Traffic: 1269 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6