how to use filter_vep to filter vcf file
0
2
Entering edit mode
12 weeks ago
htchd ▴ 20

I'm trying to filter vcf files with the following criteria: 1) Variants with alternate allele depth (i.e. supporting reads) > 5
2) Variants that are deletion > 25 BPs which are likely introns confused as deletion

filter_vep \  
  -i test_SVs.vcf.gz \  
  -o test_panel.vcf.gz \  
  --format vcf \  
  --force_overwrite \  
  --only_matched \  
  --filter "(FORMAT/AD[0:1] > 5) and ((VARIANT_TYPE is not deletion) or (VARIANT_TYPE is deletion and LENGTH <= 25))"

is this the correct way to do it?

vcf filter • 468 views
ADD COMMENT
0
Entering edit mode

are these VCFs derived from transcripts? why would an intron be missing from a genomic VCF?

ADD REPLY
0
Entering edit mode

yes, they are from RNAseq

ADD REPLY

Login before adding your answer.

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