filter varaints from vcf sample based on their presence in N number of samples ?
1
0
Entering edit mode
3.8 years ago

Is there any command in vcftools / bcftools aur any other tool to separate variants based on their presence in N number of samples. E.g I have 20 samples in a merged vcf file. I want to separate variants present in a minimum of 15 samples.. ?

vcf SNP alignment • 602 views
ADD COMMENT
1
Entering edit mode
3.8 years ago

not tested.

 bcftools view -i 'N_PASS(GT!="mis" && GT!="RR")>=15 ' input.vcf.gz
ADD COMMENT
0
Entering edit mode

I tried this but it didn't work

ADD REPLY

Login before adding your answer.

Traffic: 1825 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6