a44421a79fb36cc2036fe116b97ea3bc9590cd0c braney Fri Dec 2 09:34:39 2011 -0800 removed rcsid (#295) diff --git src/lib/spacedSeed.c src/lib/spacedSeed.c index 1544e27..169e223 100644 --- src/lib/spacedSeed.c +++ src/lib/spacedSeed.c @@ -1,66 +1,65 @@ /* spacedSeed - stuff to help with spaced seeds for alignments. */ #include "common.h" #include "spacedSeed.h" -static char const rcsid[] = "$Id: spacedSeed.c,v 1.1 2005/01/10 00:02:37 kent Exp $"; /* Seeds - the weight 9 and 11 seeds are from PatternHunter paper. * The weights 10,12,13,14,15,16,17 and 18 are from the Choi, Zeng, * and Zhang paper. The others are just guesses. */ char *spacedSeeds[] = { /* 0 */ "", /* 1 */ "1", /* 2 */ "11", /* 3 */ "1101", /* 4 */ "110101", /* 5 */ "1101011", /* 6 */ "111001011", /* 7 */ "1110010111", /* 8 */ "1110010100111", /* 9 */ "111001010011011", /* 10 */ "1101100011010111", /* 11 */ "111010010100110111", /* 12 */ "111010110100110111", /* 13 */ "11101011001100101111", /* 14 */ "111011100101100101111", /* 15 */ "11110010101011001101111", #ifdef EVER_NEEDED_IN_64_BIT_MACHINE /* 16 */ "111100110101011001101111", /* 17 */ "111101010111001101101111", /* 18 */ "1111011001110101011011111", #endif /* EVER_NEEDED_IN_64_BIT_MACHINE */ }; int spacedSeedMaxWeight() /* Return max weight of spaced seed. */ { return ArraySize(spacedSeeds)-1; } int *spacedSeedOffsets(int weight) /* Return array with offsets for seed of given weight. */ { char *seed; int *output, offset, outCount = 0, seedSize; assert(weight >= 1 && weight < ArraySize(spacedSeeds)); seed = spacedSeeds[weight]; seedSize = strlen(seed); AllocArray(output, weight); for (offset=0; offset<seedSize; ++offset) { if (seed[offset] == '1') output[outCount++] = offset; } assert(outCount == weight); return output; } int spacedSeedSpan(int weight) /* Return span of seed of given weight */ { return strlen(spacedSeeds[weight]); }