Entering edit mode
Michael Cho
▴
10
@michael-cho-4714
Last seen 10.1 years ago
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]]