管道从宏基因组数据集中发现同工的同系物

This pipeline uses a combination of blast to overcome difficulties with current functional gene annotation techniques. * Uses no E-value cutoff to determine if something is a homolog, to allow for multiple proteins to be compared without difficulty of varying conservancy * Creates a HMM and consensus sequence of a whole gene, rather than singular conserved domains (e.g. PFAM) * Runs against nucleotide databases, to avoid issues of orf-calling, but using protein sequences to for better resolution of matches * Verifies putative homologs against the STRING database to check if the target and the query are in the same COG.

FileUtility.py
HomologFinder.conf
HomologFinder.py
HomologFinderUtility.py
STRINGUtility.py
datamodel.py
.
._gene.list.txt
._db.list.txt
db.list.txt
gene.list.txt
logging.conf
multi.cog.proteins