This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Multiple Sequence Alignment (Same Directory Linux)

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

1 answer

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

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

Log in to answer this question.