[Bioc-devel] export() to 2bit file

[Hervé Pagès]

Hi Michael,

   library(rtracklayer)
   library(Biostrings)
   x <- DNAStringSet("AAA-CCC-GGG-TTT-NNN-KKK")

Then:

   > x
     A DNAStringSet instance of length 1
       width seq
   [1]    23 AAA-CCC-GGG-TTT-NNN-KKK

   > export(x, "x.2bit")

   > import("x.2bit")
     A DNAStringSet instance of length 1
       width seq                                               names 

   [1]    23 AAATCCCTGGGTTTTTNNNTTTT                           1

What about having the "export" method for TwoBitFile raise an error
(or at least issue a warning) instead of silently turning everything
that is not A, C, G, T, or N into a T?

Thanks,
H.

-- 
Hervé Pagès

Program in Computational Biology
Division of Public Health Sciences
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N, M1-B514
P.O. Box 19024
Seattle, WA 98109-1024

E-mail: hpages at fhcrc.org
Phone:  (206) 667-5791
Fax:    (206) 667-1319