This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Hmmer Alignment input parse error

Hello everybody,

I want to build a pipeline for automatically composing MSAs and afterwards creating HMMER Profiles. Now I encounter an error I can't find a solution for:

 # hmmbuild :: profile HMM construction from multiple sequence alignments  
 # HMMER 3.1b2 (February 2015); http://hmmer.org/  
 # Copyright (C) 2015 Howard Hughes Medical Institute.  
 # Freely distributed under the GNU General Public License (GPLv3).  
 # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -  
 # input alignment file:             Q30154.sto  
 # output HMM file:                  test.hmm  
 # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -  

 # idx name                  nseq  alen  mlen eff_nseq re/pos description  
 #---- -------------------- ----- ----- ----- -------- ------ -----------  
 Alignment input parse error:  
 sequence has more than one #=GS <seqname> DE accession line
 while reading Stockholm file Q30154.sto
 at or near line 11

Has anybody encountered this problem and found a solution? I was suspecting that something with the MSA wasn't alright, but it looks in my opinion like the other MSAs I created and which work flawlessly.

alignment hmmer

I have encountered the same problem, and I have not solved it yet. Have you found a solution

0 answers

No answers yet.

Log in to answer this question.