[BioC] Question about translate funciton in Biostrings package

ligang luzifer.li at gmail.com
Fri Mar 18 06:23:15 CET 2011


Dear list,

I'm using "tanslate" function in "Biostrings" package to translate DNA sequence
in proteins.

It did well when the base letter is "A/G/C/T"

But while the DNA sequence contain nucleotide ambiguity codes such as "N"/"M",
"tanslate" function did not work, for example:

translate(DNAString("AACTGTCGMCCC"))
#Error in translate(DNAStringSet(x)) : not a base at pos 9

translate(DNAString("AACTGNTCG"))
#Error in translate(DNAStringSet(x)) : not a base at pos 6


sessionInfo()
R version 2.12.1 (2010-12-16)
Platform: i386-pc-mingw32/i386 (32-bit)

locale:
[1] LC_COLLATE=Chinese_People's Republic of China.936  
LC_CTYPE=Chinese_People's
Republic of China.936    LC_MONETARY=Chinese_People's Republic of China.936
[4] LC_NUMERIC=C                                       LC_TIME=Chinese_People's
Republic of China.936    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] Biostrings_2.18.2 IRanges_1.8.9    

loaded via a namespace (and not attached):
[1] Biobase_2.10.0 tools_2.12.1  

---
LiGang



More information about the Bioconductor mailing list