建立 DNA blastdb
為了將查詢序列與參考序列進行比較,你必須建立引用的 blastdb。這是使用 makeblastdb
完成的,當安裝爆炸時包含 makeblastdb
。
makeblastdb -in <input fasta> -dbtype nucl -out <label for database>
所以如果你有一個檔案 reference.fasta
包含以下記錄:
>reference_1
ATCGATAAA
>reference_2
ATCGATCCC
你將執行以下內容:
makeblastdb -in reference.fasta -dbtype nucl -out my_database
這將建立以下檔案:
my_database.nhr
my_database.nin
my_database.nsq
注意,資料庫檔案標有 -out
引數。