Missing enzyme

k saward

unread,

Jun 7, 2026, 3:34:50 AMJun 7

to Stacks

Hi,

I want to use process_radtags for my SLAF-Seq raw reads, but my second restriction enzyme is "Hpy166II" which is not included in the process_radtags enzyme list. Would you be able to add it?

Many thanks for considering my request!

k saward

unread,

Jul 31, 2026, 5:44:25 PM (12 days ago) Jul 31

to Stacks

Hi, just trying to bump this thread. Would really appreciate any help with this enzyme. Many thanks :)

Julian Catchen

unread,

Jul 31, 2026, 5:46:16 PM (12 days ago) Jul 31

to Stacks

Hi,

The enzyme check in process_radtags just verifies that the expected cut site overhang is present in the data. Since the Hpy166II only leaves a three nucleotide overhang on the paired-end read, and the third nucleotide is non-specific, not much is gained by checking that it is present. So, you can just omit the second enzyme when you run the program and it will proceed without checking for it and everything should work fine otherwise.