Minimum and maximum depth filter - vcftools

Hello,

I would like to implement the minimum and maximum depth filters in vcftools using the --minDP and --maxDP flags. Is there a standard way to calculate the thresholds for such filters?

Thanks in advance!

Matteo

vcftools

• 1,218 views • link

updated 2.0 years ago by Pierre Lindenbaum

• written 4.1 years ago by Matteo •

vcftools is deprecated use bcftools -i 'INFO/DP>1 && INFO/DP<1000'

• 0 views • link

written 4.1 years ago by Pierre Lindenbaum

0 answers

No answers yet.

More posts like this

Estimating timing of admixture
written 3.5 years ago by Matteo •

Hi everyone, I would like to know if there is a way to estimate the number of generations since the admixture event given a hybrid …
Detecting Assortative Mating
written 3.7 years ago by Matteo •

Hi everyone, I would like to ask whether there is a method to detect assortative mating (or signatures of) in a wild population from genetic …
Filter SNPs for DP using mean + 2*SD
written 6.3 years ago by QPaps04 •

Hello all, I am trying to filter a SNP file (.vcf) for max depth. Instead of choosing a random number it has been suggested that …
How to replicate vcftool's sliding window SNP density?
written 6.9 years ago by serpalma.v •

Hello! vcftools can calculate statistics such as pi and Fst using sliding windows. I've noticed that windows in which the alternative allele is fixed are …
Calculating effective coverage/callable genome for variant calling from a bam file
written 7.0 years ago by Rubal

We have mapped whole-genome sequence data in bam format. We are testing a variety of variant callers to identify mutant sites. The variant callers we …
Depth of coverage
written 7.4 years ago by jaafari.omid

Hello dears all, Excuse me I think I have got stuck in a misunderstanding about depth. Could you please let me know what is the …
filter multisample vcf based on per-sample average read depth
written 8.7 years ago by agathejouet •

Hi, I have a multisample vcf that I would like to filter based on read depth. More precisely, I would like to modify the vcf …
VCFtools --minDP flags not working
written 9.7 years ago by beausoleilmo

I have a VCF where I'm applying filters. I want to filter by minDP and maxDP. vcftools --vcf output.vcf --minDP 4 --maxDP 100 --recode --out …
vcftools does not filter by GQ
written 10.2 years ago by AP

Hello, I am trying to filter based on GQ < 15. I do the following: ``` vcftools --vcf infile.vcf --minGQ 15 --recode --out filtered ``` …
Calculate Coverage From Bam File
written 13.3 years ago by meredith19 •

Is there a way to calculate 5x, 10x, 20x, 30x coverage from a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file ? EDIT: I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file and …