Coverage per gene
Hey,
I'm working with bedtools trying to find out the coverage of several target regions. I have my bam files and my target regions in a .bed file and I'm running this:
bedtools bamtobed -i WA01.bam | bedtools coverage -a - -b mit.bed > exons.mit.coverage
When I open the exons.mit.coverage file, in the coverage section I just have 0s in the areas that are not matching my target genes and 1s in the areas of my target genes, as I have coverage per read.
My question: How can I obtain coverage per gene and not coverage per read?
I imagine I can filter this out, selecting the reads with a good quality score per gene, but I think this should be already implemented anywhere else that I don't know
Thanks for your help! :)
1 answer
Log in to answer this question.
More posts like this
-
How to obtain zero-based coordinates read depth using bedtools coverage for a specific region?
written by vladimiralan •
-
genome coverage visualization
written by gubrins
-
Regular coverage in target capture sequencing
written by gubrins
-
Coverage of target genes
written by gubrins
-
Target capture sequencing
written by gubrins
-
bedtools bamtobed of mixed paired-end and non-paired datasets
written by tassa.saldi •
-
Difference in BAM coverage between bedtools output and IGV coverage track
written by myology •
-
Depth of Coverage - Calculated Per Probe / Bed Region
written by jnowacki
-
getting reads counts at RefSeq genes
written by sumithrasank75
-
Extract Coverage Track From Bam To Bed File
written by William
see How to calculate average coverage for all genes
I'm also going to check that link, thanks!