Identify mitochondrial variants

Identify mitochondrial variants present in single cells.

IdentifyVariants(object, ...)

# S3 method for default
IdentifyVariants(
  object,
  refallele,
  stabilize_variance = TRUE,
  low_coverage_threshold = 10,
  verbose = TRUE,
  ...
)

# S3 method for Assay
IdentifyVariants(object, refallele, ...)

# S3 method for Seurat
IdentifyVariants(object, refallele, assay = NULL, ...)

Arguments

object	A Seurat object
...	Arguments passed to other methods
refallele	A dataframe containing reference alleles for the mitochondrial genome.
stabilize_variance	Stabilize variance
low_coverage_threshold	Low coverage threshold
verbose	Display messages
assay	Name of assay to use. If NULL, use the default assay.

Value

Returns a dataframe

Examples

if (FALSE) {
data.dir <- "path/to/data/directory"
mgatk <- ReadMGATK(dir = data.dir)
variant.df <- IdentifyVariants(
  object = mgatk$counts,
  refallele = mgatk$refallele
)
}