How to add function descriptions to FASTA sequences

Short descriptions in fasta sequence help us to quickly gain insight into important information about a sequence. Automatic assignment of Human Readable Descriptions (AHRD) is popularly used to add select descriptions and Gene Ontology terms that are concise, informative and precise. This bash script run DIAMOND to search homology to three different databases (TAIR, uniprot_sprot, and uniprot_trembl), then run AHRD and, finally perform the addition of new information fasta description.

USES 


bash run_AHRD.sh fasta_input
Extract Part of a FASTA Sequences with Position by python script HERE

Dependencies


  • Install AHRD
  • Install DIAMOND (move to dist directory)
  • Download database sequences for TAIR, uniprot_sprot, and uniprot_trembl. make the DIAMOND blast database and name them uniprot_sprot, arabidopsis, and uniprot_trembl and move the file to a new directory database.
  • Download and decompress the resources.tar in the working directory.

run_AHRD.sh


No comments:

Post a Comment

Have Problem ?? Drop a comments here!