cdhit90_muscle_trim_round2-remove-gap.fa为cdhit之后的fasta,也可以直接使用下载的fasta序列，无需cdhit,生成cd90_trim_remove_gap一系列文件

unbuntu@unbuntu-virtual-machine:~/1funclib/1gsk$ makeblastdb -in cdhit90_muscle_trim_round2-remove-gap.fa -dbtype prot -out cd90_trim_remove_gap

生成PSSM矩阵

unbuntu@unbuntu-virtual-machine:~/1funclib/1gsk$ psiblast -query 1gsk.fasta -db cd90_trim_remove_gap -out out_trim -num_iterations 3 -out_ascii_pssm cd90_trim_remove_gap_pssm