Hi,
I am working with ChIP-Seq data and am interested in doing quantile normalization.
I created bedgraph files from my bam files using bamCoverage with a window size of 50 bp so that I get the number of reads for every 50 bp of genome. However I noticed that every bam file has a different number of lines (even after I split up the intervals that bamCoverage automatically merged) so I cannot use this for quantile normalization.
Now I am using bedtools multicov with a bed file I created with 50 bp windows across the genome (using bedtools makewindows) and my bam files to get read counts for all of my bam files using the same intervals but it is taking a long time to run.
I am wondering if there is a simpler, faster way to use quantile normalization with bam files for ChIP-Seq analysis.
Any help would be appreciated,
Thank you,
ej
1 answer
By default bamCoverage does not report bins with zero counts - this may be the reason why you are getting different number of lines reported. Instead of running bedtools multicov, you could run bedtools coverage on each of your samples in parallel, then merge them into one large count matrix.
Log in to answer this question.
-
deepTools bamCompare - RPGC normalisation not supported?
written by nach.mcnach •
-
Bedtools Fisher Genome File
written by megannj •
-
Mapping reads from a BAM file to a custom BED file.
written by Shaurya Jauhari •
-
what does the bed file in bedtools multicov command refer to ?
written by rania.hamdy1 •
-
Normalizing ChIP-seq and ATAC-seq data?
written by star
-
bedtools multicov displays only the last BED feature
written by ThePresident
-
how to calculate coverage in bam file over a specific bed file
written by Lila M
-
Calling peaks with macs2 if not all replicates have a corresponding input-file
written by jklopp •
-
Bedtools multicover how to submit a bam list
written by benoit.tessoulin •
-
average coverage, chip-seq
written by marina-orlova
