I have few quick question regarding macs14?
How one can treat the replicates? Do I have to run separately and then take average from the final peaks or is there any way to include it in one run?
macs14 generates two wig files under "--single-profile" option; one for control and one for treated. which file to use for plotting final peak profile? Do I have to subtract treated - control manually or does it give the final peak wig files?
Thanks in advance
By control I mean input here. Can I use treated wig/bw directly to see my enrichment of do I have to do manual subtraction of treated wig - control (input) wig file?
"I would get the bigwigs with -log10(p-value), to show enrichment over inputs." will this file be generated or do I have to use some other file to generate this?
I am asking this because my folder NAME_MACS_wiggle has two separate folders naming "control" and "treated" and each of them having one wig file; NAME_control_afterfiting_all.wig.gz, NAME_treat_afterfiting_all.wig.gz respectively.
You can use your treated file directly but its better to show the input as well so that it will be more convincing. Are you talking about ucsc tracks or a profile plot over promoters/enhancers ?
You can run MACS2 sub command like this:
macs2 bdgcmp -t treat.pileup.bdg -c local_lambda.bdg -m ppois -o pvalue.bdgwhich contains -log10(p-values) for each base pair which can be converted to bigwig (bedGraph to bigWig)
This is because its good to show both treat and input so that its more convincing that your peaks are real.
Thanks @geek_y for such detail explanation. Al thought have you ever encounter the following error while running the macs14. "error occurred: Invalid variableStep wig file format in line starting with "variableStep chrom=chr3 span=10"
what could be the reason for this?
MACS doesn't work with bigwig. It works with bed graph files.