Hi all- I would like to calculate nucleosome dyad densities from my MNase-seq data. Here is where things get complicated: I would first like to take the ratio of two normalized datasets, then calculate the dyad density from there. I have been successful in producing ratio files from two normalized datasets using deepTools bamCompare, but this can only generate BigWig or BedGraph files. I would now like to calculate dyad densities from those files. Is there a way this can be done? Are there alternative programs which would allow me to take the ratio of normalized datasets and produce a bam file as a ratio instead? Producing a BAM file would allow me to use deepTools bamCoverage "--MNase" option. Forgive me if this seems a stupid question, the concept of dyad density calculation is pretty foreign to me.
0 answers
No answers yet.
Log in to answer this question.
-
Nuclesome visualization
written by jetson •
-
Error: unrecognized arguments using deepTools bamCoverage
written by Daria •
-
Scaling vs normalization in bamCompare
written by lorenzo •
-
Split bam files from ATAC-seq PE data into nucleosome free and nucleosomal regions
written by anlo •
-
[ChIP-seq. signal plot] How to merge ChIP-seq. biological replicates for signal profile plotting?
written by chiefcat
-
Help with DeepTools
written by landrjos •
-
Creating coverage files containing paired-end fragments of particular sizes
written by mmmmcandrew
-
Sorting heatmaps by gene expression
written by mmmmcandrew
-
Extracting data from bigWig/BED using threshold?
written by mmmmcandrew
-
MNase-seq / fragment size selection
written by Constantine
You can use bamCoverage with the
--MNasefunction first and thenbigwigComparethe results to presumably end up with your desired bigWig track.Cool, thanks for the input. This is kind of where I was thinking of going. My only concern is that I would prefer to normalize the datasets and take the ratio first, but maybe this can't be helped.
There's no way to go from bigWig->BAM, making a bigWig file is innately a lossy process. I would recommend that you do the following:
bamCompare, but note the scale factors applied to each sample. You might need to use the--verboseoption, I don't recall.bamCoverage --MNaseon each sample to make bigWig files.bigwigComparewith the--scaleFactorsset according to step 1. This should produce more or less exactly what you're trying to get (at least I suspect so).