blast - Search related sequences in databases¶

Blast(query, database, gap_open=11, gap_ext=1, matrix='BLOSUM62', blast_location=None, outfmt=0, filter_low_complexity=True)¶

Runs a protein vs. protein blast search. The results are returned according to the value of the outfmt parameter.

Parameters:

• query (seq.ConstSequenceHandle) – the query sequence

• database – The filename of the sequence database. Make sure that formatdb has been run on the database and the <database>.pin file exists.

• matrix – The substitution matrix to be used. Must be one of ‘BLOSUM45’, ‘BLOSUM62’, ‘BLOSUM80’, ‘PAM30’, ‘PAM70’.

• gap_open – Gap opening penalty. Note that only a subset of gap opening penalties is supported for each substitutition matrix. Consult the blast docs for more information.

• gap_ext – Gap extension penalty. Only a subset of gap extension penalties are supported for each of the substitution matrices. Consult the blast docs for more information.

• outfmt – output format, where ‘0’ corresponds to default output (parsed blast output and 1 to raw string output).

• filter_low_complexity – Mask off segments of the query sequence that have low compositional complexity, as determined by the SEG program of Wootton & Federhen (Computers and Chemistry, 1993)

Return type:

BlastHit (with outfmt=0) or str (with outfmt=1)

ParseBlastOutput(string, seqid_thres=0, evalue_thres=inf)¶

Parses the blast output and returns a list of BlastHits setting no seqid_thres or evalue_thres, restores default behaviour without filtering

class AlignedPatch(aln, bit_score, score, evalue, seqid)¶

An aligned patch, aka. HSP

aln¶

The local alignment. Sequence offset of both sequences in the alignment are set to the starting position in the query and target sequence, respectively.

Type:

AlignmentHandle

bit_score¶

The bit score of the HSP

score¶

The score of the HSP

evalue¶

The E-value of the HSP

seqid¶

The sequence identity of the HSP

class BlastHit(identifier, aligned_patches)¶

A positive match found by BLAST.

Each blast hit consists of one or more HSPs, encoded by the AlignedPatch class.

identifier¶

The identifier of the matching sequence

aligned_patches¶

list of AlignedPatch instances holding the actual HSPs.

class BlastError(brief, details)¶