[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