[Bioperl-l] Bio::Restriction::Enzyme

Nathan S. Haigh n.haigh at sheffield.ac.uk
Fri Oct 27 14:19:02 UTC 2006


Brian Osborne wrote:
> Nathan,
>
> I don't know how this is supposed to work, there would be different ways to
> make is_prototype true. One way would be to make the enzyme with the first
> occurrence of a given restriction site the prototype (and the next enzymes
> with the same site are isoschizomers). Or, one could wait until one site had
> appeared twice, with 2 different enzymes, then make the first the prototype,
> etc. I would have done it the first way myself but I took a quick look at
> IO/withrefm.pm and it looks like it's doing it the second way. That means
> one can read an enzyme file and end up with no duplicated restriction sites,
> or prototypes and isoschizomers.
>
> Brian O.
>
>   
Hmm, I'd have done it the first way also. Doing it the second way would
mean you only ended up with something as a prototype if there were
multiple enzymes with the same restriction site - is that correct
biologically?

Nath



More information about the Bioperl-l mailing list