I want to blast several proteomes against CAZy database on my terminal, and I downloaded the CAZy database from dbCAN2:
And here is my blast code after building my own database:
$blastp -query input.fasta -db cazydatabase.fa -evalue 1e-5 -outfmt 6 -out output -max_target_seqs 1
I got over 900 sequences in my result with 1e-5. However, other people's work showed 300-400 is the appropriate significant sequences number. Could someone please help me solve this problem?
Thanks in advance!