def parseSequenceData(self):
self.sequenceData = {}
self.qualityData = {}
for record in FastqReader( self.ccsFile ):
zmw = get_zmw( record.name )
new_record = FastqRecord(zmw, record.sequence, record.quality)
self.sequenceData[zmw] = new_record
self.qualityData[zmw] = meanPQv(new_record)
def parseSequenceData(self):
self.sequenceData = {}
self.qualityData = {}
for record in FastqReader(self.ccsFile):
zmw = get_zmw(record.name)
new_record = FastqRecord(zmw, record.sequence, record.quality)
self.sequenceData[zmw] = new_record
self.qualityData[zmw] = meanPQv(new_record)
def trimClusterNames(self, clusters):
trimmed = []
for cluster in clusters:
cluster = [get_zmw(c) for c in cluster]
trimmed.append( frozenset(cluster) )
return trimmed
def trimClusterNames(self, clusters):
trimmed = []
for cluster in clusters:
cluster = [get_zmw(c) for c in cluster]
trimmed.append(frozenset(cluster))
return trimmed
