How to sub-sample .bam file to target coverage?
Hi,
I have .bam files for blood samples from which I extract signals. Let's assume my samples have a 30x coverage and my signal of interest is nicely detectable. To see how well the signal is detectable with lower coverage samples, I want to downsample to 20x, 10x, etc. Is there an easy way to do this?
So far I had a closer look at samtools view -s, but it takes only the target fraction and not a target coverage. Is there a program that lets me specify the target coverage that generalizes well for multiple samples with different start coverage?
• 1,968 views
•
link
written
by
hivemind •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Variant counts decrease across cumulative coverage tiers after genotype masking
written by itsmesmb •
-
How to calculate depth of coverage and applying read depth filter
written by Sd •
-
Dealing with variable sequencing depth.
written by a •
-
Any Tools to calculate specific gene coverage from BAM file?
written by koay •
-
Capping/downsampling read coverage in bam file
written by hivemind •
-
Lower Assembly Quality Despite Higher Coverage?
written by AXR •
-
Extract the overlap of well-covered regions across multiple samples
written by DVA
-
Downsampling bam file to 2x coverage
written by Jautis
-
Down-sample to one read per SNP loci per sample
written by aylward.megan •
-
Calculate Coverage From Bam File
written by meredith19 •
If you want to downsample from 30x to 20x then then just do
samtools --subsample 0.66 .....This assumes that unaligned reads are not in the BAM file? If they are then it may be best to filter those out first.
someXcoverage is a fluid term and subject to interpretation so it could mean different things to different people.I know... If I had one sample or even five samples this would be my solution. I have hundreds of samples, all with different coverage and I was wondering whether there is a simpler solution that for example caps the max coverage to a certain value so that I don't have to calculate the fraction for each sample.