Has anyone used Sicer to peak call using bam files and bed files. I would appreciate it if someone could help me get the command to do peak calling using treatment vs input bam files. Also, how does sicer use bed files to peak call, i thought peak calling was produced bed files. Would appreciate any help.
•
written
by
mropri
1 answer
Have you read the manual? It gives very clear commands as for how to run with a control file (or without):
sicer -t treatment.bam -c control.bam -s hg38 -w 200 -rt 1 -f 150 -egf 0.74 -fdr 0.01 -g 600 -e 1000
Those are the default parameters. You can supply the reads in either bam or BED format. Having aligned reads in BED format isn't too uncommon.
The output files are in BED-like format, you can read more about them here.
So are you still running into issues or were you just a bit confused about the input/output?
Log in to answer this question.
-
DiffBind analysis with ENCODE data using R and .bed/.bam files
written by stalo.lili •
-
Chipseq data peak calling issue
written by Adyasha •
-
What is the advantage of having BAM files per sample compared to one BAM file containing all sample…
written by L_to_the_m •
-
Collapse two or thee bed files into one bed file using bedtools merge
written by mropri
-
Merging multiple fastq files and having IndexError when calling peaks using SICER
written by dcheng1 •
-
Using genome tools to convert bed files to gff3 files
written by mropri
-
MACS2 MODEL VS NO MODEL FOR CHIRP SEQ
written by prathyushakonda •
-
Data size for ChIPseq and RNAseq data
written by adelest.ad •
-
peak calling of ChIP-seq
written by Ben
-
Why is an Input file used as a control during MACS peak calling?
written by christopher.mcnair •
