Package: BSgenome.Cfamiliaris.UCSC.canFam2 Title: Full genome sequences for Canis lupus familiaris (UCSC version canFam2) Description: Full genome sequences for Canis lupus familiaris (Dog) as provided by UCSC (canFam2, May 2005) and stored in Biostrings objects. Version: 1.4.2 Suggests: GenomicFeatures organism: Canis lupus familiaris common_name: Dog provider: UCSC provider_version: canFam2 release_date: May 2005 release_name: Broad Institute v2.0 source_url: http://hgdownload.cse.ucsc.edu/goldenPath/canFam2/bigZips/ organism_biocview: Canis_familiaris BSgenomeObjname: Cfamiliaris seqnames: paste("chr", c(1:38, "X", "M", "Un"), sep="") circ_seqs: "chrM" SrcDataFiles: chromFa.tar.gz from http://hgdownload.cse.ucsc.edu/goldenPath/canFam2/bigZips/ PkgExamples: genome$chr1 # same as genome[["chr1"]] . ## --------------------------------------------------------------------- ## Upstream sequences ## --------------------------------------------------------------------- ## The upstream sequences for canFam2 (i.e. the sequences 1000 bases ## upstream of annotated transcription starts) can easily be extracted ## from the full genome sequences with something like: . library(GenomicFeatures) txdb <- makeTxDbFromUCSC("canFam2", "refGene") gn <- sort(genes(txdb)) up1000 <- flank(gn, width=1000) up1000seqs <- getSeq(genome, up1000) . ## IMPORTANT: Make sure you use a TxDb package (or TxDb object), ## that contains a gene model based on the exact same reference genome ## as the BSgenome object you pass to getSeq(). Note that you can make ## your own custom TxDb object from various annotation resources. ## See the makeTxDbFromUCSC(), makeTxDbFromBiomart(), ## and makeTxDbFromGFF() functions in the GenomicFeatures ## package. seqs_srcdir: /fh/fast/morgan_m/BioC/BSgenomeForge/srcdata/BSgenome.Cfamiliaris.UCSC.canFam2/seqs