/ECE

Finding Extremely Conserved Elements (ECE) from multiple species alignments

Primary LanguageC

+============================================+
+ Finding ECEs in .wigFix phyloP files
+============================================+

Usage:
	python ECE_finding_phyloP.py -f <filename> [-p <phyloP>] [-m <min_ece_length>] [-o <output_filename>]

Only the input filename (should be wigFix format) is absolutely required.

Type:

	python ECE_finding_phyloP.py --help

to see the param defaults and usage.

An example usage is:
    python ECE_finding_phyloP.py -f example.wigFix -m 50 -p 3.1 -o example.ece.txt

 which finds ECEs with length >= 50 and avg. phyloP score >= 3.1 

+=================================+
+ Additional note
+=================================+

The core ECE algorithm is implemented in a pure Python version (findECE_pure.py)
and a Cython version (findECE.pyx). The latter reduces the run time by 50%.

In most cases the Cython module (findECE.so) should be directly usable.
In that case nothing needs to be done for setting up the scripts.

To test that the Cython version works, invoke the python interpreter and type:
>>> import findECE
>>> findECE.test()

However if it looks like it needs to be recompiled, do:

	make clean
	make


The main ECE script (ECE_finding_phyloP.py) will try to use the Cython version
whenever possible. When it fails, it falls back to pure python version.