[khmer] Unique reverse-complement

C. Titus Brown ctbrown at ucdavis.edu
Wed Dec 7 12:31:57 PST 2016


On Wed, Dec 07, 2016 at 08:18:31PM +0000, Bintu, Bogdan wrote:
> 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.

Hi Bogdan,

thanks for the email!

Can you try doing:

make clean
python setup.py install

after the modification?  That should work.

thanks,
--titus
-- 
C. Titus Brown, ctbrown at ucdavis.edu



More information about the khmer mailing list