Filtering no-calls from VCF file
1
0
Entering edit mode
4.6 years ago
SOHAIL ▴ 330

Hi Everyone, i want to remove all sites with more than 5% of no-call (./.) among all the samples in my VCF file, could you please guide me how i can do that variant filtration??

Thanks!

ngs vcf-tools • 2.2k views
ADD COMMENT
1
Entering edit mode
4.6 years ago

using vcffilterjs https://github.com/lindenb/jvarkit/wiki/VCFFilterJS

 java -jar dist/vcffilterjs.jar -e 'var n=0.0;for(var i=0;i< variant.getNSamples();i++ ) { if(variant.getGenotype(i).isNoCall()) n++;} n/ variant.getNSamples() < 0.05'  input.vcf
ADD COMMENT
0
Entering edit mode

Thanks @Pierre.. That worked perfectly fine for me..

However i am curious and interested to use VCFFILTERJS more, Could you please indicate any guideline for writing java scripts expressions particularly for VCFFILTER -e. I have seen some examples but i wanna know about basic functions like: variant.getNSamples() or variant.getGenotype().isNoCall() , etc.

ADD REPLY
1
Entering edit mode

this tools uses a javascript engine (nashorn) that can read the underlying java classes of the htsjdk library like https://samtools.github.io/htsjdk/javadoc/htsjdk/htsjdk/variant/variantcontext/VariantContext.html . It uses a old and well-known java technology called reflection: https://docs.oracle.com/javase/tutorial/reflect/ .

ADD REPLY

Login before adding your answer.

Traffic: 2194 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