[Bioc-sig-seq] trimLRPatterns question

Kunbin Qu KQu at genomichealth.com
Sat Oct 30 00:17:00 CEST 2010


Hi, all:

Is there a way to NOT trim a segment in the pattern from trimLRPatterns? For example:

t1<-DNAString("AAAATGGTAC")
t2<-DNAString("TGGTAC")
t3<-trimLRPatterns(Rpattern=t2, subject=t1, max.Rmismatch = rep(0,4), ranges=T, with.Rindels=TRUE)
t4<-DNAStringSet(t1, start=start(t3), end=end(t3))

> t4
  A DNAStringSet instance of length 1
    width seq
[1]     4 AAAA
> t1
  10-letter "DNAString" instance
seq: AAAATGGTAC
> t2
  6-letter "DNAString" instance
seq: TGGTAC
>


If I want to keep "TG" in t1, followed by the AAAA, which also appears in the beginning of the pattern t2, is there a way? The rep(0,4) does not help since it could find a perfect match at length 6 (the first try when starting to trim). I am trying to use "TG" as my anchor points to avoid false over-trimming. Thanks.

-Kunbin



______________________________________________________________________
The contents of this electronic message, including any attachments, are intended only for the use of the individual or entity to which they are addressed and may contain confidential information. If you are not the intended recipient, you are hereby notified that any use, dissemination, distribution, or copying of this message or any attachment is strictly prohibited. If you have received this transmission in error, please send an e-mail to postmaster at genomichealth.com and delete this message, along with any attachments, from your computer.



More information about the Bioc-sig-sequencing mailing list