Question: Searching FASTA file for specific IDs

[aw90]

Hi all,

Apologies if this is answered elsewhere, I've searched for it but can't find!

I have a list of UniProt IDs (e.g., Q95T64 Q8INK6 Q9GNK5) and want to search through them all (n=~450) to see if they exist in a fasta file of a species, outputting two lists, one list of the IDs that do appear in the species fasta, and one list of those that do not!

Grep perhaps?

Many thanks

Andy

EDIT: The list of IDs are contained within a .txt file, return seperated

[rhinoceros]

Something like:

grep '>' species.fasta | grep -f idFile.txt > HeaderLinesThatMatch.txt
grep '>' species.fasta | grep -v -f idFile.txt > HeaderLinesThatDoNotMatch.txt

Might work, might not. I'm too lazy to check how the second command behaves..