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
[email protected] and delete this message, along with any
attachments, from your computer.
_______________________________________________
Bioc-sig-sequencing mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
