Metadata-Version: 2.0
Name: moca
Version: 0.4.3
Summary: Tool for motif conservation analysis
Home-page: https://github.com/saketkc/moca
Author: Saket Choudhary
Author-email: saketkc@gmail.com
License: ISC
Keywords: moca
Platform: UNKNOWN
Classifier: Development Status :: 4 - Beta
Classifier: Intended Audience :: Developers
Classifier: License :: OSI Approved :: ISC License (ISCL)
Classifier: Natural Language :: English
Classifier: Programming Language :: Python :: 2
Classifier: Programming Language :: Python :: 2.6
Classifier: Programming Language :: Python :: 2.7
Classifier: Programming Language :: Python :: 3.3
Classifier: Programming Language :: Python :: 3.5
Classifier: Topic :: Scientific/Engineering :: Bio-Informatics
Classifier: Topic :: Scientific/Engineering :: Visualization
Classifier: Operating System :: POSIX
Classifier: Operating System :: Unix
Classifier: Operating System :: MacOS
Requires-Dist: biopython (>=1.68)
Requires-Dist: click (>=6.6)
Requires-Dist: click-help-colors (>=0.3)
Requires-Dist: coverage (>=4.2)
Requires-Dist: future (>=0.16.0)
Requires-Dist: matplotlib (>=2.0.0)
Requires-Dist: numpy (>=1.11.3)
Requires-Dist: pandas (>=0.19.2)
Requires-Dist: pyBigWig (>=0.2.8)
Requires-Dist: pybedtools (>=0.7.9)
Requires-Dist: pytest (>=3.0.5)
Requires-Dist: pytest-cov (>=2.4.0)
Requires-Dist: pytest-mpl (>=0.5)
Requires-Dist: scipy (>=0.18.1)
Requires-Dist: seaborn (>=0.7.1)
Requires-Dist: six (>=1.10.0)
Requires-Dist: statsmodels (>=0.6.1)
Requires-Dist: tqdm (>=4.7.2)

==========================================
MoCA: Tool for MOtif Conservation Analysis
==========================================

.. image:: https://img.shields.io/pypi/v/moca.svg
        :target: https://pypi.python.org/pypi/moca/

.. image:: https://img.shields.io/travis/saketkc/moca.svg
        :target: https://travis-ci.org/saketkc/moca

.. image:: https://img.shields.io/badge/install%20with-bioconda-brightgreen.svg
        :target: https://bioconda.github.io/recipes/moca/README.html

.. image:: https://coveralls.io/repos/github/saketkc/moca/badge.svg?branch=master
        :target: https://coveralls.io/github/saketkc/moca?branch=master

.. image:: https://landscape.io/github/saketkc/moca/master/landscape.svg?style=flat
        :target: https://landscape.io/github/saketkc/moca/master

.. image:: https://requires.io/github/saketkc/moca/requirements.svg?branch=master
        :target: https://requires.io/github/saketkc/moca/requirements/?branch=master


LICENSE
-------
ISC



Installation
------------


Requirements
~~~~~~~~~~~~

* pybedtools
* biopython
* pandas
* scipy
* statsmodels
* pybigwig
* seaborn
* MEME==4.10.2

NOTE: MoCA also relies on `fasta-shuffle-letters` that was introduced in MEME `4.11.0`
hence if you are using `4.10.2` make sure the `fasta-shuffle-letters` is the updated one.

For a sample script see `travis/install_meme.sh`

Using Conda
~~~~~~~~~~~
``moca`` is most compatible with the `conda`_ environment.

::

    $ conda config --add channels bioconda
    $ conda install moca


Using pip
~~~~~~~~~

::

   $ pip install moca


For development
~~~~~~~~~~~~~~~

::

    $ git clone https://github.com:saketkc/moca.git
    $ cd moca
    $ conda env create -f environment.yml python=2.7
    $ source activate mocadev
    $ python setup.py install



Workflow
--------

MoCA makes use of PhyloP/PhastCons/GERP scores to assess the quality of a
motif, the hypothesis being a 'true motif' would evolve slower as compared
to its surrounding(flanking sequences).

.. image:: https://raw.githubusercontent.com/saketkc/moca_web/master/docs/abstract/workflow.png


Usage
-----

::

    $ moca
    Usage: moca [OPTIONS] COMMAND [ARGS]...

      moca: Motif Conservation Analysis

    Options:
      --version  Show the version and exit.
      --help     Show this message and exit.

    Commands:
      find_motifs  Run meme to locate motifs and create...
      plot         Create stacked conservation plots



Motif analysis using MEME
~~~~~~~~~~~~~~~~~~~~~~~~~

MoCA can perform motif analysis for you given a bedfile containing
ChIP-Seq peaks.

Genome builds and MEME binary locations are specified through a configuraton file.
A sample configuration file is available: `tests/data/application.cfg` and should be
self-explanatory.

moca find_motifs
~~~~~~~~~~~~~~~~


::

    $ moca find_motifs -h
    Usage: moca find_motifs [OPTIONS]

      Run meme to locate motifs and create conservation stacked plots

    Options:
      -i, --bedfile TEXT            Bed file input  [required]
      -o, --oc TEXT                 Output Directory  [required]
      -c, --configuration TEXT      Configuration file  [required]
      --slop-length INTEGER         Flanking sequence length  [required]
      --flank-motif INTEGER         Length of sequence flanking motif  [required]
      --n-motif INTEGER             Number of motifs
      -t, --cores INTEGER           Number of parallel MEME jobs  [required]
      -g, -gb, --genome-build TEXT  Key denoting genome build to use in
                                    configuration file  [required]
      --show-progress               Print progress
      -h, --help                    Show this message and exit.


moca plot
~~~~~~~~~


::

    $ moca plot -h
    Usage: moca plot [OPTIONS]

      Create stacked conservation plots

    Options:
      --meme-dir, --meme_dir TEXT     MEME output directory  [required]
      --centrimo-dir, --centrimo_dir TEXT
                                      Centrimo output directory  [required]
      --fimo-dir-sample, --fimo_dir_sample TEXT
                                      Sample fimo.txt  [required]
      --fimo-dir-control, --fimo_dir_control TEXT
                                      Control fimo.txt  [required]
      --name TEXT                     Plot title
      --flank-motif INTEGER           Length of sequence flanking motif
                                      [required]
      --motif INTEGER                 Motif number
      -o, --oc TEXT                   Output Directory  [required]
      -c, --configuration TEXT        Configuration file  [required]
      --show-progress                 Print progress
      -g, -gb, --genome-build TEXT    Key denoting genome build to use in
                                      configuration file  [required]
      -h, --help                      Show this message and exit.


Example
-------

Most users will require using the command line version only:

::

    $ moca find_motifs -i encode_test_data/ENCFF002DAR.bed\
        -c tests/data/application.cfg -g hg19 --show-progress



Creating plots if you already have run MEME and Centrimo:

::

    $ moca plot -c tests/data/application.cfg -g hg19\
        --meme-dir moca_output/meme_out\
        --centrimo-dir moca_output/centrimo_out\
        --fimo-dir-sample moca_output/meme_out/fimo_out_1\
        --fimo-dir-control moca_output/meme_out/fimo_random_1\
        --name ENCODEID


.. image:: http://www.saket-choudhary.me/moca/_static/img/ENCFF002CEL.png


There is also a structured API available,
however it might be missing examples and documentation at places.

API Documentation
-----------------

http://saketkc.github.io/moca/



Tests
-----
``moca`` is mostly extensively tested. See `code-coverage`_. 

Run tests locally

::

    $ ./runtests.sh


Credits
---------

This package was created with Cookiecutter_ and the `audreyr/cookiecutter-pypackage`_ project template.

.. _`MoCA0.1.0`: https://github.com/saketkc/moca_web
.. _Cookiecutter: https://github.com/audreyr/cookiecutter
.. _`audreyr/cookiecutter-pypackage`: https://github.com/audreyr/cookiecutter-pypackage
.. _`conda`: http://conda.pydata.org/docs/using/using.html
.. _`code-coverage`: https://coveralls.io/github/saketkc/moca?branch=master


=======
History
=======
0.4.1 (2017-02-16)
------------------
* Added $\delta$ options for sample vs control ttest
* Removed unused requirements: ipaddress
* Fixed meme runner to plot

0.3.3 (2016-10-03)
------------------
* Removed pycairo dependency

0.2.9 (2016-05-31)
------------------

* Do not fail silently on MEME failing
* Support --cores to support parallel threads

0.2.8 (2016-05-30)
------------------
* Fixed MEME pipeline missing from mocacli

0.2.7 (2016-05-30)
------------------
* Fixed bug where missing wig keys were not handled in mocacli

0.2.4 (2016-05-29)
------------------

* Cleaned up unused scripts under scripts directory
* Add configuration file example


0.2.3 (2016-05-29)
------------------
* Include package_dir in setup.py
* Include requirements.txt in MANIFEST

0.2.0 (2016-05-29)
------------------

* First release on PyPI.


