makeblastdb -in PlantTFDB-all_TF_pep.fas -dbtype prot -title PlantTFDB-all_TF_pep.fas
#提取基因列表对应的蛋白序列：
python $scriptdir/get_fa_by_id.py -i ../7.all_enrich/gene.list \
-f ../7.all_enrich/pep.fa -p query
#转录因子blastp注释
blastp -query query.fa -db Ath_pep.fas -max_hsps 1 \
-num_alignments 1 -evalue 1e-5 -num_threads 2 \
-outfmt "6 qseqid sseqid qlen slen pident length mismatch gapopen qstart qend sstart send evalue bitscore qcovs stitle"\
-out tf.blast.out
