Is there a strsplit method for DNAStringSet objects? I have a DNAStringSet object generated from the readFastq function (ShortRead package). In the middle of my reads is a barcode which I would like to use to split the reads into two reads (effectively treating the split read as paired-end reads).
Hi,
I just added an strsplit method for XStringSet objects to the Biostrings packages. So you shouldn't need to convert to character vector anymore before calling strsplit(). This conversion can be expensive for an XStringSet object with tens of millions of sequences in it.
It's in Biostrings 2.43.8 (BioC devel only), which will propagate and become available via biocLite() in the next 24 hours or so.
Cheers,
H.
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
version 2.3.6
The DNAStringSet object sequences can be convert to a character vector to be used in strsplit. For instance:
library(ShortRead) sp <- SolexaPath(system.file('extdata', package='ShortRead')) rfq <- readFastq(analysisPath(sp), pattern="s_1_sequence.txt") strsplit(as.character(id(rfq)), split="_1_1_1_") strsplit(as.character(sread(rfq)), split="CGCG")Ah so simple, silly me. Thanks for the answer.