View Single Post
Old 08-13-2013, 08:53 AM   #2
rhinoceros
Senior Member
 
Location: sub-surface moon base

Join Date: Apr 2013
Posts: 372
Default

If your sequences aren't split to multiple lines you can do this with grep. I think:

grep -A 1 -f yourSeqIDFile.txt yourFastaFile.fasta > SeqsFromIDList.fasta
grep -A 1 -v -f yourSeqIDFile.txt yourFastaFile.fasta > TheOtherSeqs.fasta

might remember wrong..


If you have QIIME, you can do this with filter_fasta.py..

Last edited by rhinoceros; 08-13-2013 at 08:56 AM.
rhinoceros is offline   Reply With Quote