I have some PE ATAC-seq data and i like to do peak calling using MACS2. I found that If I like use BAM file I have to use -f BAMPE
option otherwise I have to change BAM file to BED and then use -f BED -- nomodel --shift -100 --extsize 200
. Also, I found in ENDCODE https://docs.google.com/document/d/1f0Cm4vRyDQDu0bMehHD7P7KOMxTOP-HiNoIvL1VcBt8/edit
part 2a, that use bed file but there 2 files BEDPE
and TAG file
-tagAlign file ${FINAL_TA_FILE}
-BEDPE file (with read pairs on each line) ${FINAL_BEDPE_FILE}
-Subsampled tagAlign file for CC analysis ${SUBSAMPLED_TA_FILE}
that I am not sure which one should be use as Input of MACS2. should I use TAG file as Input?