【原创】blast本地化 2018-10-30

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