This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Concatenating single gene alignments

I have run Panaroo on 160 genomes and as one of the outputs I have received single gene alignments files for each gene in all of the genomes they are present in. I wish to concatenate several of these files in a way such that within the concatenated files, sequences which originated from the same genome are concatenated together. What is a simple way of achieving this using R or the linux terminal?

panaroo multiple sequnce alignment msa r

2 answers

I like to use the pxcat program from the phyx tool suite for MSA concatenation:

https://github.com/FePhyFoFum/phyx

You could also use PhyKit, though I haven't tried it yet:

https://jlsteenwyk.com/PhyKIT/usage/index.html#create-concatenation-matrix

This is what I use:

Log in to answer this question.