Possible mistake in digest_genome.py

[DavideBrex]

Running the script, If I set the "-r" flag equal to an enzyme name (e.g. hindiii), the program will halt with this message:

Find unexpected character in restriction motif.

I think the problem is here:

HiC-Pro/bin/utils/digest_genome.py

Lines 144 to 149 in 351471b

	for nuc in list(set(cs)):
	if nuc not in ['A','T','G','C','N','^']:
	print("Find unexpected character ['{}']in restriction motif".format(nuc))
	print("Note that multiple motifs should be separated by a space (not a comma !)")
	sys.exit(-1)

Where for nuc in list(set(cs)): should be instead for nuc in list(set(cseq)):

Let me know what you think please! and thank you for the amazing pipeline.
Davide

[nservant]

Hi ! yes, indeed, I think you're right !

[nservant]

I'll fix that in the devel ! many thanks

[nservant]

done !