SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Using profile HMMs (HMMER) to annotate genomes teatime Bioinformatics 0 12-07-2011 01:34 PM

Reply
 
Thread Tools
Old 06-10-2013, 07:50 AM   #1
rhinoceros
Senior Member
 
Location: sub-surface moon base

Join Date: Apr 2013
Posts: 372
Default hmmer and pfam

Could somebody explain to me briefly how to use Pfam's hmm with hmmer.

For example, what are the steps for using the file from this site (raw hmm) with hmmalign?

hmmalign -o output thefile.hmm myseqs.fasta

Doesn't produce any errors. However, it's no alignment. All the sequences are still in the same exact order as they were in myseqs.fasta
rhinoceros is offline   Reply With Quote
Old 06-11-2013, 03:58 AM   #2
Apexy
Member
 
Location: Africa

Join Date: Apr 2011
Posts: 62
Default

Hi,
I noticed your query has note attended to. Let me attempt a learners non-specialist answer. if I remember well, there are a number of executables in the hmmer package to achieve different things. What do you plan to achieve with hmmer?
For example, the only thing I remember is to acribe functional lable to sequences in a FASTA file. I evoked hmmer like so:
hmmscan --domtblout test.tab --cut_ga --acc --cpu 10 Pfam-A.hmm test.fa

where test.tab is the output in tabular format.

Hope this helps you one step into your task.

cheers
Apexy is offline   Reply With Quote
Old 06-11-2013, 05:44 AM   #3
rhinoceros
Senior Member
 
Location: sub-surface moon base

Join Date: Apr 2013
Posts: 372
Default

I think I figured it out. There's nothing wrong with my output. hmmalign doesn't align sequences (vs sequences), more like, it aligns sequences to a hmm. Thus, the order of the sequences in the alignment is meaningless (and subsequently the same than in the input file).
rhinoceros is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:40 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO