# #——消息= FALSE --------------------------------------------------------- 库(BRGenomics)库(Biostrings ) ## ----------------------------------------------------------------------------- # 得到路径包括2文件sfile < -系统。文件(“extdata”、“dm6_chr4chrM。2位”,包= " BRGenomics ") ## ----------------------------------------------------------------------------- 进口(sfile) seq_chr4 seq_chr4 < - ## ----------------------------------------------------------------------------- 数据(“txs_dm6_chr4”)txs_pr < -促进剂(txs_dm6_chr4 0 100 ) ## ----------------------------------------------------------------------------- seq_txs_pr < getSeq (TwoBitFile (sfile),txs_pr) seq_txs_pr ## ----------------------------------------------------------------------------- RNAStringSet (seq_txs_pr) suppressWarnings(翻译(seq_txs_pr)) oligonucleotideFrequency (seq_txs_pr[1:5],宽度= 1)oligonucleotideFrequency (seq_txs_pr[1:5],宽度= 2 ) ## ----------------------------------------------------------------------------- tss_seq < - getSeq (TwoBitFile (sfile)启动子(4)txs_dm6_chr4 4) tsspwm < - PWM (tss_seq) tsspwm