maf-cull
========

This script aims to "cull" redundant alignments.  It removes
alignments whose coordinates in the top-most sequence are contained in
LIMIT or more higher-scoring alignments.

Alignments using opposite strands of the top-most sequence are not
considered to contain each other.

The alignments must have been sorted beforehand by maf-sort.  If not,
maf-cull will quit with an error message.

You can use it like this::

  maf-cull.py --limit=3 sorted-alignments.maf > culled.maf
  maf-cull.py -l3 sorted-alignments.maf > culled.maf

If you want to cull based on a sequence other than the top one, use
maf-swap to bring that sequence to the top.

This type of culling can also be done by NCBI BLAST, and it is
described in this article:

  Winnowing sequences from a database search.
  Berman P, Zhang Z, Wolf YI, Koonin EV, Miller W.
  J Comput Biol. 2000 7:293-302.

Options
-------

  -h, --help                  Print a help message and exit.
  -l LIMIT, --limit=LIMIT     Set the culling limit to LIMIT.
