Biostar Beta. Not for public use.
getorf: Remove Overlapping Open Reading Frames
1
Entering edit mode
2.5 years ago
JMac • 20
Ireland

Hi,

I am using getorf to extract open reading frames (at least 70 amino acids long) from genomes/scaffolds. I am using the command:

getorf -find 1 -minsize 210 -sequence genome.fasta -outseq orfs.fasta

My question is, does getorf report overlapping ORFs within the same reading frame?

Out of a set of alternative overlapping ORFs within the same reading frame, we only wish to extract the longest ORF.

Thanks for any help

ADD COMMENTlink
0
Entering edit mode

I have read the manual but it doesn't mention anything about alternative/overlapping ORFs.

ADD REPLYlink
0
Entering edit mode

what do your ids look like in your fasta? If they have similar patterns, such as different isoforms from Trinity for example, where you want to the longest isoform per gene cluster, I have a scrip that may help.

ADD REPLYlink
0
Entering edit mode
18 months ago
utsafar • 20

I used getorf in galaxy. It reports all ORFs in all reading frames even in complementary strand. For my case I extracted all ORFs of desired contigs and then, using microsoft excel deleted all ORFs of each contig except the longest one.

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3