Counting allele counts/frequencies from a vcf file with both cases and controls
Dear Biostars members,
I a multi-sample vcf which has been annotated with VEP. The vcf has both cases and controls. I would like to know the allele counts and frequencies for cases vs controls.
I am able to use bcftools to pull the relevant fields but was wondering if there was a way of feeding it a case control label and it giving me a split of the results based from this rather than dividing the vcf into cases and controls and then running the analysis separately?
Many thanks for your time
• 2,461 views
•
link
updated
by
Paolo •
•
written
by
tacrolimus
•
written
by
tacrolimus
1 answer
This is the aim of https://pcingola.github.io/SnpEff/ss_casecontrol/
Log in to answer this question.
More posts like this
-
Counting the specific number of a type of variant per sample in a multi-sample VCF
written by tacrolimus
-
Error occurence after merging files with bcftools: wrong number of fields ?
written by Lukáš •
-
Count number of occurrences an annotation in a multisample vcf split into multiple chunks
written by tacrolimus
-
Extract samples with alternative genotype from a multi-sample VCF
written by tacrolimus
-
Finding unique variants by sample ID in a multisample vcf
written by tacrolimus
-
Faster filtering of vcfs using split-vep over filter_vep
written by tacrolimus
-
Creating a genotype matrix
written by tacrolimus
-
Creating a genotype matrix for BeviMed
written by tacrolimus
-
Filtering large VCFs by allele freqency using bcftools split-vep function
written by tacrolimus
-
Filter a VCF using fields stored in a different VCF
written by Sean
You can create 2 one-sample-ID-per-line files of case IDs and control IDs and use them with the
-Sparameter available in many bcftools sub-tools.