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); }