另一种blast2go的思路
blast2GO是用序列相似性得到go的信息,而网上很多软件要么针对某个蛋白的,要么就是只能网上运行的,还有就是只有windows的,烦- - 但如果用名字对应会怎么样,尝试了一番记录如下: 我们的序列先进行blast Nr库得到的是比对的蛋白序列,名字有gi号和refseq号; NCBI FTP里面有个gene2go的文件,发现是geneid对应go号; 又从FTP里面发现了gene2refseq文件,是geneid对应gi号; 好了,由此我们可以把这三者进行了串联,结果发现了350多条对应关系,而blast2go的结果有303,\(^o^)/~! 但是我将两者取了交集之后,发现只有5个!只有5个!只有5个!!!瞬间伐开心。。。 所以我在想是不是有几种情况要考虑: 首先是过滤,我们blas...阅读全文