sequences.py 文件源码

def _get_merged_reads_sequence(self, merged_reads: MergedReads):
        if len(merged_reads.reads) != len(merged_reads.prefix_lengths)+1:
            raise ValueError(
                "Invalid `MergedReads` data structure, not enough information "
                "on read prefix lengths."
            )

        sequence_parts = []
        for read, prefix_len in zip_longest(merged_reads.reads,
                                            merged_reads.prefix_lengths):
            sequence = self.get_sequence(read)
            if prefix_len:
                sequence_parts.append(sequence[:prefix_len])
            else:
                sequence_parts.append(sequence)

        return b"".join(sequence_parts)