Questions tagged [macs2]
The macs2 tag has no usage guidance.
15
questions
1
vote
1
answer
72
views
What is the required input peak file format for IDR (Irreproducible Discovery Rate)
I have managed to install IDR (https://github.com/nboley/idr) correctly (idr --version returns 2.0.3), but I can't get it to run on my peak files (*.bed) without crashing.
My peak files are generated ...
- 323
1
vote
0
answers
16
views
Questions about ChIP-seq reads and MACS
I had a few questions raised by the below diagram in relation to how ChIP-seq reads look like (source). I would be very grateful for your insight on them!
From the below diagram, it seems that, in a ...
- 111
1
vote
2
answers
68
views
How to extract DNA sequence from browser track input files (BigWig, bed etc.) files
A BigWig file for something like Chip-seq, or ATAC data has the following format for each line i:
chrom_i - start_i - end_i - value_i
I need to convert for each line the (start_i -> end_i) ...
- 323
0
votes
2
answers
50
views
Separating peaks of chip-seq with specific length
My data contain several chip-seq results. I have the peaks called by MACS2.I wanted to only look at those peaks that their size is e.g 500bp to 1000bp. How can I separate those peaks efficiently?
I ...
- 115
0
votes
3
answers
196
views
differential analysis of chip-seq data
I have several sets of chip-seq data. I called the peaks using Macs2. I am pretty new to the field and I will appreciate any help. I wanted to annotate the peaks and see which peaks are shared between ...
- 115
0
votes
2
answers
373
views
How to make Venn diagram for the Macs peak calling output of two data sets?
I have two outputs of macs2 peak calling for two of my data sets. I wanted to plot the Venn diagram to see how many peaks are shared. I mainly work on Unix/Linux. Do you know any way that I can have ...
- 115
1
vote
2
answers
5k
views
merging two/or more bed file into one bed file
I am trying to merge two bed files (more in future) to one. my bed files are something like :
.
I need to merge them in a way to have the shared chromosome location.
Is there a way to do that ?
- 115
1
vote
1
answer
767
views
How to get the intersection of peaks after peak calling using MACS2?
I have following 5 narrow peak files after peak calling.
K14_peaks.narrowPeak
K15_peaks.narrowPeak
K16_peaks.narrowPeak
K3_peaks.narrowPeak
K8_peaks.narrowPeak
I ...
- 341
0
votes
1
answer
92
views
1
vote
1
answer
243
views
When running ChIPQC bioconductor package after peak calling, should I use sorted BAM files or unsorted BAM files?
I have called peaks using MACS2 and now I'm trying to run ChIPQC bioconductor package. Can someone please tell me should I use sorted BAM files or unsorted BAM files? I have created sorted BAM files ...
- 341
1
vote
1
answer
2k
views
How to visualize called narrowPeak files in UCSC Genome browser or IGV?
I have called peaks using MACS2. Then I got a narrowPeak file like this.
...
- 341
0
votes
1
answer
141
views
ATAC-seq macs2 peak splitting in sliding windows
This question has also been asked on Biostars
I used macs2 to call peaks for atac-seq data. now my goal is to split the peaks into 50 bp windows with 25 bp steps and then calculate the Tn5 integration ...
- 11
1
vote
2
answers
1k
views
Convert paired-end BAM into a single-end BAM and keep all the reads
I want to call peaks using MACS2 on my paired-end ATAC-seq BAM and I am interested in cutting sites. So I make use of the -BAM option (opposed to the -BAMPE) option. However this throws away all the ...
- 111
1
vote
2
answers
821
views
Biological replicates on Chip-seq Transcription factor data
I have 2 biological replicates for chip-seq transcription factor data and I have applied macs2 peak calling for each replicate separately.
How can i make use of the biological replicates to extract ...
- 11
1
vote
1
answer
184
views
Score predefined ChIP-seq peaks with MACS2 or equivalent
I have performed peak calling on a number of separate ChIP-seq experiments and would like to harmonise the peaks from each of the experiments in order to convert my data into a matrix for further ...
- 2,133