Last chance! 50% off unlimited learning
Sale ends in
pairwiseAlignment
. This function uses default of type="overlap", gapOpening=-1, and gapExtension=-1 to align the patternSeq against subjectSeqs. One can adjust these parameters if prefered, but not recommended. This function is meant for aligning a short pattern onto large collection of subjects. If you are looking to align a vector sequence to subjects, then please use BLAT or see one of following blatSeqs
, findAndRemoveVector
pairwiseAlignSeqs(subjectSeqs = NULL, patternSeq = NULL, side = "left",
qualityThreshold = 1, showStats = FALSE, bufferBases = 5, doRC = TRUE,
returnUnmatched = FALSE, returnLowScored = FALSE, parallel = FALSE, ...)
BiocParallel
. Defaults to FALSE. If no parallel backend is
registered, then a serial version is ran using SerialParam
.pairwiseAlignment
primerIDAlignSeqs
, vpairwiseAlignSeqs
,
doRCtest
, findAndTrimSeq
, blatSeqs
,
findAndRemoveVector
subjectSeqs <- c("CCTGAATCCTGGCAATGTCATCATC", "ATCCTGGCAATGTCATCATCAATGG",
"ATCAGTTGTCAACGGCTAATACGCG", "ATCAATGGCGATTGCCGCGTCTGCA",
"CCGCGTCTGCAATGTGAGGGCCTAA", "GAAGGATGCCAGTTGAAGTTCACAC")
subjectSeqs <- DNAStringSet(xscat("AAAAAAAAAA", subjectSeqs))
pairwiseAlignSeqs(subjectSeqs, "AAAAAAAAAA", showStats=TRUE)
pairwiseAlignSeqs(subjectSeqs, "AAATAATAAA", showStats=TRUE,
qualityThreshold=0.5)
Run the code above in your browser using DataLab