Entering edit mode
7.5 years ago
hivemind
▴
20
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.
I have encountered the same problem, and I have not solved it yet. Have you found a solution