Question: GO annotation on diamond blasted transcriptome (nr db)
2
gravatar for jfo
10 months ago by
jfo40
jfo40 wrote:

Hello!

So far, Approaches to GO annotation of transcriptomes of non-model organisms (besides B2G) point to the use of Trinotate and Genescf. I would like to ask if there is an easy way to directly add GO annotation to a diamond blast result (against NCBI nr)? Since B2G is really slow on a basic subscription, I would like to find a way to do GO annotation.

Also, here's my confusion right now: can I just use the GO annotation from Trinotate pipeline or do I really need to use the nr result for this since more sequences have annotations in nr? If using only trinotate (swissprot+eggnog+pfam), would I be missing some GO terms that may be recovered if I've used nr? Or will it be safe to just use interproscan GO result to represent the transcriptome instead?

Finally, what is the best way to convert the genbank protein IDs (XX#####.#) and RefSeq IDs (XX_#*) from diamond blast of a non-model organism, to say Gene name (e.g. FOLH, DUPD1)?

Thanks.

rna-seq blast go diamond • 321 views
ADD COMMENTlink written 10 months ago by jfo40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1857 users visited in the last hour