Entering edit mode
2.5 years ago
addamyeo
•
0
To specify from an earlier post, I need to filter out the passed SNVs, such that only the failed variants remain. The command I used to filter out the failed SNV is:
bcftools isec -f 'PASS,.' -e 'TYPE!="snp"' [file1.vcf.gz] [file2.vcf.gz] -p dir
Thanks!
dup of How to filter out failed SNVs to create new files using bcftools isec.