Questions tagged [deeptools]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
1
vote
0answers
24 views

Metagene analysis with deepTools; 5UTR-CDS-3UTR and Exon-Intron-Exon

I have some PCR-CLIP data. PAR-CLIP identifies targets for RNA binding proteins. I am interested to build metagene profile so that transcript wide binding pattern can be compared. What I am mostly ...
0
votes
1answer
43 views

Pol_II genome wide occupancy analysis

I would like to know how to perform this global Pol II occupancy which is shown here: Paper Figure 1G So in this figure they are analysing the Pol II occupancy in Control,PAF1 KD and LEO1. I have ...
2
votes
2answers
282 views

How to convert BAM file to bigWig in Python?

I am manipulating some BAM files using pysam package which seems to be very fast and handy. However, I ran into problem when I am trying to generate bigWig files ...
1
vote
1answer
75 views

merge the matrix from computematrix of deeptools

I have several matrix from computematrix of deeptools. I need to merge two of them using "computeMatrixOperations cbind -m input1.mat.gz input2.mat.gz -o output.mat.gz" but I am running to error "/...
1
vote
1answer
43 views

Clarification regarding (deepTools) bamCoverage output

Im using bamCoverage on a file and ive set the bin size to 1 but also normalizing using RPKM for defined regions using a bedgraph file. The output from bamCoverage however is in binned segments each ...
0
votes
0answers
31 views

Setting scale of profile plot from zero for chip seq profile plot

Figure from this paper which shows a profile plot in figure4 first figure. .This i guess is made using deeptools perhaps. Meanwhile When i tried to do similar plot the scale doesn't start from 0. So ...
2
votes
1answer
24 views

How to calculate the number of peaks that are upstream/downstream of some other peaks

I have 3 histone marks,I have used Macs2 for peak calling and diffBind to analyze the peaks. I was wondering if you know any way to calculate the peak numbers of one specific histone mark that occur ...
1
vote
1answer
138 views

plot correlation between several bam files

I have several mab files from medip-seq. they supposed to be replicates and I have to draw their correlation. First I used multiBamSummary from deeptools to merge all the bam files based on bi/bed. ...
0
votes
1answer
54 views

multiBamSummary output chromosome region annotate issue

I was looking how to take readcount from ATAC seq data , this post which suggested to use bedtools and featurecounts to take out readcounts .But multiBamSummary seems to be pretty straight forward. ...
2
votes
1answer
107 views

RPGC normalisation creates artefacts at centromere

I am studying ChIP-Seq data in HeLa cells and I've started using the RPGC normalisation of deepTool's bamCoverage. MACS2 also uses this normalisation for its peak calling. I am seeing a large number ...
1
vote
0answers
20 views

scaling algorithm of deeptools

I am using deeptools to generate a metagene plot. I am wondering how the deeptools scale genes if the size of genes is shorter than bin size. For example, if my bin size is 1000 bp and the gene is ...
2
votes
1answer
171 views

Different outputs of computeMatrix depending on binsize

I'd like to output a figure in +/-1kb around a center using computeMatrix as following command: ...
5
votes
1answer
163 views

Dealing with absence of coverage

In my snakemake workflow, I use deeptools bamCoverage (in a make_bigwig rule), then ...
3
votes
1answer
91 views

IGV faster loading remote bigWig files

I noticed remote bigWig files (http URL) take sometimes about 1 minute to load on IGV. Once they loaded, jumping from one region to another can also take several seconds. My understanding is that ...
6
votes
1answer
658 views

What unit do I get on the y-axis of a metagene profile plot?

I start with a sorted and indexed bam file ("mapped.bam") representing the mapping of small reads on a reference genome, and a bed file ("genes.bed") containing the coordinates of a set of features of ...