When using MACS2 with -f BAM, it is possible to have a shift other than 0?
Hi,
In the manual for MACS2, there is the following statement in the information about the --shift option: "Note, you can't set values other than 0 if format is BAMPE or BEDPE for paired-end data. Default is 0"
Is this also the case if the format is set as BAM or can I have different shift sizes with BAM? The command I would like to use to analyse my ATAC-seq data is as follows:
macs2 callpeak <input file> -f BAM --no-model --shift -100 --extsize 200 --keep-dup-all --call-summits
Many thanks, Lucy
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
call peak for all replicates in all conditions using macs2
written by QX •
-
MACS2 not working in paired end mode
written by Filago •
-
Using MACS2 parameters
written by Nishat •
-
Different peak sites from .xls and .narrowPeak, macs2 output
written by sh.kazempour94 •
-
ATAC-seq peak calling?
written by star
-
MACS2 option for ATAC seq?
written by star
-
shift direction in ATAC/DNase data?
written by star
-
ATAC-seq macs2 peak calling bampe mode extension shift
written by nanoide
-
MACS2 produces different results using BAMPE and BEDPE
written by Marius •
-
Help understanding MACS2 --extsize and --shift
written by datascientist28
My understanding is that you do not need any shifting parameters when using
--nomodel. Is this ATAC-seq?Yes this is ATAC-seq data.
Did you solve your problem or do you need further assistance? You raised this question in different posts, that is why I am asking.
Thanks for the message. I think I have decided to go with -f BAMPE and no shifting.