How many INDELs and SNPs should I expect in vcf file and what the best parameters to filter false positives? · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

How many INDELs and SNPs should I expect in vcf file and what the best parameters to filter false positives?

I am analyzing human samples (WGS and WES) and the total number of variants I get after variant calling are 6 millions (WGS, over 1 million WES). I think the number is too high and there are many false positives. I am using Illumina Dragen Bio-IT plaftorm (which has GATK best practices implemented) for variant calling and I tested it on Coriell NA12878 sample the accuracy is 99.97 for SNPs and 97% for INDELs. Primary QC metrics looks fine and samples are normalized to the sampe coverage 35 for WGS and 100 for WES.

variants wes wgs

• 1,406 views • link

updated 5.4 years ago by Emily

• written 5.4 years ago by esimonova.me •

1 answer

Sounds about right. A diploid human genome is 6 Gb, on average we have one variant every 1 kb, so 6 million variants per sample is to be expected.

• 0 views • link

written 5.4 years ago by Emily

Log in to answer this question.

More posts like this

Seeking Advice on Multisample Variant Calling Pipeline with BWA-MEME and GATK DRAGEN Mode
written 2.0 years ago by George •

Hello everyone, I am in the process of building a pipeline for multisample variant calling and have some questions regarding tool usage and their safety. …
Variant quality scores with different variant caller: HaplotypeCaller vs DeepVariant. How to interp…
written 4.8 years ago by Jordi •

Hi, I am trying to compare variant calling outputs of GATK's HaplotypeCaller and DeepVariant. Their raw output is very different; for example, in a WGS …
VG Toolkit for human WGS analysis
written 4.9 years ago by venkat •

Hello! I have been developing a computational pipeline with VG Toolkit for human WGS analysis. I have built a genome graph with GRCH38 and variants …
Needs recommendation - Variant Call from RNA seq data of tumor-normal pairs.
written 5.1 years ago by field654 •

Dear folks, I understand that the WGS and WES are more frequently applied for variant call analysis. I've seen a few reports using RNA-seq for …
Freebayes - too many variants produced on WES sample
written 5.3 years ago by esimonova.me •

I am currently testing different variant callers on my sample (NA12878) and the majority of other variant callers generate 30,000 variants in this sample, while …
Precision of variant calling (INDELs)
written 5.3 years ago by esimonova.me •

The same sample (NA12878) was sequenced as WES and WGS, and analysed by the same pipeline (in case of WES it was restricted to the …
GATK VQSR filter
written 5.4 years ago by esimonova.me •

Hi, I have sequenced (WES) some paired blood and saliva samples along with NA12878 standard sample. I have used a bioinformatics pipeline to generate raw …
Optimal Number of SNPs in NA12878 WES Data
written 5.6 years ago by hkarakurt

Hello everyone, I downloaded one of the NA12878 WES data, SRR098401, from NCBI FTP server with sratoolkit. I used BWA-MEM to align it to hg19 …
Tool: elPrep 4.0.0, a high-performance drop-in replacement tool for GATK4/Picard/SAMtools for proce…
written 7.8 years ago by Charlotte.Herzeel

Dear colleagues, We are happy to announce the release of elPrep 4.0.0, an open-source, drop-in replacement tool for GATK4/Picard/SAMtools for preparing SAM/BAM files for variant …
Variant analysis on 1 or 2 samples: are the final steps of the GATK pipeline unneccessary?
written 8.2 years ago by BioinfGuru

Hello all, I'm a bit confused as to what steps are necessary, and what steps are not going to add much benefit. I have 2 …