Dear community,
I didn't find a direct answer to my question. I am planning to annotate a huge bunge of genomes and metagenomes with profile HMMs like Pfam.
Because of the large size of the data set, I won't be able to optimize individual alignments of the HMMs and I am looking for a good set of parameters for hmmsearch (inclusion parameters) and postprocessing thresholds for e-values and coverage values, i.e. length of the match with respect to the length of the domain model. Could you please cite or tell me infos about finding values that are a good compromise!? And how would it be when processing some metagenomes?
3 answers
The various databases (PFAM,SMART,TIGR,PANTHER) use their own parameters and in some cases several different pre-selecting steps for running the hmms. Have you consider just running [?]InterPro scan[?]? It will let you select which application/databases you want to use, including the HMM databases you've listed.
The only downside is that it's pretty slow. But if you turn off most of the applications you don't need, it shouldn't be that bad.
You also can download the software and run it yourself instead of on their server.
Log in to answer this question.
-
How pfam is made?
written by GenesisBio •
-
hmmsearch vs hmmscan - e-value, speed, & output differences?
written by jabaron.phd •
-
How to parse HMMsearch domtblout, resolve overlaps & return PDA?
written by Anand Rao
-
E-value results from hmmsearch are not accurate. [HMMER]
written by danzanzu •
-
Annotating transposons in bacterial genomes
written by wanderingstefan •
-
Dubious Pfam Domain Annotation Of A Protein
written by Pappu
-
How Can I Efficiently Parse Hmmer Results To Find Listed Domains And Their Architecture?
written by James Ashmore •
-
Realignment Of A Protein Domain After Global Multiple Sequence Alignment
written by miquelduranfrigola
-
Gene Annotation By Protein Domain Match Criteria
written by Burlappsack •
-
How Do I Choose The Proper Data Representation When Designing A Webportal
written by Mariam Rizkallah •
I don't know the answer but I am pretty much looking forward to one. Great question!
Thanks for the support. Let's hope that someone may provide the answers. I wanted to add that I will combine PFAM and SMART HMMs with custom build HMMs and I would like to have some general literature or review on how one should choose thresholds depending on the nature of modelled domains/proteins, without making it too difficult.
Hello, I am facing the exact same problem and I am struggling quite a lot to find any reference with an explicit set of parameter I could use. Did you get any relevant answer for this post? Thanks a lot