Question: Finding longest non overlapping open reading frames
0
gravatar for Mr.RB
6 months ago by
Mr.RB20
Mr.RB20 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 6 months ago • written 6 months ago by Mr.RB20
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: 418 users visited in the last hour