This program aligns proteins from the open reading frame to SwissProt. This is done using BLASTP.
This program runs hmmscan to find protein domains. This takes the protein sequeces found in the open reading frame's and searches those against the profile-HMM database.
This program uses the 'TransDecoder. Predict' program to predict the likely coding regions from the open reading frame's that were identifies from Transdecoder.LongOrfs. The homology searhes will aslo be incuded in the open reading frame retention criteria.
This program aligns the output of TransDecoder.Predict against SwissProt, and sets a filter for evaluate < 1e-10 with the output in tabular format.
- bash
Victoria R. Liebsch
7 November 2021