biopython v1.71.0 Bio.AlignIO.MafIO

Bio.AlignIO support for the “maf” multiple alignment format.

The Multiple Alignment Format, described by UCSC, stores a series of multiple alignments in a single file. It is suitable for whole-genome to whole-genome alignments, metadata such as source chromosome, start position, size, and strand can be stored.


You are expected to use this module via the Bio.AlignIO functions(or the Bio.SeqIO functions if you want to work directly with the gapped sequences).

Coordinates in the MAF format are defined in terms of zero-based start positions (like Python) and aligning region sizes.

A minimal aligned region of length one and starting at first position in the source sequence would have start == 0 and size == 1.

As we can see on this example, start + size will give one more than the zero-based end position. We can therefore manipulate start and start + size as python list slice boundaries.

For an inclusive end coordinate, we need to use end = start + size - 1. A 1-column wide alignment would have start == end.

Link to this section Summary


Iterate over a MAF file handle as MultipleSeqAlignment objects

Link to this section Functions

Iterate over a MAF file handle as MultipleSeqAlignment objects.

Iterates over lines in a MAF file-like object (handle), yielding MultipleSeqAlignment objects. SeqRecord IDs generally correspond to species names.