This is a test version of Biostars. For the public version, visit https://www.biostars.org.
codeml analysis for multiple alignments

Hello Everyone,

I'd want to do a codeml analysis on many OG files created by orthofinder (including tree files). . I used "muscle" to construct the protein alignments, then "pal2nal" to convert all of the transcript into codon alignments. I was wondering if I should write several Codeml control files for each OG (over 1400) or if Codeml could be executed with a single ".ctl" file.!?

Suggestions appreciated. Thanks

paml codeml orthofinder
• 2,004 views • link
written by sbissi102 • sbissi102

2 answers

When I had such a task, I wrote a script that automatically made a lot of .ctl files, one file for each orthogroup, and then run Codeml with them.

• 0 views • link
written by shelkmike shelkmike

Thank you, In fact, I can't seem to locate the appropriate script.

• 0 views • link
written by sbissi102 • sbissi102

They explain something of the sort in this tutorial https://bioinformaticsworkbook.org/dataAnalysis/ComparativeGenomics/Finding_Positive_Selection_With_Codeml.html#gsc.tab=0

However, there's been an explosion of tools to automate paml and similar eg. https://pubmed.ncbi.nlm.nih.gov/31888452/

• 0 views • link
written by Roberto • Roberto

Thank you, Roberto. I attempted to follow the first instruction, but I am unable to generate the.clt files for the various orthogroups! In truth, I am unable to write the script.

• 0 views • link
written by sbissi102 • sbissi102

Log in to answer this question.

More posts like this