I'm a complete novice in Biostrings and am starting to learn about
DNAStringSet. I was wondering how do you do a subseq operation on just one
of the DNAStrings in the set.
library('Biostrings')
a <- DNAStringSet('TTGGAACT')
subseq(a, start=2, end=3) <- NULL # works fine
x <- c(DNAStringSet('TTGGAACT'), DNAStringSet('AAGGTT'))
names(x) <- c('a', 'b')
# now I'm interested in doing this operation on just the a string and not b
subseq(x[['a']], start=2, end=3) <- NULL # doesn't work
a1 <- x[['a']]
subseq(a1, start=2, end=3) <- NULL # does work
Welcome any thoughts or suggestions!
Thanks,
Andrew
[[alternative HTML version deleted]]
_______________________________________________
Bioc-sig-sequencing mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
