Dear all,
I have many pairs of protein sequences, an example shown below.
MDDDIAALVVDNGSGMCKAG DDDIAALVDNSGSMCKAG
TTAEREIVRDIKEKLCY TTAEIVRKEKLCYVA
RMQKEITAPSTMKIKI KEITALPSTMKIKII
... ...
I need to perform blastp to compare the two sequences in each pair. What's command for that?
To my knowledge, blastn -subject -query can be used for two nucleotide sequence comparison. So how to compare pr sequences?
THANK YOU very much!