blast+安裝及簡單使用


軟件下載

wget https://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/2.9.0/ncbi-blast+-2.9.0-src.tar.gz
tar zxvf ncbi-blast+-2.9.0-src.tar.gz
cd ncbi-blast+-2.9.0-src
./configure
make
make install

BLAST+的一般用法如下:
格式化數據庫

makeblastdb -in db.fasta -dbtype prot -parse_seqids -out dbname

參數說明:
-in:待格式化的序列文件
-dbtype:數據庫類型,prot或nucl
-out:數據庫名
蛋白序列比對蛋白數據庫(blastp)

blastp -query seq.fasta -out seq.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8

參數說明:
-query: 輸入文件路徑及文件名
-out:輸出文件路徑及文件名
-db:格式化了的數據庫路徑及數據庫名
-outfmt:輸出文件格式,總共有12種格式,6是tabular格式對應BLAST的m8格式
-evalue:設置輸出結果的e-value值
-num_descriptions:tabular格式輸出結果的條數
-num_threads:線程數
核酸序列比對核酸數據庫(blastn)以及核酸序列比對蛋白數據庫(blastx)
與上面的blastp用法類似:

blastn -query seq.fasta -out seq.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8
blastx -query seq.fasta -out seq.blast -db dbname -outfmt 6 -evalue 1e-5 -num_descriptions 10 -num_threads 8

參考BGI參數,結果以xml格式輸出時,只輸出了5條比對結果

blastx  -evalue 1e-5  -num_threads  8 -db animal.fa -query All-Unigene.fa.9 -out All-Unigene.fa.9.blast.nr -outfmt 5 -max_target_seqs 5


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM