Hello,
I am new to exome sequencing data anslysis and want to ask questions regarding what to do in my situation. I have spent quite a time to figure this out by myself, but since there's no one around me to direct me, I couldn't get much. I am given vcf files suspected to be from GATK UnifiedGenotyper on case and control samples (A1 and A2 are our case and B1 is our control), namely case[or control].indel.raw.vcf, case[or control].snp.raw.vcf, case[or control].var.raw.vcf. Now, I need to identify 1) rare variants ((SNPs or indels) with frequency less than 0.01% in EXaC or GenomAD) present only in the two cases and not in the control. 2) PolyPhen/SIFT or other scores for the identified rare variants.
My questions are 1) GATK manual says that, since UnifiedGenotyper would produce many false positives, these files need to go through a lot of filtering processes. However, I can't find in the manual what kind of filters I need to apply using what kind of tools. 2) Can you please give me the pipeline to identify rare variants and PolyPhen/SIFT scores from the vcf files?
Thank you very much for your time.
•
written
by
thejustpark
2 answers
Hello thejustpark,
about your question concerning filter false positive variants I recommend reading this blog post first. After that you can try your first steps with this tutorial. But keep in mind that there is no gold standard for doing hard filtering, as it depends on so many things.
The easiest way for you to get population frequencies is to use ensembl's VEP. For missense variants it also report Polyphen/SIFT scores. There are other ways to do this. The term you are looking for is "Variant Annotation". Tools that can do this are for example SnpEff and bcftools.
fin swimmer
Log in to answer this question.
-
WGS based variant analysis for biomarkers
written by halo22
-
filter_vep not properly filtering
written by jeni •
-
Functional Score Annotation Bacteria (VCF)
written by campero189 •
-
What is SIFT-score?
written by mangfu100
-
How to differentiate between protective and damaging variants
written by praj_bonjour •
-
Is It Nonsense To Have Sift/Polyphen2 Scores For Nonsense Variants ?
written by Raony Guimarães
-
Tool: Annovar: Functional Annotation Of Genetic Variants From High-Throughput Sequencing Data
written by Istvan Albert
-
Why Is Gwas Underpowered To Detect Association In Rare Variants?
written by michealsmith
-
Tools To Identify Compound Heterozygous Variants?
written by Ahdf-Lell-Kocks •
-
Missense Mutations Predictions Using Sift And Polyphen-2
written by Kamila001
