Multiple Sequence Alignment (Same Directory Linux)
1
0
Entering edit mode
4.1 years ago

Hello,

I'm used to doing smaller multiple sequence alignments (MSA) for 10-20 at once .fasta or .faa files, however I have over 100 files stored in the same directory, is there a specific command to do a MSA for all my files present in the same directory? Below I put an example of something I was trying to do but I want an all vs all alignment, not individually. Preferably with mafft or clustalw :)

Thanks in advanced

for i in *.faa; do mafft --quiet $i > ${i%..faa}.aligned.faa; done
msa assembly sequence alignment • 2.3k views
ADD COMMENT
3
Entering edit mode
4.1 years ago
Joe 21k

This isn't how you do MSA, even with a small number of files. You need all the sequences to be in the same file, or you need to provide all the fasta's at once via STDIN (if the tool supports it).

cat *.faa > allseqs.fa
mafft allseqs.fa > output.msa
ADD COMMENT
0
Entering edit mode

Thank you for the clarification, I am used to using mega or geneious for this. Best regards

ADD REPLY

Login before adding your answer.

Traffic: 1953 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6