def main():

    actions = (
        ("breakpoint", "find scaffold breakpoints using genetic map"),
        ("ld", "calculate pairwise linkage disequilibrium"),
        ("bed", "convert MSTmap output to bed format"),
        ("fasta", "extract markers based on map"),
        ("anchor", "anchor scaffolds based on map"),
        ("rename", "rename markers according to the new mapping locations"),
        ("header", "rename lines in the map header"),
        # Plot genetic map
        ("blat", "make ALLMAPS input csv based on sequences"),
        ("dotplot", "make dotplot between chromosomes and linkage maps"),
    )
    p = ActionDispatcher(actions)
    p.dispatch(globals())