filtering heterozygous genotypes

I have merged vcf file. I need to keep heterozygous genotypes with heterozygosity < 0.2.

How can I do through bcftools?

I found this command to remove heterozygous genotypes:

bcftools view -e 'GT="het"' merged_annotated_snps.vcf > merged_annotated_polymorphic_snps.vcf

How can I modify above command to retain genotypes with heterozygosity < 0.2 ?

Also my goal is to keep polymorphic alleles.

Thanks a lot

bcftools heterozygous

• 1,244 views • link

written 1.5 years ago by Bushra Bibi

1 answer

How can I modify above command to retain genotypes with heterozygosity < 0.2 ?

you meant Allele Frequency ? if INFO/AF is defined

bcftools view -i 'AF>=0.2'  in.vcf

otherwise, set AF using bcftools +fill-tags

• 0 views • link

written 1.5 years ago by Pierre Lindenbaum

No i mean heterozygosity such as 'GT="het"'

• 0 views • link

written 1.5 years ago by Bushra Bibi

More posts like this

bcftools not filtering snps
written 1.2 years ago by Bushra Bibi

Dear all, I have called variants through freebayes. Now I am filtering SNPs and Indels. Opening SNPs file I found other variants also. I used …
vcf file filtration w.r.t. missingness per snps and unmapped snps
written 1.9 years ago by Bushra Bibi

Hi, I have wheat hexaploid merged vcf file and I have to filter vcf w.r.t missingness per snp 50%, unmapped snps to remove, and heterozygous …
merging tetraploid vcf files through bcftools
written 2.4 years ago by Bushra Bibi

I have 7 tetraploid vcf files (generated through gatk Hapotypecaller) that I tried to merge through bcftools but It gave me following error: **Incorrect number …
how to extract variants from the vcf.gz files linked below?
written 2.5 years ago by Mahad •

I have tried several commands to extract the variant (non-homozygous) genotype from the two vcf.gz files found below, but no variant genotypes are coming up. …
How to get only polymorphic genotypes from combined vcf file filtering out homozygous, heterozygous…
written 2.6 years ago by Bushra Bibi

Hi everyone, I have a combined vcf file of multiple samples. I have to remove homozygous, heterozygous and missing genotypes and keep only polymorphic genotypes. …
Filter on Allele Balance using BCFTools
written 5.0 years ago by sumudu_rangika •

Hi All, I need to filter my variants based on the following criteria. 1) Include SNP sites with **at least one heterozygous with allele balance(AB) …
filtering on F_PASS(GT="het") in bcftools gives contradictory results
written 5.5 years ago by scottbrainard •

I want to remove sites from a VCF file with an overabundance or an underabundance of heterozygous genotypes. I have 8940 markers to start with: …
bcftools to remove non-variants sites
written 5.9 years ago by shiy05 •

Dear all, I am used to using vcftools to remove non-variants sites. The command I used: ``` vcftools --vcf input.vcf --mac 1 --recode --out input.mac …
Bcftools filtering only 10% missing and heterozygous calls
written 6.4 years ago by User000

Hello, I would like to keep only 90% of missinfg and heterogygous genotype on all my samples. According to [link][1] I can filter missing like …
Pileup Haploid Base Calling, Creating Consensus Sequence
written 14.2 years ago by Lee Katz

<p>Hi all, I couldn't find if this has been asked here. I realize that samtools assumes that you are working with a diploid and will …