本地blast步骤

# 本地化基因组数据库构建
mkdir /opt/biosoft/ncbi-blast-2.9.0+/db/
cd /opt/biosoft/ncbi-blast-2.9.0+/db/
# 构建核酸数据库
#makeblastdb -in ~/00.incipient_data/data_for_genome_assembling/assemblies_of_Malassezia_sympodialis/Malassezia_sympodialis.genome_V01.fasta -dbtype nucl -title Malassezia_sympodialis_V01.genome -parse_seqids -out #Malassezia_sympodialis_V01.genome -logfile Malassezia_sympodialis_V01.genome.log
# 构建蛋白质数据库
makeblastdb -in /home/train/Ficus_erecta.protein.fasta -dbtype prot -title Ficus_erecta.protein -parse_seqids -out Ficus_erecta.protein -logfile protein.log

# 本地化运行 blast 程序
mkdir -p ~00.blast_to_genome
cd ~00.blast_to_genome
# 在https://www.ncbi.nlm.nih.gov/gene/数据库搜索"3-phytase a"，点击Send to，选择File，在Format中选择UI List，下载文件保存为gi.ncbi_gene_phyA.list
#cp ~/00.incipient_data/data_for_functional_annotation/gi.ncbi_gene_phyA.list ./
#ncbi_acc_or_gi_to_fasta.pl gi.ncbi_gene_phyA.list 10 > phyA.nucl.fasta
#perl -p -e 's/>(\S+).*/>$1/; s/\|/_/g' ~/00.incipient_data/data_for_functional_annotation/phyA.nucl.fasta > phyA.nucl.fasta
#blastn -query phyA.nucl.fasta -db Malassezia_sympodialis_V01.genome -out blastn.out -evalue 1e-5 -outfmt 7 -num_threads 4
# 在https://www.uniprot.org/uniprot/搜索"3-phytase a" AND reviewed:yes，并下载数据
#perl -p -e 's/>(\S+).*/>$1/; s/\|/_/g' ~/00.incipient_data/data_for_functional_annotation/uniprot-_3-phytase+a_-filtered-reviewed_yes.fasta > phyA.pep.fasta

blastp -query phyA.pep.fasta -db Ficus_erecta.protein -out blastp.out -evalue 1e-5 -outfmt 7 -num_threads 4