What allele frequency threshold do you recommend?
Hi all!
I have a result from Varscan mpileup2snp and mpilepu2indel with default parameters including variant allele frequency =0.01 (1%). (--min-var-freq)
I am wondering, what is the recommendation for allele frequency threshold? This 1% looks very low and creates results with unwanted records.
Thanks in advance!
Best,
Agata
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Missing SAM header when piping the samtools mpileup output to a VarScan call
written by melissachua90 •
-
Varscan output variant file - is it the "true" file?
written by melissachua90 •
-
varscan can't call a SNV
written by tan •
-
annovar mpileup2cns and mpileup2indel give different results
written by Assa Yeroslaviz
-
Set minimum variant frequency while calling variants with samtools
written by kspata
-
varscan not taking opstions
written by fabbri.marco •
-
The different results for the same tumor sample between VarScan2 somatic and mpileup2snp
written by rayliu0115 •
-
Different output from pileup and mpileup commands VarScan
written by agata88
-
Varscan mpileup2snp parameters
written by AW
-
varscan2 mpileup2snp output contains indels
written by AW
Is this for somatic or germline variants in a diploid organism?
Its germline, human.
So theoretically the variant frequency for a heterozygous variant would be ~50%, right? Then indeed 1% is very unlikely. Obviously it will rarely be exactly 50%, but with sufficient coverage it should be quite close to that. It's a trade-off between sensitivity and specificity, but I don't think it's unreasonable to expect at least 25%.
Thanks, I agree, especially when we have probe, WES and WGS sequencing where the de-duplication should be performed. I thought that maybe even higher in this situation, 35%?. In primer targeted sequencing (multipleks) I think it should stay 1%, do you agree?
You can indeed be more stringent and as such reduce false positives but you'll miss variants which are by chance only represented in 34% of the reads. It's a trade-off.
Multiplex PCR-based target enrichment prior to NGS? I think you'll commonly get high coverage there and can be more stringent. Unless you have SNPs in your primer I don't think you should expect a divergent frequency of variant alleles.