This is a test version of Biostars. For the public version, visit https://www.biostars.org.
getorf: Remove Overlapping Open Reading Frames

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

genome software orfs genomics
• 3,354 views • link
written by JMac • JMac

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

• 0 views • link
written by JMac • JMac

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.

• 0 views • link
written by st.ph.n st.ph.n

1 answer

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.

• 0 views • link
written by utsafar utsafar

Log in to answer this question.

More posts like this