BOL

Extract fasta with Ids. $ perl -ne 'if(/^>(\S+)/){$c=grep{/^$1$/}qw(id1 id2)}print if $c' fasta.file #extract #fasta #sequence #ids