call_resistance: CMV Resistance Genotyping from command line

Description Usage Arguments Value

View source: R/call_resistance.R

Description

Command line version of the website application Takes as input a VCF, varscan tab or fasta file. The program assumes variant files are generated relative to Merlin strain. Fasta files if not Whole Genome, or not aligned / assembled relative to Merlin are Processed using MAFFT & snp-sites. In this case the output files are returned to your working directory.

Usage

1
2
3
4
5
6
call_resistance(
  infile = system.file("testdata", "A10.vcf", package = "cmvdrg"),
  all_mutations = FALSE,
  inc_anecdotal = FALSE,
  outdir = ""
)

Arguments

infile

the input fasta, vcf or varscan.tab file

all_mutations

when FALSE only recognised resistant variants present are returned.

inc_anecdotal

include anectodat database entries in returned results?

outdir

for fasta input files intermediate alignment fasta & vcf files are generated, this defines the directory they are saved to. "out.fasta" "out.vcf"

Value

A data.frame containing resistance information for variants identified


ucl-pathgenomics/cmvdrg documentation built on Dec. 8, 2020, 2:36 a.m.