Question: Ka/Ks (Or Dn/Ds) Calculation From Sequenced Transcriptome Contigs
gravatar for tom.newport
7.8 years ago by
tom.newport10 wrote:

I have around 25,000 contigs each for two species. I'm trying to calculate Ka and Ks (or Dn and Ds) for potential orthologs between species, and potential paralogs within species. I have a list of contigs which I want to calculate Ka and Ks for, which looks like:

Species1_Contig12, Species1_Contig98,0

(The number signifies the comparison being made - 0 = species 1 to species 1, 1=species 1 to species 2, 2 = species 2 to species 2)

I also have a FASTA file with the consensus sequences for each contig. Names are the same in the FASTA file as in the comparison file, and those sequences have been annotated against other sequences by BLASTx, which suggests which reading frame I should be looking at.

I would like to calculate Ka and Ks between the pairs in the file, but I have no idea how to do it - PAML has been suggested, but the data isn't in the right format.

Can anyone offer me any pointers on how to go about this?

paml • 4.2k views
ADD COMMENTlink modified 13 months ago by Adrian Pelin2.4k • written 7.8 years ago by tom.newport10
gravatar for Raygozak
7.8 years ago by
State College, PA, Penn State
Raygozak1.3k wrote:

You can look into the Los Alamos tool ( for conservation analysis, and for the alignment, this might help you: . You want to look into the codon suite to see if it helps you.

ADD COMMENTlink modified 7.8 years ago by Manu Prestat4.0k • written 7.8 years ago by Raygozak1.3k
gravatar for Adrian Pelin
13 months ago by
Adrian Pelin2.4k
Adrian Pelin2.4k wrote:

If you can generate a translational alignment for each of your suspected orthologous group, you can use codeml

ADD COMMENTlink written 13 months ago by Adrian Pelin2.4k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2056 users visited in the last hour