gene_midas_data: Get all MIDAS data from a gene

Description Usage Arguments Value

View source: R/midas_data.r

Description

Returns a list similar to the one by read_midas_data but only with sites from a given gene.

Usage

1
gene_midas_data(Dat, gene, depth_thres = 1, freq_thres = 0.5)

Arguments

Dat

A list. See output from read_midas_data

gene

Name of the gene to be extracted, must exist in 'gene_id' column of the 'info' element from Dat.

depth_thres

Minimum depth trheshold.

freq_thres

Freq threshold.

Value

A tibble that results of merging the info, dat and freq elements of Dat. Plus allele assignments.


surh/HMVAR documentation built on Aug. 18, 2021, 1:21 a.m.