ICAMS (version 2.0.10)

GetExomeKmerCounts: Generate exome k-mer abundance from a given reference genome

Description

Generate exome k-mer abundance from a given reference genome

Usage

GetExomeKmerCounts(k, ref.genome, exome.ranges, filter.path,
  verbose = FALSE)

Arguments

k

Length of k-mers (k>=2)

ref.genome

A ref.genome argument as described in ICAMS.

filter.path

If given, homopolymers will be masked from genome(sequence). Only simple repeat masking is accepted now.

verbose

If TRUE generate progress messages.

exome.range

A keyed data table which has exome 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.

Value

Matrix of the counts of exome k-mer across the ref.genome