When using MACS2 with -f BAM, it is possible to have a shift other than 0?

0

Entering edit mode

6.1 years ago

Lucy ▴ 140

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

macs2 macs • 1.9k views

ADD COMMENT • link updated 6.1 years ago by Biostar 20 • written 6.1 years ago by Lucy ▴ 140

0

Entering edit mode

My understanding is that you do not need any shifting parameters when using --nomodel. Is this ATAC-seq?

ADD REPLY • link 6.1 years ago by ATpoint 85k

0

Entering edit mode

Yes this is ATAC-seq data.

ADD REPLY • link 6.1 years ago by Lucy ▴ 140

0

Entering edit mode

Did you solve your problem or do you need further assistance? You raised this question in different posts, that is why I am asking.

ADD REPLY • link 6.1 years ago by ATpoint 85k

0

Entering edit mode

Thanks for the message. I think I have decided to go with -f BAMPE and no shifting.

ADD REPLY • link 6.0 years ago by Lucy ▴ 140

Login before adding your answer.