库（org.hs.eg.db）EID < - select（org.hs.eg.db，“brca1”，“EntrezID”，“符号”）[[“EntrezID”]库（TXDB.hsapiens.ucc。HG19.knowngene）TXDB < - TXDB.HSAPIENS.CUCSC.HG19.KCONSC.HG19.KCNOKNGENE TXID < - 选择（TXDB，EID，“TXNAME”，“GENEEID”）[[TXNAME“] CD < - CDSBY（TXDB，BY =”Tx“，使用.names = true）brca1cds < - cds [cds（cds）％在％txid]库（gviz）plottracks（list（genomeaxistrack（），AnnotationTrack（BRCA1CD））库（BSGenome.hsapiens.ucsc.hg19）基因组< - bsgenome.hsapiens.ucsc.hg19 Tx_seq < - AstractTranscriptSeqs（基因组，BRCA1CDS）TX_Seq