Hi everyone,
I have a multisample .vcf file which I've annotated with annovar (gene annotation) and I've used SnpSift CaseControl on it to assign samples as cases and controls. For each variant SnpSift reports (as extra columns in the .vcf) the number of times this variant was found in cases and in controls and it does some statistical testing. However, instead of doing this now on a variant level, I want to check if certain GENES have more variants in the controls than in the samples. This can be again an extra column in the .vcf, or something seperate. I annotated the .vcf with Annovar, so the gene information is present in the .vcf.
Is there a tool (like SnpSift) to do this? Or can I better write my own script? (if so, what would be the easiest way?)
1 answer
Short answer: Use "SnpEff count"
Long answer: Split your file in two: cases.vcf and controls.vcf, then use "SnpEff count", it will tell you the number of variants and bases on each possible genomic interval (gene, transcript, exon, intron, etc.).
http://snpeff.sourceforge.net/SnpEff_manual.html#utils (scroll down to "SnpEff count")
Usage:
java -Xmx4g -jar snpEff.jar count -v hg19 test.vcf > test_count.txt
Log in to answer this question.
-
Tool: Open-source tool for cancer variant interpretation from VCF files
written by Punith •
-
Error occurence after merging files with bcftools: wrong number of fields ?
written by Lukáš •
-
SnpSift case control command for multiple samples
written by kumari.indu31 •
-
Counting allele counts/frequencies from a vcf file with both cases and controls
written by tacrolimus
-
Parsing through VCF annotated by Funcotator
written by abbysue •
-
bam processing before variant calling
written by user230613
-
EPACTS association analysis
written by hellbio
-
EPACTS association analysis
written by hellbio
-
which one would be the better way for RNA-seq Experiment design
written by Manvendra Singh
-
small case-control exome sequencing study
written by User6891
I think you will need to extract that information into a separate table and run an independent statistical test on that.