[BioC] XstringViews and Fasta
Droit Arnaud
Arnaud.Droit at ircm.qc.ca
Thu Jul 9 20:13:47 CEST 2009
Hello,
I am using XstringViews object to transform bed file to fasta file :
as(chr, "XStringViews")
TopSeq<<-Views(chr, start=start, end=end,names=as.character(i))
write.XStringViews(TopSeq, file=paste("BED","_TEST_QUEST",".",method,as.character(nTop),"_",i,".fa",sep=""), format="fasta", width=10000)
Example :
>TopSeq
Views on a 242951149-letter DNAString subject
subject: CGTATCCCACACACCACACCCACACACCACACCC...NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
views:
start end width
[1] 162681889 162682089 201 [TATGTGAAATTTACATTCCATGA...TCCTTGAAGTTAGGAGCAACTT]
I cannot get only fasta sequence into variable. I am wondering if it possible to get sequence into variable.
Thanks
Arnaud Droit.
More information about the Bioconductor
mailing list