Question: how to calculate coverage in bam file over a specific bed file
gravatar for Lila M
2.1 years ago by
Lila M 470
Lila M 470 wrote:

Hi every body, I have bam files (in different conditions) and a bed file with specific genes of interest. I would like to calculate the coverage for each bam file over my gene list. What I am doing is, for each condition:

bedtools multicov -bams *cond1.bam -bed exon.bed > multicov_condition1_exon
bedtools multicov -bams *cond2.bam -bed exon.bed > multicov_condition2_exon
bedtools multicov -bams *cond3.bam -bed exon.bed > multicov_condition3_exon

The final idea is to identified in which condition the coverage is higher. Is that a good option to analyze my data?


ADD COMMENTlink modified 2.1 years ago by Devon Ryan90k • written 2.1 years ago by Lila M 470
gravatar for Devon Ryan
2.1 years ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

From deepTools: multiBamSummary BED-file --bamfiles *.bam --BED exon.bed -o multicov.npz. The result in a compressed numpy matrix that you can load with numpy.load in python.

ADD COMMENTlink written 2.1 years ago by Devon Ryan90k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1053 users visited in the last hour