This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CDS sequence from mRNA and protein sequence

Hello guys,

Here's for the (probably) dumb question of the day, but I can't seem to find fitting solutions in a reasonable amount of time, so I might be missing something here.

I have an mRNA multifasta and a protein multifasta (IDs are correspondent between the two). What I'm after, is to get the CDS sequences from the mRNA file, using the proteins as model for that. Or any other way to get the CDS sequences from those mRNAs. I do not have genome coordinates or any type of GFF file.

Any easy way that you might know?

fasta cds protein mrna code

You can try blastx (with the protein fasta as db), I assume your mRNA sequences are without introns?

BLASTX should do the trick but not sure if you would be able to extract the exact coding sequences using it, however, you can certainly use the coordinates.

1 answer

Simply run some software on it get the longest ORF (plenty of tools around for that) from the mRNA, should (in theory) end up to the proteins. Additionally just double-check them against your protein set (eg calculate MD5 key and compare both).

potential software to use:

and many other around I guess

I am interested in which software would be the way to go!

OK, I'll add some to my 'answer'

My general advise would be to go for one that does nothing more than to look for the longest orf, so preferentially none that use coding potential or such but simply (and bluntly) gets the longest one from a potential start to a stop (by limiting to the forward strand even)

Log in to answer this question.