get_pos_from_end: Retrieve variant distance from end of read

In mmaitenat/ideafix: ideafix, DEAmination FIXing

Description

get_pos_from_end retrieves, for each C:G > T:A variant present in vcf_filename, the median value of the position of the variant (starting from the end of the read) in those reads that support the alternate allele. This data is retrieved from vcf_filename.

Usage

get_pos_from_end(vcf_filename, samplename, ...)

Arguments

vcf_filename	character string naming the path to the input vcf, i.e. the vcf file containing the variants to classify. This file must have been generated with Mutect2, either in tumor only or tumor/normal mode with strand bias annotation enabled.
samplename	character string naming the sample in vcf_filename. This must match the name given to the sample when running Mutect2. It can be obtained by calling get_samplename on vcf_filename.
...	additional arguments not to be checked.

Value

Numeric vector containing the median distance from end of read of each C:G > T:A variant in vcf_filename, ordered by appearance in the file.

mmaitenat/ideafix documentation built on Sept. 18, 2021, 7:55 a.m.