EddyRivasLab / hmmer

HMMER: biological sequence analysis using profile HMMs
http://hmmer.org
Other
317 stars 70 forks source link

hmmconvert can't read H2 "Nucleic" save files #90

Closed traviswheeler closed 7 years ago

traviswheeler commented 7 years ago

ID h80 TITLE hmmconvert can't read H2 "Nucleic" save files AFFECTS 3.0 and earlier. FIXED_IN 3.1b1 STATUS CLOSED XREF J6/61 REPORTED_BY Martin Hartmann martinha@interchange.ubc.ca OPENED_DATE SRE, Tue May 18 14:16:19 2010 CLOSED_DATE SRE, Tue May 18 14:16:20 2010 DESCRIPTION
H2 save files use the tag "nucleic" on their ALPH line. Modern Easel/H3 expects "DNA" or "RNA". p7_hmmfile.c::read_asc20hmm() was passing the token to esl_abc_EncodeType(), which then failed to recognize "nucleic".

Fix is for read_asc20hmm() to decode the tag itself: "nucleic" mapped to eslDNA, "amino" mapped to eslAMINO.