importVCF: Import a VCF directly into MTseeker

Description Usage Arguments Value Examples

View source: R/importVCF.R

Description

Import a VCF directly into MTseeker

Usage

1
importVCF(x, gzip = TRUE, biscuit = FALSE, verbose = TRUE)

Arguments

x

character, list or dataframe for file name of VCF(and BAM)

gzip

Is the input gzipped? (FIXME: autodetect this)

biscuit

Is the input VCF from biscuit? (FALSE)

verbose

Print all messages during progress? (TRUE)

Value

an MVRangesList

Examples

1
2
3
4
5
6
7
library(MTseekerData)
VCFdir <- system.file("extdata", package="MTseekerData")
VCF <- file.path(VCFdir, list.files(VCFdir, pattern="*.vcf.gz$"))
mvr <- importVCF(VCF)[[1]]
mvr$FILTER <- ifelse(totalDepth(mvr) >= 20, "PASS", ".")
mvr$PASS <- mvr$FILTER == "PASS"
filterMT(mvr)

trichelab/MTseeker documentation built on March 8, 2021, 6:20 p.m.