Hi Everyone, I am trying to go from a bam file to a bedgraph file where I have fixed intervals of 100bp.
chr1 0 100 4
chr1 100 200 5
chr1 200 300 6
...
I have used genomeCoverageBed to compute the coverage, but can I use it to output counts in fixed window sizes of 100bp?? I would like to use that because it has scaling option and that would make one less step downstream.
If not what can I use to include scaling factor and fixed window size? igvtools count?? or deeptools??
Thanks
•
written
by
Varun Gupta
1 answer
If you make a pre-binned bed for your genome, then plain old coverageBed should take that + bam file and give you bin counts or means. But, coverageBed does not take scaling factors. It writes to stdout, you could pipe it through awk to scale it.
It is also fairly quick in R, if you want to use R -- the code is a little complicated to set up, but you only have to do it once! -- see the rtracklayer package for reading BAMs into GRanges objects, and the GenomicRanges package for tileGenome() and viewApply() functions.
Log in to answer this question.
-
On applying spike-in-derived scaling factors to ChIP-seq data: adjust for background signal (input)…
written by kalavattam •
-
Bedtools coverage -hist "all" in chr column
written by Smriti •
-
partition single bed file into non-overlapping regions
written by jrose •
-
Genome coverage within specified regions
written by ilovesuperheroes1993 •
-
Output bedGraph from genomeCoverageBed is missing "chr"
written by m98
-
genomeCoverageBed for per nucleotide coverage
written by Varun Gupta
-
Merge bed files, keep the peak count
written by rioualen
-
How do I replace a value used in script with a range of values in a file?
written by josephmccarter1
-
How to identify consecutive (adjacent) genomic intervals
written by A. Domingues
-
Split bed file per sequence length
written by cg2827 •
