gravatar for Mensur Dlakic

8 hours ago by


This is a crude way, but it should work.

Download and unpack Swiss-Prot:


gunzip uniprot_sprot.fasta.gz

Save your IDs into a file, say ID.txt. I did so for the 5 ID numbers you provided. Then simply (f)grep your ID file against the database, and either count the hits or save them in a file.


fgrep -f ID.txt uniprot_sprot.fasta | wc -l

It counts 3 matches.

Printing the output:

fgrep -f ID.txt uniprot_sprot.fasta

>sp|Q8C050|KS6A5_MOUSE Ribosomal protein S6 kinase alpha-5 OS=Mus musculus OX=10090 GN=Rps6ka5 PE=1 SV=2
>sp|P47811|MK14_MOUSE Mitogen-activated protein kinase 14 OS=Mus musculus OX=10090 GN=Mapk14 PE=1 SV=3
>sp|Q15672|TWST1_HUMAN Twist-related protein 1 OS=Homo sapiens OX=9606 GN=TWIST1 PE=1 SV=1

Source link