This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to obtain Coding sequence from Genomic sequence?

Hi Guys,

I have some assembled genomic sequences and I know the exact frames I need to translate to obtain coding sequence and putative protein sequences, but I do not know exactly where to chop the translated frame and hence the coding region and splice sites for intron and exon boundaries. Is there a way to get coding region from given frames? Please share your knowledge.

Thank you!

coding genome
• 3,065 views • link
updated by Ram Ram • written by MAPK MAPK

Please show what you have; if you know the exact frame, then that's where to cut!

• 0 views • link
written by karl.stamm karl.stamm

Thank you for your reply, Karl. When I say assembled the genomic sequences are still unannotated and there are multiple scaffolds. Suppose I have three scaffolds for any particular protein coding gene and if I have to translate the scaffolds in two different frames each. After merging all the translated frames, I will have an unreasonably long peptide sequence/coding region. That is because the frame gets translated further beyond the 'GT' and 'AG' boundary and hence includes the non-coding regions as well. In this case I think I need to have transcriptomic sequences to unambiguously infer the coding regions. Please clarify if otherwise. Thanks again!

• 2 views • link
updated by Ram Ram • written by MAPK MAPK

1 answer

There is an easy ways to do it .

It's Coding Potential Calculator (CPC).

It takes input as fasta sequence.

Once you submit your fasta file, you get CP score for each sequences. >1 is potential protein coding and < -1 is potential non-coding RNA

CPC can reliably discriminate the coding and non-coding transcripts in ~98% accuracy.

• 0 views • link
updated by Ram Ram • written by Manvendra Singh Manvendra Singh

Log in to answer this question.

More posts like this