How to calculate depth of coverage and applying read depth filter
I want to calculate coverage for a CRAM file. I want to see how the coverage varies by chromosome in different genomic regions and apply filters for 1X, 10X, 20X, and 30X to compare the results. For example, applying a 10X filter means calculating and providing genomic regions with at least 20 read depth, excluding regions with below 20 read depth, and then providing coverage results. Is there a way to do this using samtools depth, samtools coverage, or any other tools?
1 answer
You can do that with https://github.com/brentp/mosdepth as a starting point and then postprocess the output.
Log in to answer this question.
More posts like this
-
Pacbio Longread Depth of coverage
written by hellbio
-
Calculating coverage/depth for nanopore dRNAseq reads with reference transcriptome/genome
written by kenneditodd •
-
Any Tools to calculate specific gene coverage from BAM file?
written by koay •
-
allele frequency from different read depths
written by rthapa •
-
subset of paired end reads from bam file or fastq file of different depth
written by rthapa •
-
Lower Assembly Quality Despite Higher Coverage?
written by AXR •
-
Depth of Coverage without Read Overlap
written by drkennetz
-
Coverage analysis of NGS data by Depth
written by nilakshafreezon
-
coverage per base using samtools depth and mpileup
written by tonja.r
-
Calculate Coverage From Bam File
written by meredith19 •
