CpelAsm is a julia package for haplotype allele-specific methylation, based on the method in [CITE], that draws ideas from statistical physics and information theory to identify allele-specific methylation events at the haplotype level.
- julia > v1.0.0
- git.
CpelAsm
and dependencies can be installed via the following command:
julia -e 'using Pkg; Pkg.add(PackageSpec(url="https://github.com/jordiabante/CpelAsm.jl.git"))'
In a julia
session run
(v1.0) pkg> test CpelAsm
The package includes a small toy example for illustrative purposes.
The example consists of two alleles a1
and a2
. The former has a
mean-methylation level (MML) of 0.4, while that of the later is 0.6.
Nevertheless, given the symmetry of the problem, both alleles have
the same Shannon entropy. Thus differential analysis only identifies
differences in terms of MML. The output bedGraph files can be found
in out_path
. To run the toy example run the following commands in
a julia
's REPL:
using CpelAsm
dir = "/path/to/CpelAsm.jl/test/"
b1 = "$(dir)/bam/example.a1.bam"
b2 = "$(dir)/bam/example.a2.bam"
fa = "$(dir)/fasta/example.fa"
vcf = "$(dir)/vcf/example.vcf"
out = "$(dir)/out/"
run_analysis(b1,b2,b1,vcf,fa,out;g_max=25,cov_ths=3,cov_b=20.0,win_exp=10,n_null=100,n_max=10)
- Jordi Abante
This project is licensed under the MIT License - see the LICENSE.md file for details.
[CITE]