This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Remove A List of Sequences from the MSA File

Hello Everyone,

I have an MSA file containing hundreds of sequences aligned by BWA. Now I would like to remove a list of sequences from the MSA file, is there a quick way to do this? or Could the MSA file be split into two?!

I would appreciate your response!

bamtools samtools multiple bwa bcftools alignment sequence

1 answer

You can prepare a list with the sequences you want to keep (so the reverse of the sequences you want to remove) and use seqtk. As explained in their github (https://github.com/lh3/seqtk):

 seqtk subseq in.fq name.lst > out.fq

Log in to answer this question.