How To Convert Protein (Plus Dna) Into Cdna Sequence?

Question

[Michael Kuhn]

I have a eukaryotic protein sequence and the corresponding contig containing the gene. I would like to recover the cDNA sequence (i.e. only the exons spliced together) . Are there any tools to do this? tblastn only gives me the fragments, but doesn't piece them together.

Answer 1

[Haibao Tang]

I have used both exonerate and NAP-AAT. They work reasonably well, especially when the protein is derived from this contig.

I would like to mention a recent paper on this by Gotoh, which reviews some of the softwares, as well as his new software.

[Michael Kuhn]

Thanks a lot for the answer! Exonerate works fine for my purposes with this command line: "exonerate --model protein2genome query.fa contig.fna --ryo ">%ti (%tab - %tae)n%tcsn" | less"

Answer 2

[Michael Schubert]

Take a look at Wise2, although it's not a perfect match as it only aligns protein sequences with contigs and you have to extract the exons yourself.

[Yannick Wurm]

is there a modern equivalent to blastwise?

[Michael Schubert]

I have never used blastwise.