atifrahman / HAWK

Hitting associations with k-mers

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Using a different k-mer size

Asutu opened this issue · comments

commented

Hi,

is it possible to use a different k-mer size with HAWK? I see this definition in kmer.h:

#define KMER_LENGTH 31

and was wondering if by changing the value and compiling HAWK again will suffice or if there are any other internals that depend on a k-mer of 31.

Thanks,
Pedro

Edit:

other lines that seem to me to need to be changed are int kmerLength=31; in files kmersearch.cpp (L23), kmersummary.cpp (L25), bonf_fasta.cpp (L25) and hawk.cpp (L27). I can change all these lines for a different k-mer size, but my initial question remains: "changing the value and compiling HAWK again will suffice or if there are any other internals that depend on a k-mer of 31."

Edit 2:

and apparently there are also a few cases in hawk.cpp where one would need to change the calls to getKmer() in order to pass the correct value.