BOL: Bash oneliner to extract all ids from a multifasta file

X

BOL wishing you a very and Happy New year

Alternative content

BioScripts
biogeek
Bash oneliner to extract all ids from a multifasta file

Bash oneliner to extract all ids from a multifasta file

By biogeek 1768 days ago

#List of ids - one per line in allIds.txt

$ awk 'BEGIN{while((getline<"allIds.txt")>0)l[">"$1]=1}/^>/{f=!l[$1]}f' seq.fa

# You can play with this f=!l[$1 ] if wanted to extract or not extract the ids