您好,登錄后才能下訂單哦!
今天小編給大家分享一下blast如何批量注釋輸出結(jié)果指定需要的列信息的相關(guān)知識點,內(nèi)容詳細,邏輯清晰,相信大部分人都還太了解這方面的知識,所以分享這篇文章給大家參考一下,希望大家閱讀完這篇文章后有所收獲,下面我們一起來了解一下吧。
blast批量注釋輸出結(jié)果指定需要的列信息
makeblastdb -in PlantTFDB-all_TF_pep.fas -dbtype prot -title PlantTFDB-all_TF_pep.fas #提取基因列表對應(yīng)的蛋白序列: python $scriptdir/get_fa_by_id.py -i ../7.all_enrich/gene.list \ -f ../7.all_enrich/pep.fa -p query #轉(zhuǎn)錄因子blastp注釋 blastp -query query.fa -db Ath_pep.fas -max_hsps 1 \ -num_alignments 1 -evalue 1e-5 -num_threads 2 \ -outfmt "6 qseqid sseqid qlen slen pident length mismatch gapopen qstart qend sstart send evalue bitscore qcovs stitle"\ -out tf.blast.out
以上就是“blast如何批量注釋輸出結(jié)果指定需要的列信息”這篇文章的所有內(nèi)容,感謝各位的閱讀!相信大家閱讀完這篇文章都有很大的收獲,小編每天都會為大家更新不同的知識,如果還想學(xué)習(xí)更多的知識,請關(guān)注億速云行業(yè)資訊頻道。
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點不代表本網(wǎng)站立場,如果涉及侵權(quán)請聯(lián)系站長郵箱:is@yisu.com進行舉報,并提供相關(guān)證據(jù),一經(jīng)查實,將立刻刪除涉嫌侵權(quán)內(nèi)容。