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