[BioC] XstringViews and Fasta
Arnaud.Droit at ircm.qc.ca
Thu Jul 9 21:12:37 CEST 2009
Thank you very much for your help.
On 09/07/09 3:00 PM, "Patrick Aboyoun" <paboyoun at fhcrc.org> wrote:
Could you provide more description, or code, for your situation? In particular what does "I cannot get only fasta sequence into variable. I am wondering if it possible to get sequence into variable." mean? In your code below, TopSeq is an XStringViews object. If you use the DNAStringSet constructor, you can create a DNAStringSet by DNAStringSet(TopSeq). Is that what you are looking to create in R?
Droit Arnaud wrote:
> 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 :
> Views on a 242951149-letter DNAString subject
> subject: CGTATCCCACACACCACACCCACACACCACACCC...NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
> start end width
>  162681889 162682089 201 [TATGTGAAATTTACATTCCATGA...TCCTTGAAGTTAGGAGCAACTT]
> I cannot get only fasta sequence into variable. I am wondering if it possible to get sequence into variable.
> Arnaud Droit.
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor