get snps from specific individuals from VCF
Hi,
I have a VCF called by GATK with 150 samples, and have around 20k SNPs in total.
I want to filter the VCF, and only keep the 20 samples out of 150, and also only keep the SNPs for these 20 samples. This means I only want to extract the SNPs of 20 samples from a VCF called from 150 samples. I tried with VCF tools --keep, as well as Bcftools view -S, they both can exclude the other 130 samples from the VCF, but the number of SNPs remain exactly the same...
Is there any solution to exclude samples as well as the SNPs called from them?
Many thanks!
• 966 views
•
link
written
by
puddingmeow516 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
GATK SelectVariants behavior with gvcfs
written by Gianfilippo •
-
filter VCF based on FORMAT
written by Matteo Ungaro •
-
Does a merged VCF file (using bcftools merge) contain unique SNPs or INDELs from all merged VCF fi…
written by mohsamir2016 •
-
SNP numbers don't add up to the number of record in a vcf
written by Timotheus •
-
generate FASTA/MSA from VCF file
written by puddingmeow516 •
-
SNPs not filtered from a GATK vcf file
written by evelyn
-
Editing vcf file using bcftools
written by janhuang.cn
-
Variants called differ dramatically between runs
written by Ram
-
SNP calls all seem to have quality of zero
written by novice
-
BCFTools view equivelent of GATK -T SelectVariants -sn sample_x --excludeNonVariants
written by William
If there are multiple sites with missing/unknown GT (
./.) for the samples that remain, you may want to add an expression with with thebcftools view -Sor pipe the output of theview -Sto aview -e 'COUNT(GT="mis" =20)'to exclude sites where all 20 GTs are missing.