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
How to get the description of blast hit using blastdbcmd?
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
