maf2spec | R Documentation |
The function converts a MAF file into a spectrum matrix which has mutation counts catalogued by contexts. It uses the fact that DNA has two complementary strands and folds the A/G mutations into to T/C.
It is a wrapper on bedtools(https://github.com/arq5x/bedtools2/releases).
maf2spec( bedtools_path = "bedtools", maf, ref_genome, context_length = 1, output_file )
bedtools_path |
Path of bedtool executable, e.g. "~/bedtools/bin/bedtool" |
maf |
A maf file |
ref_genome |
An fa file path of reference genome sequence |
context_length |
The length of extension of each side of the mutation. By default it is 3 (trinucleotide context). |
output_file |
An output file that contains the context matrix |
Returns a matrix of mutation counts in different context, in each of the samples
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.