Remove Noncoding Regions From Cdna Sequences
1
1
Entering edit mode
9.9 years ago
Katy ▴ 10

I am trying to create a cDNA file to calculate Ka/Ks values using PAML. and I am wondering if any of you out there can introduce an automated, quick and reliable method of removing the noncoding regions! clustal does not like stop codons! I appreciate your help :) K

cdna utr alignment • 2.5k views
ADD COMMENT
0
Entering edit mode

Have you got a template protein for the cDNA?

ADD REPLY
0
Entering edit mode

Where do your cDNA sequences come from? your own sequencing? downloaded from public databases? do you have annotations for these sequences? and which species are they from? If you answer these questions, we might be able to help.

ADD REPLY
0
Entering edit mode

Hi, I don't have a template protein for it...but I can use softwares to translate the sequence. This is how I generated my cDNA: I have Brassica nigra BACS, sequenced them, ran a BLAST analysis against Arabidopsis genome and predicted the genes using GlimmerHMM and training data from arabidopsis...then extracted the mRNA features from the gff file. Is it more clear now?

ADD REPLY
0
Entering edit mode

you're stuck at the extraction mRNA step or already got them?

ADD REPLY
0
Entering edit mode

I already have them! need to remove the noncoding region away!

ADD REPLY
0
Entering edit mode

I have the mRNA sequences for each gene...now i need to remove the noncoding sequences!

ADD REPLY
0
Entering edit mode

So you want every sequence to start with ATG and cut off the stop codon? I think you probably need coordinates for the UTRs and stop codons, then use a perl script to feed in the coordinates and sequences to do the trimming.

ADD REPLY
1
Entering edit mode
9.9 years ago
Vitis ★ 2.4k

Ensembl API to grab only the coding regions? You'll have to write some Perl code to do that. I'm not a good programmer but I'm sure the computational folks here will give you an answer with example code.

ADD COMMENT

Login before adding your answer.

Traffic: 2298 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6