Metadata-Version: 2.0
Name: moca
Version: 0.3.2
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: Pillow (==3.3.1)
Requires-Dist: PyYAML (==3.11)
Requires-Dist: biopython (==1.68)
Requires-Dist: certifi (==2016.8.31)
Requires-Dist: cffi (==1.8.3)
Requires-Dist: click (==6.6)
Requires-Dist: click-help-colors (==0.3)
Requires-Dist: coverage (==4.2)
Requires-Dist: cryptography (==1.5.2)
Requires-Dist: cycler (==0.10.0)
Requires-Dist: enum34 (==1.1.6)
Requires-Dist: functools32 (==3.2.3.post2)
Requires-Dist: future (==0.15.2)
Requires-Dist: idna (==2.1)
Requires-Dist: ipaddress (==1.0.16)
Requires-Dist: matplotlib (==1.5.3)
Requires-Dist: mmtf-python (==1.0.2)
Requires-Dist: msgpack-python (==0.4.7)
Requires-Dist: numpy (==1.11.1)
Requires-Dist: pandas (==0.19.0)
Requires-Dist: patsy (==0.4.1)
Requires-Dist: py (==1.4.31)
Requires-Dist: pyBigWig (==0.2.8)
Requires-Dist: pyasn1 (==0.1.9)
Requires-Dist: pybedtools (==0.7.8)
Requires-Dist: pycairo (==1.10.0)
Requires-Dist: pycparser (==2.14)
Requires-Dist: pyparsing (==2.1.8)
Requires-Dist: pysam (==0.9.1)
Requires-Dist: pytest (==3.0.3)
Requires-Dist: pytest-cov (==2.3.1)
Requires-Dist: pytest-mpl (==0.5)
Requires-Dist: python-dateutil (==2.5.3)
Requires-Dist: pytz (==2016.6.1)
Requires-Dist: reportlab (==3.3.0)
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.8.4)

==========================================
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://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


Python rewrite of `MoCA0.1.0`_

LICENSE
-------
ISC



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


Current Version
~~~~~~~~~~~~~~~
0.3.1.dev0


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

* bedtools>=2.25.0
* biopython>=1.66
* pandas>=0.18
* scipy>=0.17
* statsmodels>=0.6
* pybigwig>=0.2.8
* seaborn>=0.7.0
* 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 env create -n mocaenv python=2.7
    $ source activate mocaenv
    $ 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:

::

    $ mocacli 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.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.


