VCFtools filter missing genotypes conditional on a factor (population, species, etc.) · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

VCFtools filter missing genotypes conditional on a factor (population, species, etc.)

I wonder if it's possible to filter missing genotype values with VCFtools conditional on a factor. As you can see on VCFtools website:

GENOTYPE VALUE FILTERING

--max-missing <float>

But I have a file that contains individuals coming from different populations. I would like to filter missing genotypes per population. As you can see in this image, the titles circled in red are the individuals. The last 3 letters (ful and for) are the species name. I would like this top be my conditional factor.

If VCFtools is not able to do this, how can I execute this?

enter image description here

vcftools filter ngs rad

• 2,656 views • link

written 9.6 years ago by beausoleilmo

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Problem after filtering vcf files
written 3.8 years ago by yoser4 •

Hello everyone I have a vcf file (snp) of 99 samples, after filtering with vcftools, the code is as follows: ``` vcftools --gzvcf my.vcf.gz \ …
Calculating Fst from low coverage data using genotype likelihoods
written 5.3 years ago by Jautis

Hello, I have low coverage (~1x) sequencing data from ~50 individuals from two populations and would like to calculate Fst for each biallelic variant. I …
How to order Admixture (ancestry) plot in R
written 9.6 years ago by beausoleilmo

I created an ancestry plot in R from an Admixture analysis. But my order in is very bizarre. This is the command I used for …
How common are Deletions and Indels?
written 9.6 years ago by beausoleilmo

I extracted some information of a VCF file to output summary statistics. I have 96 individuals that were sequenced from a RADseq protocol. I ran …
Can Scaffolds overlap?
written 9.6 years ago by beausoleilmo

I was wondering if Scaffolds could overlap and why. Here is a table showing my scaffolds (Under CHR). The minPos is the minimum position of …
Association mapping with plink: what type of phenotypic data?
written 9.6 years ago by beausoleilmo

I want to know what kind of phenotypic data I can put in PLINK. I have this phenotypic file (pheno1.txt): ![enter image description here][1] I …
Filtering Imputation output: how to filter on a VCF INFO field
written 10.0 years ago by mmagoo6 •

I used the [Michigan Imputation Server][1] to impute data that I have, and got three files per chromosome as output: .dose.vcf.gz, .dose.vcf.gz.tbi, and .info.gz. I …
Using VCFtools to find which population codes individuals in a vcf file belong to
written 10.1 years ago by severalorks

EDIT: The population code for each individual is listed elsewhere and is found on the 1000Genomes site, so the question's essentially been answered I've been …
Converting Vcftools output to R readable format
written 11.2 years ago by pifferdavide

I have used the following Vcftools code to calculate Fst distances from 1000 Genomes Chr 1 data. c:/path --vcf ALL.chr1.phase3_shapeit2_mvncall_integrated_v5.20130502.genotypes.vcf --weir-fst-pop POP1.txt --weir-fst-pop POP2.txt --out …
Filtering Vcf File For Fixed Substitutions
written 13.6 years ago by Rubal7

<p>Hello Everyone,</p> <p>I have 2 VCF files for 2 separate populations and I would like to look for candidate functional sites that are fixed in …