How to get the description of bl
2018-10-17 本文已影响0人
OmicsAcademy
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
