ORF translation of a huge number of DNA sequences
1 answer
You can split your file into smaller files with GenomeTools (gt splitfasta -numfiles 60 seqs.fasta), faSplit by Jim Kent or fasta-splitter by Kirill Kryukov, then loop or parallel through the files.
Log in to answer this question.
More posts like this
-
Software Tool to Translate DNA to Protein
written by sil_bioinfo •
-
How to work with Jmodeltest2 on XSEDE on Cipres
written by ouassal.keltoum •
-
Need help : How can I get the Emboss sixpack output in GFF file
written by shiv •
-
Need Help: How can I perform Emboss-sixpack on multi-sequences fasta file
written by shiv •
-
Alingment of ~23 Mb sequences
written by l.souza
-
Software for cleaning FASTA files
written by l.souza
-
Big Data Translation
written by l.souza
-
Tool for translation of large set of DNA sequences to protein
written by onspotproductions
-
Emboss sixpack replacement
written by onspotproductions
-
Trouble Calling Emboss Program From Python
written by charles.bridges
Do you have to extract the ORF, or your dna sequences are ATG-STOP already? Eitherway, Biopython can make quick work of it. Also, see this post.
I need to extract the ORF! Can I do this with Biopython?
You will need to know a little more about your sequences. Are the ORFs in the forward or reverse direction? If reverse, you'll need to take the reverse complement, and then find 'ATG' in a window search of 3 nucleotides. Do you happen to know the approximate length of the ORFs? You will have to then search for stop codons. With finding 'ATG', you can then translate in the forward direction from the first position.
They are in the forward direction. They have something about 7000nt.
see Emboss sixpack replacement : download & install the standalone version.