150 GO注释:blast了nr数据库并通过blast2go获得了GO注释号;也通过blast了swiss-prot数据库获得了GO注释号。这两个结果需要合并吗?

老师,你好,

首先,我在本地将我的序列blast了nr数据库,并导入blast2go软件进行了GO注释,并获得了go注释号。另外,我还将数据在本地比对swiss-prot数据库并进一步获得了GO注释号。

1. 我想在想讲这个结果导入WEGO进行作图,请问上面两份结果需要怎么处理,是需要合并吗,怎么合并?

我是参考一片文献进行操作的,他的描述如下:CDS sequences with SSR were aligned against NCBI non-redundant (nr) and SWISS-PROT protein databases (http://www.uniprot.org), using BLASTx (Altschul et al., 1997) with a cutoff E-value 1E-5. Blast results were imported into Blast2Go (Conesa et al., 2005) for GO term mapping. Mapping results were submitted to WEGO (Ye et al., 2006for GO classification.

他的结果如下:

814 and 745 CDSs with SSR in genomes of C. exilicauda and M. martensii, respectively, were subjected to GO analysis based on sequence alignment and the results were used in GO enrichment analysis. In total, 261 (C. exilicauda) and 197 (M. martensii) CDSs were assigned with 1562 (C. exilicauda) and 1390 (M. martensii) GO terms of known
function. The percentages of CDSs with SSR assigned to each subcategory are shown in
Fig. 4. Comparison of function distribution between these two species demonstrated that only CDSs of C. exilicauda were assigned to auxiliary transport protein (GO: 0015457) and rhythmic process (GO: 0048511). Only CDSs of M. martensii were assigned to virion (GO: 0,019,012).  

attachments-2019-06-Oyd56EB65d09953f02998.jpg

2. 但奇怪的是:我把我从blast了nr数据库获得的GO注释号的几个文件导入WEGO后,发现缺少细胞组分这一块(如果我只放单个文件进去就会有CP,MF和BP),这又是什么问题?

attachments-2019-06-tBL8vr0S5d0995b6b7c07.jpg













请先 登录 后评论

最佳答案 2019-07-04 08:03

这个问题比较复杂,以前没有做过。得检查数据,是不是那一步出错了; 

请先 登录 后评论

其它 1 个回答

Mayomics

请走过路过的朋友帮忙解答下

请先 登录 后评论