- if count > maxMultiReadCount:
- totalReadCounts["multiDiscard"] += 1
- else:
- multiInsertList.append(getRDSEntry(readData, rdsEntryName, chrom, readsize, weight=count))
- totalReadCounts["multi"] += 1
+ multiInsertList.append(getRDSEntry(readData, rdsEntryName, chrom, readsize, weight=count))
+
+ if dataType == "RNA":
+ for entry in spliceReadDict.keys():
+ (readData, rdsEntryName) = spliceReadDict[entry]
+ chrom = samfile.getrname(readData.rname)
+ spliceInsertList.append(getRDSSpliceEntry(readData, rdsEntryName, chrom, readsize))
+ totalReadCounts["splice"] += 1
+
+ for entry in multispliceReadDict.keys():
+ (readData, count, rdsEntryName) = multispliceReadDict[entry]
+ chrom = samfile.getrname(readData.rname)
+ multispliceInsertList.append(getRDSSpliceEntry(readData, rdsEntryName, chrom, readsize, weight=count))
+ totalReadCounts["multisplice"] += 1