Perl Script To Retrieve All Orf?
Hello,
I have a draft assembly. Does anyone know of scripts to retrieve all ORF in protein format from a a fasta file of contigs?
Adrian
• 5,155 views
•
link
updated
by
Abhishek •
•
written
by
Adrian Pelin
•
written
by
Adrian Pelin
1 answer
Hi
From EMBOSS toolkit:
getorf -sequence genome.fasta -outseq genome.ORFs -minsize 180 -find 1 &
Cheers!
• 528 views
•
link
written
by
Rahul Sharma
Log in to answer this question.
More posts like this
-
Extracting sequence and annotation of a particular gene
written by Adrian Pelin
-
QUAST alignment contig retrieval
written by bpz •
-
Annotating/Finding "NNNNNN" scaffold gaps in genome assembly
written by Adrian Pelin
-
Importance of the reference genome, it's quality when mapping reads
written by Adrian Pelin
-
Using R to calculate diversity between hundreds of orthologous sequences
written by Adrian Pelin
-
Determining abnormalities in synteny blocks
written by Adrian Pelin
-
bedtools: extracting no coverage regions
written by Adrian Pelin
-
Bwa Map Single And Pe Reads To Ref?
written by Adrian Pelin
-
Cnv With Cufflinks?
written by Adrian Pelin
-
Mate-Pair/Paired-End Contamination
written by Adrian Pelin
Unless this is a prokaryote, getting all open reading frames from a draft assembly is not an informative analysis (splicing, low gene density), and also for bacteria it is of very limited use. Instead, look for gene prediction , e.g. on BioStar: gene-prediction
I work on microsporidia. They have very little introns (up to 20 genes with introns, some none at all) and small genomes. They are Eukaryotes.
Then you can use getorf as suggested by R@hul, you should still attempt to do a proper gene prediction.
Please share the fully functional perl script to translate CDNA to ORF (protein) selecting the longest one only. I have Active Perl installed.