How to get the description of blast hit using blastdbcmd?

weixin_34279579發表於2018-10-17
blastn -db nt -remote -query query.fa -entrez_query "Homo sapiens [Organism]" -outfmt 6 -out blast_res_2_nt.table
cat blast_res_2_nt.table |perl -e 'while(<>){@ele=split; print "$ele[0]\t$ele[1]\n" if !exists $rec_id{$ele[0]};$rec_id{$ele[0]}++}'  |cut -f2 |sort -u  |blastdbcmd  -db nt  -outfmt "%a|%g|%t" -entry_batch -  -out  blast_res_2_nt.tab_description.tab  >> blast_res_2_nt.tab_description.log 2>&1

相關文章