rareVariantVis: Interactive visualization of rare variants on the chromosome

In RareVariantVis: A suite for analysis of rare genomic variants in whole genome sequencing data

Description

Reads file containing table of rare variants from one chromosome and provides adequate visualization. Input file can be obtained from function chromosomeVis. Function outputs visualization html figure in current working directory. Figure illustrates variants (dots) in their genomic coordinates (x axis). Ratio of alternative reads and depth (y axis) gives information about type of variant: homozygous alternative (expected ratio 1) and heterozygous (expected ratio 0.5). Zoom to the figures is possible, by marking the region of interest with mouse left click. Right click induces zoom out and return to the original plot. Pointing on variants provides basic information about the variant - gene name and position on chromosome.

Usage

rareVariantVis(input, outputFile, sample, chromosomes, append)

Arguments

input	Name of the input file (string) containing the table with rare variants generated by chromosomeVis. It can also be the variable with the table itself.
outputFile	Name of the output file (string) with visualisation.
sample	Name of the sample (used only for entitling the charts).
chromosomes	Vector of chromosome names (strings) to be included in the visualisation (all chromosomes by default). Chromosomes included in the parameter but absent in the table are ommited.
append	Logical value indicating whether charts should be appended to the output file without destroying it (FALSE by default).

Value

comp1

function returns html file with visualization of rare variants

Author(s)

Adam Gudys and Tomasz Stokowy

Examples

file = system.file("extdata", "RareVariants_CoriellIndex_S1.txt",
    package = "RareVariantVis")
rareVariantVis(input=file, "RareVariants_CoriellIndex_S1.html", "CorielIndex")

RareVariantVis documentation built on Nov. 8, 2020, 11:08 p.m.