Generate stranded k-mer abundance from a given genome and gene annotation file
GetStrandedKmerCounts(
k,
ref.genome,
stranded.ranges,
filter.path,
verbose = FALSE
)
Matrix of the counts of each stranded 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 stranded ranges information. It has four columns: chrom, start, end and strand. It should use one-based coordinate system.
If given, homopolymers will be masked from genome(sequence). Only simple repeat masking is accepted now.
If TRUE
generate progress messages.