- mmGenome.createIndices()
-
-
-def buildMmusculusDB(db=geneDB, downloadDir="%s/download" % cisRoot):
- genePath = "%s/seq_gene.md" % downloadDir # ftp://ftp.ncbi.nih.gov/genomes/M_musculus/mapview/seq_gene.md
- goDefPath = "%s/GO.terms_and_ids" % downloadDir # ftp://ftp.geneontology.org/pub/go/doc/GO.terms_and_ids
- goPath = "%s/gene2go" % downloadDir # ftp://ftp.ncbi.nih.gov/gene/DATA/gene2go.gz
- # chromosomes are from ftp://hgdownload.cse.ucsc.edu/goldenPath/mm9/chromosomes
- # but ignoring all random chromosomes
- chromDict = {"1": "%s/chr1.fa" % downloadDir,
- "2": "%s/chr2.fa" % downloadDir,
- "3": "%s/chr3.fa" % downloadDir,
- "4": "%s/chr4.fa" % downloadDir,
- "5": "%s/chr5.fa" % downloadDir,
- "6": "%s/chr6.fa" % downloadDir,
- "7": "%s/chr7.fa" % downloadDir,
- "8": "%s/chr8.fa" % downloadDir,
- "9": "%s/chr9.fa" % downloadDir,
- "10": "%s/chr10.fa" % downloadDir,
- "11": "%s/chr11.fa" % downloadDir,
- "12": "%s/chr12.fa" % downloadDir,
- "13": "%s/chr13.fa" % downloadDir,
- "14": "%s/chr14.fa" % downloadDir,
- "15": "%s/chr15.fa" % downloadDir,
- "16": "%s/chr16.fa" % downloadDir,
- "17": "%s/chr17.fa" % downloadDir,
- "18": "%s/chr18.fa" % downloadDir,
- "19": "%s/chr19.fa" % downloadDir,
- "X": "%s/chrX.fa" % downloadDir,
- "Y": "%s/chrY.fa" % downloadDir,
- "M": "%s/chrM.fa" % downloadDir
- }
-
- print "Creating database %s" % db
- createDBFile(db)
-
- print "Adding gene entries"
- loadGeneEntries(db, genePath, chromDict)
-
- print "Adding gene features"
- loadGeneFeatures(db, genePath, chromDict)
-
- print "Adding gene annotations"
- loadGeneAnnotations(db)