+def buildSpurpuratusDB(db=geneDB):
+ chromoPath = "%s/download/Spur2.1_Nmasked.txt" % cisRoot
+ chromoOutPath = "/S_purpuratus/"
+
+ print "Creating database %s" % db
+ createDBFile(db)
+
+ print "Loading genomic sequence"
+ loadChromosome(db, chromoPath, chromoOutPath)
+
+ print "Creating Indices"
+ createDBindices(db)
+
+ print "Finished creating database %s" % db
+
+def createDBFile(db):
+ spGenome = Genome("spurpuratus", version="2.1", dbFile=db)
+ spGenome.createGeneDB(db)
+
+