Hi,
We use a custom re-formatting of IPI FASTA headers, similar but
different to that performed by the TPP util script, where:
>IPI:IPI00000001.2|SWISS-PROT:O95793-1|TREMBL:Q59F99 <description>
becomes:
>IPI00000001.2 <description> SWISS-PROT:O95793-1|TREMBL:Q59F99
This causes problems with a recent change to RefreshParser.cxx (svn
rev 4594). When there is no / only one non-IPI accession listed, the
translated header has no | symbol, and the entire header is pulled out
as the protein name, leading to problems.
The following change at line 616 of RefreshParser.cxx fixes this issue
for IPI re-formatted headers which may not contain a | after
reformatting. Just thought it may be useful to anyone who has come up
against a similar problem.
- if (!strncmp(szHdr_, ">IPI", 4) && strchr(szHdr_, '|')
<strchr(szHdr_,' '))
+ if (!strncmp(szHdr_, ">IPI", 4) && strchr(szHdr_, '|')
<strchr(szHdr_,' ') && strchr(szHdr_, '|') != NULL)
DT
|
