[Bioperl-l] seq_word and pattern counts
torsten.seemann at infotech.monash.edu.au
Tue Feb 28 21:45:16 UTC 2006
> Does anyone know if Bio::Tools::SeqWords
> will do DNA pattern searches and honor ambiguity symbols
> like exist in some restriction enzyme pattern definitions,
> e.g. GGnnCC*
Examination of the code
suggests that all it does is count N-mers of any set of letters,
and does so in a case-insensitive way ie. CAT, Cat, cat are counted as
the same N-mer.
So no it does not handle ambiguity symbols in any special manner.
What would you like it to do?
If a N-mer has 1 "N" in it, does it count towards the 4 possible N-mers
it could be?
If it has 2 "N"s in it, does it count toward all 16 possible
And so on?
Victorian Bioinformatics Consortium, Monash University, Australia
Phone: +61 3 9905 9010
More information about the Bioperl-l