You can use blastdbcmd
to extract FASTA sequences from the BLAST database. Then you can append your extra sequence to the newly created FASTA file, and remake the BLAST DB using makeblastdb
. Assuming you have a protein BLAST database in the current working directory, and a sequence you want to add to it in new.fa
:
$ blastdbcmd -db proteindb -dbtype prot -out db.fa -entry all
$ cat new.fa >> db.fa
$ makeblastdb -in db.fa -out newproteindb -dbtype prot
(Replace prot
with nucl
if you have a nucleotide database)
No comments:
Post a Comment