nr注释

blastp -db ~/super/nr-blast/nr -query ~/geneomic/BRAKER-master/new/BRAKER-master/doc/braker.aa -outfmt 5 -num_alignments 5 -max_hsps 1 -out lh_nr.xml
python ~/super/xml2txt.py -i lh_nr.xml -o lh_nr.csv
##保留最佳hit
python ../blast_besthit.py -bin lh_trembl.txt

nt注释

blastn -db nt -query ~/geneomic/BRAKER-master/new/BRAKER-master/doc/braker.codingseq -outfmt 5 -num_alignments 5 -max_hsps 1 -out lh_nt.xml
python ~/super/xml2txt.py -i lh_nt.xml -o lh_nt.csv

kegg注释

https://www.genome.jp/tools/kaas/ 网页注释得到ko号，下载 q00001.keg文件

##python3
python ~/super/kegg/kegg_trans.py q00001.keg gene.list.txt gene_anno.txt pathway_anno.txt

gene.list.txt文件获取

paste <(grep '^>' braker.aa | sed 's/>//' | sort) <(grep '^>' braker.codingseq | sed 's/>//' | sort) > gene.list.txt

gene_anno.txt pathway_anno.txt为输出文件，可以不改动

参考：https://www.jianshu.com/p/48716fa7321b

interproscan注释

本地构建

axel -n 50  #网页地址

百度网盘下载

网页版打开，选择Aria下载