Patrick,
This adapter tool looks extremely useful for my purposes: removing
adapters from smRNA reads to estimate the short template lengths.
Forgive me if the answer to this is obvious, but everything seems
to
work with trimLRPatterns, except that it doesn't seem to allow the
Rpattern or Lpattern to slide along the sequence (at least using
the
default settings--see below). Rather it looks only for exact
matches,
that leave no overhang. Thus:
Rpattern <- "CTGTAGGCACCA"
trims:
[6] 34 GCTGGAACCCAGGGTGTTGTACCTGTAGGCACCA
nicely, to:
[6] 22 GCTGGAACCCAGGGTGTTGTAC
but a sequence where resulting in an Rpattern overhang (here ~2nt):
[90] 34 GTAAGACCATACTTGGCCGAATGCCTGTAGGCAC
is not trimmed at all:
[90] 34 GTAAGACCATACTTGGCCGAATGCCTGTAGGCAC
:
What can I do to allow for flexibility at the overhanging end?
Again, thanks very much.
Joe
On 14 Jan 2009, at 19:17, Patrick Aboyoun wrote:
I just checked in a trimLRPatterns function to the Bioconductor svn
repository for BioC 2.4. Its signature is
trimLRPatterns(Lpattern = NULL, Rpattern = NULL, subject,
max.Lmismatch = 0, max.Rmismatch = 0,
with.Lindels = FALSE, with.Rindels = FALSE,
Lfixed = TRUE, Rfixed = TRUE, ranges = FALSE)
As you can infer from the arguments, this function allows the
user to
set the # of mismatches (if with.*indels = FALSE) / edit distance
(if
with.*indels = TRUE) for the left and right flanking "patterns". It
also allows for IUPAC ambiguity letters in these flanking regions
if
*fixed = FALSE. When ranges = FALSE, trimLRPatterns returns the
trimmed
strings. When ranges = TRUE, it returns the ranges that you can
use to
trim the strings. Here are some examples:
Lpattern <- "TTCTGCTTG"
Rpattern <- "GATCGGAAG"
subject <- DNAString("TTCTGCTTGACGTGATCGGA")
subjectSet <- DNAStringSet(c("TGCTTGACGGCAGATCGG",
"TTCTGCTTGGATCGGAAG"))
trimLRPatterns(Lpattern = Lpattern, subject = subject)
11-letter "DNAString" instance
seq: ACGTGATCGGA
trimLRPatterns(Lpattern = Lpattern, Rpattern = Rpattern, subject =
subjectSet)
A DNAStringSet instance of length 2
width seq
[1] 18 TGCTTGACGGCAGATCGG
[2] 0
trimLRPatterns(Lpattern = Lpattern, Rpattern = Rpattern, subject =
subjectSet,
+ ranges = TRUE)
IRanges object:
start end width
1 1 18 18
2 10 9 0
This functionality will be available on bioconductor.org (and
downloadable via biocLite) in the next day or so, but you can
also grab
Biostrings from svn directly if you need it sooner. It will also
feed
its way into Biostrings documentation and training material
before the
next release of Bioconductor in May.
Patrick
Patrick Aboyoun wrote:
David,
Following up on Martin's comments, I am putting the finishing
touches on a function called trimLRPatterns for the Biostrings
package. Its purpose is to trim left and/or right flanking
patterns from sequences, so it can strip 5' and/or 3' adapters
from your reads. The signature for this function is
trimLRPatterns(Lpattern=NULL, Rpattern=NULL, subject,
max.Lnedit=0, max.Rnedit=0,
with.Lindels=FALSE, with.Rindels=FALSE,
Lfixed=TRUE, Rfixed=TRUE,
rangesOnly = FALSE)
I will be checking this function into the BioC 2.4 code line,
which requires using R-devel, sometime today or tomorrow. I
will send out an e-mail to this group when I check it in and
show a simple example of its usage. I talked with Martin and
he will wrap this functionality in the ShortRead layer so you
don't have to leave the ShortRead class system when removing
adapters from your reads.
Cheers,
Patrick
