DiffBlocks

class geneblocks.DiffBlocks(s1, s2, blocks)

Class to generate and represent DiffBlocks.

Usage:

>>> DiffBlocks.from_sequences(s1, s2)

static from_sequences(s1, s2, blast_over=500, max_complexity=100000000.0)

Create DiffBlocks by comparing two sequences.

Parameters:

• s1 – Two sequences, either “ATGC…” strings or Biopython records.

• s2 – Two sequences, either “ATGC…” strings or Biopython records.

• blast_over – A blast will be triggered to accelerate homology finding if len(s1) + len(s2) > blast_over.

• max_complexity – If len(s1) * len(s2) is over max_complexity, no analysis is done and s1 is just labeled as a “change” of s2 (useful internally during the recursions of this method).