blast如何批量注释输出结果指定需要的列信息

今天小编给大家分享一下blast如何批量注释输出结果指定需要的列信息的相关知识点，内容详细，逻辑清晰，相信大部分人都还太了解这方面的知识，所以分享这篇文章给大家参考一下，希望大家阅读完这篇文章后有所收获，下面我们一起来了解一下吧。

blast批量注释输出结果指定需要的列信息

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

以上就是“blast如何批量注释输出结果指定需要的列信息”这篇文章的所有内容，感谢各位的阅读！相信大家阅读完这篇文章都有很大的收获，小编每天都会为大家更新不同的知识，如果还想学习更多的知识，请关注亿速云行业资讯频道。