[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