Sounds cool. I will definitely try. Thank you, Alex! Xiaoyong
• 0 views
•
link
written
by
xiaoyonf •
How to create a new bigwig file defined by a subset of bed file?
Hello, I have a ChIP-seq bigwig file and a bed file that is a subset of the binding peaks defined from the bigwig file. I want to create a new (sub) bigwig file to use the subset of bed file. How can I do it? Many thanks! Xiaoyong
• 5,041 views
•
link
updated
by
Alex Reynolds
•
written
by
xiaoyonf •
•
written
by
xiaoyonf •
1 answer
Convert from bigwig to wiggle:
$ bigWigToWig data.bw data.wig
Convert from wiggle to bed:
$ wig2bed < data.wig > data.bed
Filter the bed file based on your peak subset, and make a bedgraph file from the subset:
$ sort-bed peak_subset.unsorted.bed > peak_subset.bed
$ bedops --element-of 1 data.bed peak_subset.bed | awk 'BEGIN{ OFS="\t"; }{ print $1, $2, $3, $5 }' > data_subset.bedgraph
Convert the bedgraph file to a bigwig file, e.g. for reference genome hg19:
$ fetchChromSizes hg19 > hg19.sizes
$ bedGraphToBigWig data_subset.bedgraph hg19.sizes data_subset.bw
• 0 views
•
link
written
by
Alex Reynolds
Log in to answer this question.
More posts like this
-
Normalizing chip-seq signals with different number of peaks
written by Rajendra KC •
-
ChIP-seq heatmap input bed file
written by buffealo •
-
Motif match chip-seq peaks
written by ta_awwad
-
How to calculate the input-subtract signal for a specific enhancer region
written by biofandd •
-
How to do differential peak analysis after scaling BAM files
written by xiaoyonf •
-
A possible script request for a task probably cannot be done in Excel
written by xiaoyonf •
-
Which tool can plot the enrichment of specified motif across ChIP-seq bed files?
written by xiaoyonf •
-
Bedops element or merge operation_which one is better to use in predicting regulated genes?
written by xiaoyonf •
-
bigwig to peak
written by Nicolas Rosewick
-
bigwig : peak distance from specific genomic region
written by Nicolas Rosewick
