This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Getting frequency of sites fixed within the sample (i.e. divergence sites) from VCF file

Hi there,

I am trying to extract fixed sites within the sample from a VCF file. By fixed sites, I mean those that differ from the reference genome but that are fixed within the sample. 

    REF    ALT ind_1   ind_2  ind_3
1    A      C    1/1    1/1     1/1
2    G      T    1/1    0/1     0/0
3    C      G    1/1    1/1     1/0
4    G      C    0/1    1/1     1/0
5    A      G    1/1    1/1     1/1

For instance, this is was a simplified version of a VCF file. In this case, sites 1 and 5 belong to this category of sites that contribute to divergence. Is there any tool on vcftools or R package that I can use for this purpose?

Thanks in advance.

next-gen snp sequence
• 2,364 views • link
updated by finswimmer finswimmer • written by JGuVa • JGuVa

not clear to me. You want the variants where all the genotypes are homozygous for the ALT allele ?

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Yes, exactly, that is what I need.

• 1 views • link
written by JGuVa • JGuVa

In addition to Pierre: Is your data in this simplified format or a normal vcf?

• 1 views • link
written by finswimmer finswimmer

My file is a normal VCF, I presented it like that just for the sake of the explanation.

• 1 views • link
written by JGuVa • JGuVa

2 answers

using vcffilterjdk http://lindenb.github.io/jvarkit/VcfFilterJdk.html

java -jar dist/vcffilterjdk.jar -e 'return variant.getGenotypes().stream().allMatch(G->G.isHomVar());' in.vcf
• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Using bcftools:

$ bcftools view -i 'COUNT(GT="AA")=N_SAMPLES' input.vcf

or

$ bcftools view -e 'GT[*]!="AA"' input.vcf

fin swimmer

• 1 views • link
written by finswimmer finswimmer

Log in to answer this question.

More posts like this