SYNOPSIS

macs2diff <-t1 tfile1> [-c1 cfile1] <-t2 tfile2> [-c2 cfile2] [-n name] [-g genomesize] [options]

DESCRIPTION

Example: macs2diff --t1 CTCF_GM12878.bam --c1 Control_GM12878.bam --t2 CTCF_K562.bam --c2 Control_K562.bam -g hs -n testdiff -B -q 0.01

macs2diff -- Differential Analysis for ChIP-Sequencing

OPTIONS

--version

show program's version number and exit

-h, --help

show this help message and exit.

--t1=TFILE1

ChIP-seq treatment file for the first condition. REQUIRED.

--c1=CFILE1

Control file for the first condition. If c1 is missing while c2 is specified, t1 will be paired with c2. At least one of c1 or c2 should be available.

--t2=TFILE2

ChIP-seq treatment file for the second condition. REQUIRED

--c2=CFILE2

Control file for the second condition. If c2 is missing while c1 is specified, t2 will be paired with c1. At least one of c1 or c2 should be available.

-n NAME, --name=NAME

Analysis name, which will be used to generate output file names. DEFAULT: "NA"

-f FORMAT, --format=FORMAT

Format of tag file, "AUTO", "BED" or "ELAND" or "ELANDMULTI" or "ELANDEXPORT" or "SAM" or "BAM" or "BOWTIE". The default AUTO option will let MACS decide which format the file is. Please check the definition in 00README file if you choose ELAND/ELANDMULTI/ELANDEXPORT/SAM/BAM/BOWTIE. DEFAULT: "AUTO"

-g GSIZE, --gsize=GSIZE

Effective genome size. It can be 1.0e+9 or 1000000000, or shortcuts:'hs' for human (2.7e9), 'mm' for mouse (1.87e9), 'ce' for C. elegans (9e7) and 'dm' for fruitfly (1.2e8), Default:hs

-s TSIZE, --tsize=TSIZE

Tag size. This will overide the auto detected tag size. DEFAULT: Not set

--bw=BW

Band width. This value is only used while building the shifting model. DEFAULT: 300

-q QVALUE, --qvalue=QVALUE

Minimum FDR (q-value) cutoff for peak detection. DEFAULT: 0.01

-p PVALUE, --pvalue=PVALUE

Pvalue cutoff for peak detection. When set (e.g. -q 0.05 or -q 1e-5), qvalue cutoff will be ignored. Default is not set.

-m MFOLD, --mfold=MFOLD

Select the regions within MFOLD range of highconfidence enrichment ratio against background to build model. The regions must be lower than upper limit, and higher than the lower limit. DEFAULT:10,30

--nolambda

If True, MACS will use fixed background lambda as local lambda for every peak region. Normally, MACS calculates a dynamic local lambda to reflect the local bias due to potential chromatin structure.

--slocal=SMALLLOCAL

The small nearby region in basepairs to calculate dynamic lambda. This is used to capture the bias near the peak summit region. Invalid if there is no control data. If you set this to 0, MACS will skip slocal lambda calculation. *Note* that MACS will always perform a d-size local lambda calculation. The final local bias should be the maximum of the lambda value from d, slocal, and llocal size windows. DEFAULT: 1000

--llocal=LARGELOCAL

The large nearby region in basepairs to calculate dynamic lambda. This is used to capture the surround bias. If you set this to 0, MACS will skip llocal lambda calculation. *Note* that MACS will always perform a d-size local lambda calculation. The final local bias should be the maximum of the lambda value from d, slocal, and llocal size windows. DEFAULT: 10000.

--auto-bimodal

Whether turn on the auto pair model process. If set, when MACS failed to build paired model, it will use the nomodel settings, the '--shiftsize' parameter to shift and extend each tags. Not to use this automate fixation is a default behavior now. DEFAULT: False

--nomodel

Whether or not to build the shifting model. If True, MACS will not build model. by default it means shifting size = 100, try to set shiftsize to change it. DEFAULT: False

--shiftsize=SHIFTSIZE

The arbitrary shift size in bp. When nomodel is true, MACS will use this value as 1/2 of fragment size. DEFAULT: 100

--keep-dup=KEEPDUPLICATES

It controls the MACS behavior towards duplicate tags at the exact same location -- the same coordination and the same strand. The default 'auto' option makes MACS calculate the maximum tags at the exact same location based on binomal distribution using 1e-5 as pvalue cutoff; and the 'all' option keeps every tags. If an integer is given, at most this number of tags will be kept at the same location. Default: auto

-B, --bdg

Whether or not to save p/qvalue (depending on whether pvalue or qvalue is used) score tracks at every bp into a bedGraph file. DEFAULT: False

--minlen=MINLEN

The minimum length for differential calling. By default, it will be decided as the middle of fragment sizes of treatment 1 andtreatment 2. Must be an integar larger than 0.

-a NPROCESSES

Number of CPUs MACS can use, upto 4. DEFAULT: 1

--verbose=VERBOSE

Set verbose level. 0: only show critical message, 1: show additional warning message, 2: show process information, 3: show debug messages. DEFAULT:2