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]]

