Question: When using MACS2 with -f BAM, it is possible to have a shift other than 0?
0
gravatar for Lucy
14 months ago by
Lucy40
Lucy40 wrote:

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 • 427 views
ADD COMMENTlink modified 14 months ago by Biostar ♦♦ 20 • written 14 months ago by Lucy40

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

ADD REPLYlink written 14 months ago by ATpoint26k

Yes this is ATAC-seq data.

ADD REPLYlink written 14 months ago by Lucy40

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 REPLYlink modified 14 months ago • written 14 months ago by ATpoint26k

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

ADD REPLYlink written 13 months ago by Lucy40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2517 users visited in the last hour