get_lsingle_vcf: Function to parse Phred-scaled likelihoods

Description Usage Arguments Value

Description

This function parses the data column (GT:AD:DP:GQ:PL) to get the PL values -10log10(x), and converts them to a vector of likelihoods: L0,L1,L2. It's called within get_likelihood_vcf. newer version of workflow of GATK suggested calling the BAM for each individual first then combined them, so the format column can have different length for each variant like the example Lizhen provided.

Usage

1
get_lsingle_vcf(Mij, n_pl)

Arguments

Mij

a single data unit of format'GT:AD:DP:GQ:PL'.

n_pl

an integer to indicate the PL field.

Value

a vector of three values: l0, l1, l2


Struglab/RVS documentation built on May 9, 2019, 3:11 p.m.