Hello Bioconductor curators,

How I can merge two subseqs of one DNAStringSet into one DNAStringSet?

seq.list=list()
seq.list[1]=paste(rep("AT",50),collapse="")
seq.list[2]=paste(rep("GC",50),collapse="")
example = DNAStringSet(unlist(seq.list))

length(example) = 2
width(example) = 100
start1 = 1
end1 = 20
start2 = 40
end2 = 60
first.subseq = subseq(example,start1,end1)
second.subseq = subseq(example,start2,end2)

c(first.subseq,second.subseq) doesn´t paste the sequences together

Is there an easy method? Or do i have to loop over all entries?

Thank you in anticipation

Kristian Ullrich
Institute of Plant Biochemistry



	[[alternative HTML version deleted]]

