void fixPrimers(struct mahoney *el)
/* Filter out bogus looking primers. */
{
if (!isDna(el->fPrimer, strlen(el->fPrimer)))
{
static struct hash *uniqHash = NULL;
if (uniqHash == NULL)
uniqHash = hashNew(0);
if (!hashLookup(uniqHash, el->fPrimer))
{
hashAdd(uniqHash, el->fPrimer, NULL);
verbose(1, "Suppressing 'primer' %s\n", el->fPrimer);
el->fPrimer = cloneString("");
}
}
if (!isDna(el->rPrimer, strlen(el->rPrimer)))
errAbort("Not a primer: %s", el->rPrimer);
}
boolean seqIsDna(bioSeq *seq)
/* Make educated guess whether sequence is DNA or protein. */
{
return isDna(seq->dna, seq->size);
}
