- seq = string.join(seqArray, "")
- if seqLen < 500000:
- print "Added contig %s to database" % chromID
- xtGenome.addSequence(("xtropicalis", chromID), seq, "chromosome", str(seqLen))
- xtGenome.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
- xtGenome.addChromosomeEntry(chromID, outFileName, "file")