dctrud
unread,Oct 28, 2009, 11:37:20 AM10/28/09Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
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