RefreshParser IPI Issue

1 view
Skip to first unread message

dctrud

unread,
Oct 28, 2009, 11:37:20 AM10/28/09
to spctools-discuss
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

Reply all
Reply to author
Forward
0 new messages