Blast本地化快速上手教程

(*圖片來自NCBI官網(wǎng))https://blast.ncbi.nlm.nih.gov/Blast.cgi



在線提交地址:https://blast.ncbi.nlm.nih.gov/Blast.cgi

官方文檔:https://www.ncbi.nlm.nih.gov/books/NBK279690/

NCBI使用指南


步驟01:下載庫文件

????使用BLAST+自帶的update_blastdb.pl腳本下載nr和nt等庫文件,直接運行下列命令即可自動下載,或者fasta手動格式化數(shù)據(jù)庫。


步驟02:格式化數(shù)據(jù)庫

? ? ? ? 示例:makeblastdb -in demo.fasta -dbtype prot -parse_seqids -out dbname

????????參數(shù)說明:

????????-in:待格式化的序列文件

????????-dbtype:數(shù)據(jù)庫類型,prot或nucl

????????-out:數(shù)據(jù)庫名

步驟03:序列比對

????????蛋白序列比對蛋白數(shù)據(jù)庫(blastp)

????????blastp -query demo.fasta -outdemo.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8

????????參數(shù)說明:

????????-query: 輸入文件路徑及文件名

????????-out:輸出文件路徑及文件名

????????-db:格式化了的數(shù)據(jù)庫路徑及數(shù)據(jù)庫名

????????-outfmt:輸出文件格式,總共有12種格式

????????-evalue:設(shè)置輸出結(jié)果的e-value值

????????-num_descriptions:tabular格式輸出結(jié)果的條數(shù)

????????-num_threads:線程數(shù)

????????核酸序列比對核酸數(shù)據(jù)庫(blastn、blastx)

????????blastn -query demo.fasta -outdemo.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8

????????blastx -query demo.fasta -outdemo.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8

????????參數(shù)說明:

????????-query: 輸入文件路徑及文件名

????????-out:輸出文件路徑及文件名

????????-db:格式化了的數(shù)據(jù)庫路徑及數(shù)據(jù)庫名

????????-outfmt:輸出文件格式,總共有12種格式

????????-evalue:設(shè)置輸出結(jié)果的e-value值

????????-num_descriptions:tabular格式輸出結(jié)果的條數(shù)

????????-num_threads:線程數(shù)

步驟04:參考資料

1.https://www.ncbi.nlm.nih.gov/books/NBK279690/

2.http://blast.ncbi.nlm.nih.gov/

愛上生物大數(shù)據(jù)??公眾號:aisbio

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容