blast

blast是生物信息分析中比較常用的軟件.

sequenceserver

最近想搭建一個局部的web界面的blast,搜到了sequenceserver.
軟件有一些bug,但是它會讓你熟悉blast的整個流程,有參考價值.

blastn

查詢序列和數(shù)據(jù)庫序列相同,為:

>a
AGCTAGCT
>b
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
>c
AATAAGCACCAGCAGTTGAAATACAGTCGGTTGAAAACTTGCCGGGTCGGGTGAAATGCG

blastn命令

makeblastdb -in $fa -parse_seqids -dbtype nucl
blastn -db $db -query $fa -evalue 1e-5 -outfmt 6 -max_target_seqs 1 -num_threads 1 > $out

結(jié)果為:

c c 100.00 60 0 0 1 60 1 60 2e-30 111

只有c序列比對上了;為什么a,b序列都沒有比對上?

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

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容