                                                                    'molecular data matrices for phylogenetic inference '
                                                                    'based on GenBank records.')

    # Add version argument
    ver = pkg_resources.require("MatPhylobi")[0].version
    parser.add_argument('-v', '--version', action='version',
                        version='%(prog)s {version}'.format(version=ver))  # TODO add also citation to the version

    # Prepare subparsers for different actions
    actions = parser.add_subparsers(title="program action", help="either run an analysis or update an older run",
                                    dest='action')
    actions.required = True

    parser_analyze = actions.add_parser('analyze')