Hi,

I'm trying to grab sequence from a GRangesList.  Each GRanges element has
several sequences which I am concatenating together and using to write a
fasta file (with the name of the GRanges as the name of the fasta sequence).
I can't seem to figure out a better / faster way than to loop like this:

seq<-getSeq(Hsapiens,myGRangesList[[i]],as.character=FALSE)
writeFASTA(paste(seq,collapse=""), file="myOutput.fasta",
desc=names(myGRangesList[i]), append=TRUE)

Any thoughts?
Thanks,
Michael

	[[alternative HTML version deleted]]

