Displaying 1 result from an estimated 1 matches for "wholegenomematch".
2011 Apr 15
1
Whole genome searching of 100bp "D" sequence
...uot;chr14", "chr15", "chr16",
"chr17", "chr18", "chr19", "chr20", "chr21", "chr22", "chrX", "chrY")
#access/group the by subset of annotation
annotGr <- exons(txdb)
#annotGr <- cds(txdb)
wholeGenomeMatch <- function()
{
for(i in chrList)
{
#matches on plus strand
matches.plus.strand <- matchPattern(query.plus, Hsapiens[[i]],
fixed=FALSE, max.mismatch=5)
mp <- as.matrix(matches.plus.strand)
#matches on minus strand
matches.minus.strand <- matchPattern(query.m...