-
Notifications
You must be signed in to change notification settings - Fork 10
/
README
59 lines (42 loc) · 1.81 KB
/
README
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
BSmooth is by Kasper D. Hansen, Ben Langmead, Rafael A. Irizarry
Statistics pipeline is by Kasper D. Hansen
Alignment pipeline is by Ben Langmead
The BSmooth alignment pipeline is licensed under the GPLv3 license. See
`LICENSE' file for details.
Installation
============
BSmooth
-------
To install the BSmooth alignment pipeline, permanently set the BSMOOTH_HOME
environment variable to the main BSmooth directory. I.e., the directory with
'bin', 'example' and 'lib' subdirectories inside.
Building merman
---------------
The merman sources are included with BSmooth. To use merman with BSmooth, you
must first build the merman binary. This requires standard GNU build tools
such as GNU make and g++.
To build merman, cd to the $BSMOOTH_HOME/merman directory and run 'make'.
Finally, add the $BSMOOTH_HOME/merman directory to your PATH.
Installing Bowtie 2
-------------------
Follow the standard instructions for installing those tools; see their
respective manuals for details. Once the tool is installed, make sure that
the directory containing the key binaries ('bowtie2', 'bowtie2-build') is
included in your PATH.
Pipeline flow
=============
The following 4 steps (or possibly 3 steps, depending on aligner) are needed to
get from raw input reads to CpG-level measurement summaries. These summaries
constitute the input to the R pipeline for smoothing and calculating DMRs.
Step 1: Build bisulfite genome index
(not necessary if using Merman)
bswc_bowtie2_index.pl
Step 2: Align
bs_merman_align.pl
bswc_bowtie2_align.pl
Step 3: Sort evidence directory
bsev_sort.pl
Step 4: Tabulate sorted evidence directory
bsev_tabulate.pl
For an example, see example/Makefile, which uses example/sim.pl to create a
small, simulated problem, then goes through the steps listed above.