Hello everyone,
I have ~10 fasta files and each file contain more than 50 sequences and I want to get the information about 6-frame translation and ORFs for each sequences, I found Emboss-sixpack to do this work for me but when I went through the manual I got to know that it takes only single sequence as input file. Can you please suggest me with other options of this (may be I've missed) and is there any alternate option to do this thing without splitting the files ..
Thanks in advance
•
written
by
shiv •
This could be a workflow combining your favorite solution of:
How To Split A Multiple Fasta or this aw(k)esome code by Pierre: A: Is there a way to split single .txt file with multiple fasta sequences into indi
Bash Loop For Job Submission and here: A: Bash Loop For Job Submission (needed to be fully parameterized, otherwise sixpack asks for user input, and it doesn't seem to read or write to STDIN/STDOUT)
This should work without having to install any additional software on linux and mac
awk '/^>/ {if(x>0) close(outname); x++; outname=sprintf("_%d.fa",x); print > outname;next;} {if(x>0) print >> outname;}' *.fasta
for f in _*.fa
do
sixpack -sequence $f -outfile $f.sixpack.out -outseq $f.sixpack.fa
done
If you want to have the output in a single file, use cat to combine them.
I think the last line of this answer might be what you want to do. Here's an example to concatenate files:
cat *.fasta > bigfasta.fasta
Also, I'm pretty sure the EMBOSS suite is in a public Galaxy server out there for easy access.
Hi,
Thanks, Michael.. I am trying to implement your suggestion (this awk command) with Python, but I am getting syntax error every time, Can you please help me in this. Below is code :
import subprocess
cmd = "awk '/^>/ {if(x>0) close(outname); x++; outname=sprintf("_%d.faa",x); print > outname;next;} {if(x>0) print >> outname;}' {}".format(FASTA_FILE_PATH)
subprocess.call(cmd, shell=True)
Why call awk from python? Anyway you have to adapt the quotes. E.g. escape them or use alternative quotes. Like qq in Perl. Don’t know if python has similar functionality
Of course it has something similar... https://stackoverflow.com/questions/29559905/does-python-have-an-equivalent-of-perls-qq
but simply not as powerful 💪
Log in to answer this question.
