This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Assigning GO ids to blast hits

I have blasted proteome against uniprot_swissprot proteins that gave a list of blast hits. Now, I want to use Uniprot GOs (goa_uniprot_gcrp.gaf.gz) to find GO of the blast hits and assign them to the protein. Can someone help me with the script In Linux or R?

Thanks

linux go blast
• 1,276 views • link
updated by jeevii • jeevii • written by Mousumi Akter • Mousumi Akter

Hello,

One suggestion, you can use the [uniprot IDs] protein IDs present in your blast hits and upload those IDs in uniprot and get the desired columns of annotation [protein names,Gene names,GO,pathway.lineage,etc...]

For better clarity,if possible could you please post a blast hit columns?

Hope it helps!

Thanks

• 0 views • link
written by jeevii • jeevii

Here are the first two rows with all columns

A1_oases_trim_CL10016Contig1.p1 sp|Q7PC79|XPOT_ARATH 27.451 255 152 5 1 255 409 630 7.41e-21 91.7 A1_oases_trim_CL10016Contig1.p1 sp|Q8H3A7|XPOT_ORYSJ 28.077 260 145 6 1 255 401 623 8.94e-21 91.7

• 0 views • link
written by Mousumi Akter • Mousumi Akter

Hello,

yes from the given entries, extract uniprot IDs Q7PC79 and Q8H3A7 [present in column 2] and submit in uniprot DB to get annotation details

• 0 views • link
written by jeevii • jeevii

0 answers

No answers yet.

Log in to answer this question.

More posts like this