另一種blast2go的思路

這是一個建立於 的文章，其中的資訊可能已經有所發展或是發生改變。

blast2GO是用序列相似性得到go的資訊，而網上很多軟體要麼針對某個蛋白的，要麼就是只能網上啟動並執行，還有就是只有windows的，煩- -

但如果用名字對應會怎麼樣，嘗試了一番記錄如下：

我們的序列先進行blast Nr庫得到的是比對的蛋白序列，名字有gi號和refseq號；

NCBI FTP裡面有個gene2go的檔案，發現是geneid對應go號；

又從FTP裡面發現了gene2refseq檔案，是geneid對應gi號；

好了，由此我們可以把這三者進行了串聯，結果發現了350多條對應關係，而blast2go的結果有303，\(^o^)/~！

但是我將兩者取了交集之後，發現只有5個！只有5個！只有5個！！！瞬間伐開心。。。

所以我在想是不是有幾種情況要考慮：

首先是過濾，我們blast的參數都是設定好的，當然可以過濾一部分，程式裡面取的是在資料庫裡面有的且最好的，但blast2go肯定做了另一種過濾；

其次是資料庫，blast2go可能取了不僅僅是NCBI裡面的gene2go，還有其他的，還沒搞明白；

最後是程式，是不是blast2go用了層次關聯排除了一部分，挖出了一部分，我簡直是醉了！！！

還有就是2016年了，做個小總結，雖然有幾個月沒更新（心煩的事情），但最終暫時解決了，明年會怎樣？