Reading the source and example scripts is the best way to get familiar with bx-python. It has a lot of useful functionality so it's well worth the time. Here's a script that:

• Indexes the MAF file
• Retrieves all alignments in a region (organism, chromosome, start, end)
• Slices the alignments by the region of interest
• Builds a dictionary mapping organism+chromosome to sequence

Hope this helps with getting started:

import os
import sys
import subprocess

from bx.align import maf

def main(in_file, org, chrom, start, end):
    index_file = in_file + ".index"
    if not os.path.exists(index_file):
        build_index(in_file)

    region_name = "%s.%s" % (org, chrom)
    start = int(start)
    end = int(end)
    index = maf.Indexed(in_file, index_file)
    for align in index.get(region_name, start, end):
        region_align = align.slice_by_component(region_name, start, end)
        seqs_by_org = dict()
        for component in region_align.components:
            seqs_by_org[component.src] = component.text
        print seqs_by_org

def build_index(in_file):
    cl = ["maf_build_index.py", in_file]
    subprocess.check_call(cl)

if __name__ == "__main__":
    main(*sys.argv[1:])

Usage is:

% python maf_retrieve_region.py chr22.maf hg18 chr22 14430574 14430584
{'otoGar1.scaffold_92626.1-103747': 'tctttagaag',
 'tupBel1.scaffold_3803.1-85889': 'TTTTTAGGAG',
 'hg18.chr22': 'tcttaaaaag',
 'panTro2.chrUn': 'ttttaaaaag'}