Question: Finding longest non overlapping open reading frames
0
gravatar for Mr.RB
4 weeks ago by
Mr.RB10
Mr.RB10 wrote:

I found predORF from systemPipeR to predict the longest non-overlapping reading frames in a genome, however these are the longest ORFs per strand. I want to have the longest ORF considering both strands, for example:

AAAA**ATG**XXXXXX**TGA**AAAA     CDS = XXXXXX = 6 nt  <- pick this one
AAAATAGX**AGT**XXACT**GTA**A     CDS = XXACT  = 5 nt

Is there any package in R that can do this, I suppose it isn't that easy to code this myself since it requires quite some comparisons when both strands need to be considered.

ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by Mr.RB10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 148 users visited in the last hour