- seq = string.join(seqArray, "")
- if seqLen < 250000:
- print "Added contig %s to database" % chromID
- drGenome.addSequence(("drerio", chromID), seq, "chromosome", str(seqLen))
- drGenome.addChromosomeEntry(chromID, chromID, "db")
- else:
- outFileName = "%s%s.bin" % (chromOutPath, chromID)
- outFile = open("%s%s" % (cisRoot, outFileName), "w")
- outFile.write(seq)
- outFile.close()
- print "Added contig file %s to database" % outFileName
- drGenome.addChromosomeEntry(chromID, outFileName, "file")
+ #print "Adding gene entries"
+ #loadGeneEntries(db, genePath)