Generate custom k-mer abundance from a given reference genome
GetCustomKmerCounts(k, ref.genome, custom.ranges, filter.path, verbose = FALSE)
Matrix of the counts of custom k-mer across the ref.genome
Length of k-mers (k>=2)
A ref.genome
argument as described in
ICAMS
.
A keyed data table which has custom ranges information. It
has three columns: chrom, start and end. It should use one-based coordinate
system. You can use the internal function in this package
ICAMS:::ReadBedRanges
to read a BED file in 0-based coordinates and
convert it to 1-based coordinates.
If given, homopolymers will be masked from genome(sequence). Only simple repeat masking is accepted now.
If TRUE
generate progress messages.