powered by
This function filters hypermutated samples and returns the formatted mutations with the appropriate trinucleotide context
format_muts(mutations, this_genome, filter_hyper_MB = NA)
a data frame called mutations which has been formatted with an extra column for trinucleotide context
A data frame with the following columns: chr, pos1, pos2, ref, alt, patient
autosomal chromosomes as chr1 to chr22 and sex chromosomes as chrX and chrY
the start position of the mutation in base 1 coordinates
the end position of the mutation in base 1 coordinates
the reference allele as a string containing the bases A, T, C or G
the alternate allele as a string containing the bases A, T, C or G
the patient identifier as a string
The reference genome object of BSgenome
The number of mutations per megabase for which a sample is considered hypermutated. Hypermutated samples will be removed in further analyses.
# \donttest{ data(cll_mutations) this_genome = BSgenome.Hsapiens.UCSC.hg19::Hsapiens formatted_mutations = format_muts(cll_mutations[1:10,], filter_hyper_MB = 30, this_genome = this_genome) # }
Run the code above in your browser using DataLab