How to filter delly variants?
1
How can I filter delly output to only deletions and ones that PASSED the quality check? I want to use the delly output from delly call.
calling
DNA-seq
delly
variants
• 2.1k views
•
link
updated 3.7 years ago by
trausch
★
1.9k
•
written 3.7 years ago by
gt
▴
30
Bcftools works best for filtering VCF/BCF files. For instance
bcftools query -f '%CHROM\t%POS\t%INFO/END\t%ID\t%INFO/SVTYPE\t%QUAL\t%FILTER\n' delly.bcf
generates a table with all the information you want. Alternatively, you can also filter the VCF/BCF directly using:
bcftools view -i 'FILTER="PASS" & INFO/SVTYPE="DEL"' delly.bcf
Login before adding your answer.