inst/vignette-source/test-03-align.R

In SVAlignR: Recovering Structure of Long Molecules from Structural Variation Data

library(SVAlignR)
data(longreads)
seqs <- longreads$connection[1:15]
pad <- c(rep("0", 9), rep("", 6))
names(seqs) <- paste("LR", pad, 1:length(seqs), sep = "")
seqs <- seqs[!duplicated(seqs)]

try( makeSubsMatrix(match = "13") )
try( makeSubsMatrix(mismatch = "13") )

mysub <- makeSubsMatrix(match = 2, mismatch = -6)

## test the 'alignCluster' method
try( alignCluster(13) )

ab <- alignCluster(seqs, mysub)
image(ab)