Off topic:ATAC-seq peak calling?
0
0
Entering edit mode
5.2 years ago
star ▴ 350

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 BEDPEand 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?

ChIP-Seq macs2 • 1.8k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2581 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6