建立 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 引數。