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?
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點，簡書系信息發(fā)布平臺，僅提供信息存儲服務(wù)。
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點，簡書系信息發(fā)布平臺，僅提供信息存儲服務(wù)。
