Entering edit mode
Is there a Bioconductor package that is able to convert a reference sequence + HGVS SVN specification into an edited sequence? Couldn't find one, but I might be just looking at the wrong packages.
Something like:
edited_sequence <- get_variant_seq("ENST00000357654.9:c.1045G>T")
and it returns a fasta string with the ENST00000357654.9 transcript sequence edited at position 1045.
Or something like this:
seq1 <- "ACGT"
mut1 <- "c.2C>G"
edited_sequence <- get_variant_seq(sequence = seq1, mutation = mut1)
and it return AGGT
.
I know GATK, bcf-tools, etc can do this, but I'd like to use R without external command line tools.