flip=False, verbose=False, stripSpace=False, cachePages=100000,
propertyList=[]):
- delimiter = "|"
+ writeLog("%s.log" % outdbname, verstring, string.join(sys.argv[1:]))
+ geneDict = {}
dataType = "DNA"
if genedatafilename is not None:
dataType = "RNA"
genedatafile = open(genedatafilename)
-
-
- forcePair = False
- if forceID is not None:
- forcePair = True
- else:
- forceID = 0
-
- maxBorder = 0
- index = 0
- insertSize = 100000
-
- writeLog("%s.log" % outdbname, verstring, string.join(sys.argv[1:]))
-
- geneDict = {}
- mapDict = {}
- if dataType == "RNA":
for line in genedatafile:
fields = line.strip().split("\t")
blockCount = int(fields[7])
totalLength += exonLengths[index]
geneDict[uname] = (sense, blockCount, totalLength, chrom, chromstarts, exonLengths)
- mapDict[uname] = []
genedatafile.close()
fields = line.split()
readsize = len(fields[5])
pairedTest = fields[0][-2:]
+ forcePair = False
+ if forceID is not None:
+ forcePair = True
+ else:
+ forceID = 0
+
paired = False
if pairedTest in ["/1", "/2"] or forcePair:
print "assuming reads are paired"
paired = True
-
print "read size: %d bp" % readsize
if init:
rds.insertMetadata([("readsize", readsize)])
infile.close()
- trim = -4
+ maxBorder = 0
if dataType == "RNA":
+ trim = -4
maxBorder = readsize + trim
infile = open(filename, "r")
mInsertList = []
sInsertList = []
index = uIndex = mIndex = sIndex = lIndex = 0
+ delimiter = "|"
+ insertSize = 100000
for line in infile:
lIndex += 1
if stripSpace: