blast 真的是初心了～
想要啥样的输出格式就能有啥样的输出格式，是不是棒棒哒？
第一步 makeblastdb
这里的db指的是你的subject库，即你从NCBI上下载的序列集，也可以是nr，nt数据库。
-dbtype prot/nucl 蛋白或核酸
makeblastdb -in HarmHassOBPs_prot.fasta -dbtype prot
第二步blast
-outfmt 后面可以直接跟数字
0 = Pairwise,
1 = Query-anchored showing identities,
2 = Query-anchored no identities,
3 = Flat query-anchored showing identities,
4 = Flat query-anchored no identities,
5 = BLAST XML,
6 = Tabular,
7 = Tabular with comment lines,
8 = Seqalign (Text ASN.1),
9 = Seqalign (Binary ASN.1),
10 = Comma-separated values,
11 = BLAST archive (ASN.1),
12 = Seqalign (JSON),
13 = Multiple-file BLAST JSON,
14 = Multiple-file BLAST XML2,
15 = Single-file BLAST JSON,
16 = Single-file BLAST XML2,
18 = Organism Report
Options 6, 7 and 10 can be additionally configured to produce a custom format specified by space delimited format specifiers.
blastx -query Hass.unigenes.fa -db HarmHassOBPs_prot.fasta -evalue 1e-5 -outfmt "6 qseqid stitle sseqid pident length qstart qend sstart send evalue score" -out HassOBP181030.2
