DottedAlignment {QSutils} | R Documentation |
Given an alignment, it takes the first sequence as reference, and depicts all equivalences in the alignment as dots, leaving only the variants with respect to the reference.
DottedAlignment(hseqs)
hseqs |
DNAStringSet or AAStringSet with haplotype sequences. |
A character string vector of the alignment, with dots in the conserved positions.
Mercedes Guerrero-Murillo and Josep Gregori
filepath<-system.file("extdata","ToyData_10_50_1000.fna", package="QSutils") lst <- ReadAmplSeqs(filepath,type="DNA") strs <- DottedAlignment(lst$hseqs) # Create a data frame to visualize the result. data.frame(Hpl=strs,stringsAsFactors=FALSE)