Calling peaks with MACS following Bedtools

I ran the following bedtools job commanding the following:

bedtools bamtobed -I bam_file.bam

Now, how do I get from this to calling peaks in MACS within Signac (I have generated a _summits.bed, peak_annotation.tsv, _peaks.xls, and .narrowPeak file). Please provide code as the tutorial is not clear enough for me.

bedtools scatac-seq signac

• 653 views • link

updated 3.1 years ago by Ram

• written 3.8 years ago by bioinformatics.girl •

0 answers

No answers yet.

More posts like this

Diffbind or DESEQ2 Chipseq data
written 10 months ago by Irene •

Hi everyone, I have a methodological question regarding differential analysis of ChIP-seq data. From what I understand, the general workflow is to first define peak …
Joint RNA and ATAC analysis: 10x multiomic - Fragment Path
written 2.9 years ago by cgp09741 •

I am trying to follow the tutorial for Signac 1.10.0... I am getting an error when it comes to the part of the code "create …
.BAM to .BED to Coverage Plots
written 3.8 years ago by bioinformatics.girl •

Hello all, I have been stumped for a good number of weeks concerning scATAC. For analysis, I ran CellRanger to generate the /outs file which …
MACS2
written 5.4 years ago by snehu.sambare •

I ran MACS2 on a paired-end CTCF bam file with the following options: > -t treatment -c control -f BAMPE -g 2.9e+9 -p 1e-5 -B. …
MACS: peak calling from snATAC data
written 5.4 years ago by abiuma •

Hi, I followed this tutorial to create the object Vis (https://satijalab.org/signac/articles/pbmc_vignette.html). From the object created here, I have been trying the peak calling from signac …
Generalizing snakemake MACS2 code - WildcardError Wildcards in input files cannot be determined fro…
written 6.3 years ago by msimmer92

I am practicing Snakemake. Wanted to make a MACS2 callpeak code, and ran into problems when I wanted to add a level of complexity (from …
Running two bedtools commands in series produces error "has non positional records"
written 8.1 years ago by Lina F

Hello all, I have three samples, generated Illumina PE reads for each of them, and then mapped the reads against the reference to produce three …
Creating reproducible peak regions from multiple ChIP-seq studies
written 9.2 years ago by James Ashmore

I have called peaks across 10 different ChIP-seq studies. The same protein was assayed in each study. I am looking to create a 'master' peak …
Mpileup Segmentation Fault
written 13.9 years ago by AW

<p>I am interested in calling SNPs with mileup using de novo contigs generated from a Trinity run and reads aligned with Rsem.</p> <p>I ran the …
Macs 1.4 Reports Many Negative P-Values
written 14.0 years ago by Konrad Rudolph •

<p>I am calling peaks on ChIP-seq data (BED files mapped with <a href='http://bio-bwa.sourceforge.net/'>BWA</a>) using MACS 1.4. I observed that there were a large number of …