Entering edit mode
9.8 years ago
biolab
★
1.4k
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!
thank you Pierre, the shell command is really cool.