How to annotate Diamond output files
0
0
Entering edit mode
4.0 years ago

I have done the aligning of my protein sequences and of short reads against a protein reference database using diamond. However I would like to annotate my alignments but I would prefer to do this on bash. I found this https://resources.qiagenbioinformatics.com/manuals/clcmgm/current/index.php?manual=Annotate_with_DIAMOND.html to be useful but this is a seperate software.

Question: Is it is possible to annotate this alignment created from Diamond within bash?

alignment • 1.7k views
ADD COMMENT
0
Entering edit mode

Hi,

First of all, what do you mean with "annotate". If you are simply talking about adding names/descriptions to your alignments, you can easily do it with grep, awk, sed, etc, depending on what format of your outputs (e.g. if you got any tabular file of your input/query IDs versus your hits/subject IDs). Also, the proteins database is a well documented/curated one (like nr, refseq, swissprot) or is it a cusomized set of yours?

Cheers

ADD REPLY

Login before adding your answer.

Traffic: 1798 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6