Skip to content

Latest commit

 

History

History
49 lines (31 loc) · 1.46 KB

README.md

File metadata and controls

49 lines (31 loc) · 1.46 KB
Raw

VERY DEPRECATED

Use eggNOG emapper.py scripts instead

This is a process for sorting a whole proteome into Hmm profiles using hmmscan. outputs:

  1. A flat file of the best hit protein - HMM profile matches with annotations (tophit)
  2. A flat file of all protein - HMM profile matches with annotations (annotated)
  3. A flat file of all protein - HMM matches (all)
  4. A flat file of proteins with no HMM hit (nohit)
  5. The input fasta annotated with the EggNOG hmm group annotations. This is output to the proteomes/[species] folder

Instructions 1.Place the annotation and hmm files for a phylogenetic level in hmms/

ex. euNOG_hmm.tar.gz euNOG.annotations.tsv.gz

HMM profiles come from http://eggnogdb.embl.de/#/app/downloads

2.From the main directory run: bash masterscripts/startPress.sh [level]

ex. bash masterscripts/startPress.sh euNOG

This step takes about 10 TACC minutes

  1. Make a directory for the species that you want to run in proteomes/ ex. mkdir proteomes/arath

  2. Place the species' fasta in its folder in proteomes/ ex. proteomes/arath/uniprot-proteome%3AUP000006548.fasta

  3. After the hmms are pressed, from the main directory run: bash masterscripts/startHmmscan.sh [species] [proteome] [level]

ex. bash master_scripts/startHmmscan.sh arath proteomes/arath/uniprot-proteome%3AUP000006548.fasta euNOG

This step takes up to 20 TACC hours depending on proteome/hmm profile count

tophit + nonhits are combined to create look ups for the othology mass spec analysis