[khmer] Unique reverse-complement
Bintu, Bogdan
bbintu at fas.harvard.edu
Wed Dec 7 12:18:31 PST 2016
Dear Dr. Brown,
I am a graduate student using your khmer package.
I am trying to construct a khmer.Countgraph which is not degenerate between a sequence and its reverse complement.
More explicitly, in python when running:
import Khmer
ht = khmer.Countgraph(17,10,1)
ht.consume('GGCTTGCTCCAGCCATG')
ht.get('CATGGCTGGAGCAAGCC')
I would like the last line to return 0.
As suggested in a different thread I tired modifying NO_UNIQUE_RC=0 to NO_UNIQUE_RC=1 in lib/Makefile and then run python setup.py install.
But this did not solve my issue.
Any help would be appreciated.
Kind regards,
Bogdan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.idyll.org/pipermail/khmer/attachments/20161207/15f8de6c/attachment.htm>
More information about the khmer
mailing list