saveAlignments {enhancerHomologSearch} | R Documentation |
Save enhancer homologs to file in phylip format.
saveAlignments(al, output_folder = tempdir(), motifConsensus = NULL)
al |
output of alignment. |
output_folder |
output folder. |
motifConsensus |
Transcription factor binding consensus. |
The I/O status.
al <- readRDS(system.file("extdata", "al.rds", package="enhancerHomologSearch")) tmpfolder <- tempdir() library(MotifDb) motifs <- query(MotifDb, "JASPAR_CORE") consensus <- sapply(motifs, consensusString) consensus <- DNAStringSet(gsub("\\?", "N", consensus)) saveAlignments(al, output_folder=tmpfolder, motifConsensus=consensus)