Question: Ka/Ks (Or Dn/Ds) Calculation From Sequenced Transcriptome Contigs
1
gravatar for tom.newport
6.8 years ago by
tom.newport10
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
Species2_Contig16_Species1_Contig24,1

(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 • 3.7k views
ADD COMMENTlink modified 6 weeks ago by Adrian Pelin2.3k • written 6.8 years ago by tom.newport10
1
gravatar for Raygozak
6.8 years ago by
Raygozak1.3k
State College, PA, Penn State
Raygozak1.3k wrote:

You can look into the Los Alamos tool (http://www.hiv.lanl.gov/content/sequence/SNAP/SNAP.html) for conservation analysis, and for the alignment, this might help you: http://www.expasy.org/genomics/sequence_alignment . You want to look into the codon suite to see if it helps you.

ADD COMMENTlink modified 6.8 years ago by Manu Prestat3.9k • written 6.8 years ago by Raygozak1.3k
0
gravatar for Adrian Pelin
6 weeks ago by
Adrian Pelin2.3k
Canada
Adrian Pelin2.3k wrote:

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

ADD COMMENTlink written 6 weeks ago by Adrian Pelin2.3k
Please log in to add an answer.

Help
Access

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