Generate exome k-mer abundance from a given reference genome
GetExomeKmerCounts(k, ref.genome, exome.ranges, filter.path,
verbose = FALSE)
Length of k-mers (k>=2)
A ref.genome
argument as described in
ICAMS
.
If given, homopolymers will be masked from genome(sequence). Only simple repeat masking is accepted now.
If TRUE
generate progress messages.
A keyed data table which has exome ranges information. It has three columns: chrom, start and end.
Matrix of the counts of exome k-mer across the ref.genome