Assigning GO ids to blast hits

Question

[Mousumi Akter]

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

[gv]

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

[Mousumi Akter]

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

[gv]

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