[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