# HG changeset patch
# User devteam
# Date 1426789671 14400
# Node ID 34c3b846e0b6a319c5b73bedb697cfce8ca9c02f
Uploaded
diff -r 000000000000 -r 34c3b846e0b6 macros.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml Thu Mar 19 14:27:51 2015 -0400
@@ -0,0 +1,28 @@
+
+
+
+ vcflib
+
+
+
+
+
+
+
+
+
+
+
+@misc{Garrison2015,
+ author = {Garrison, Erik},
+ year = {2015},
+ title = {vcflib},
+ publisher = {GitHub},
+ journal = {GitHub repository},
+ url = {https://github.com/ekg/vcflib},
+}
+
+
+
+ is a part of VCFlib toolkit developed by Erik Garrison (https://github.com/ekg/vcflib).
+
diff -r 000000000000 -r 34c3b846e0b6 test-data/vcffilter-test1.vcf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/vcffilter-test1.vcf Thu Mar 19 14:27:51 2015 -0400
@@ -0,0 +1,26 @@
+##fileformat=VCFv4.0
+##fileDate=20090805
+##source=myImputationProgramV3.1
+##reference=1000GenomesPilot-NCBI36
+##phasing=partial
+##filter="DP > 10"
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##FILTER=
+##FILTER=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##ALT=
+##ALT=
+#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA00001 NA00002 NA00003
+20 14370 rs6054257 G A 29 PASS AF=0.5;DP=14;NS=3;DB;H2 GT:GQ:DP:HQ 0|0:48:1:51,51 1|0:48:8:51,51 1/1:43:5:.,.
+20 17330 . T A 3 q10 AF=0.017;DP=11;NS=3 GT:GQ:DP:HQ 0|0:49:3:58,50 0|1:3:5:65,3 0/0:41:3:.,.
+20 1230237 . T . 47 PASS AA=T;DP=13;NS=3 GT:GQ:DP:HQ 0|0:54:.:56,60 0|0:48:4:51,51 0/0:61:2:.,.
diff -r 000000000000 -r 34c3b846e0b6 test-data/vcflib.vcf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/vcflib.vcf Thu Mar 19 14:27:51 2015 -0400
@@ -0,0 +1,31 @@
+##fileformat=VCFv4.0
+##fileDate=20090805
+##source=myImputationProgramV3.1
+##reference=1000GenomesPilot-NCBI36
+##phasing=partial
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##FILTER=
+##FILTER=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##ALT=
+##ALT=
+#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA00001 NA00002 NA00003
+19 111 . A C 9.6 . . GT:HQ 0|0:10,10 0|0:10,10 0/1:3,3
+19 112 . A G 10 . . GT:HQ 0|0:10,10 0|0:10,10 0/1:3,3
+20 14370 rs6054257 G A 29 PASS NS=3;DP=14;AF=0.5;DB;H2 GT:GQ:DP:HQ 0|0:48:1:51,51 1|0:48:8:51,51 1/1:43:5:.,.
+20 17330 . T A 3 q10 NS=3;DP=11;AF=0.017 GT:GQ:DP:HQ 0|0:49:3:58,50 0|1:3:5:65,3 0/0:41:3:.,.
+20 1110696 rs6040355 A G,T 67 PASS NS=2;DP=10;AF=0.333,0.667;AA=T;DB GT:GQ:DP:HQ 1|2:21:6:23,27 2|1:2:0:18,2 2/2:35:4:.,.
+20 1230237 . T . 47 PASS NS=3;DP=13;AA=T GT:GQ:DP:HQ 0|0:54:.:56,60 0|0:48:4:51,51 0/0:61:2:.,.
+20 1234567 microsat1 G GA,GAC 50 PASS NS=3;DP=9;AA=G;AN=6;AC=3,1 GT:GQ:DP 0/1:.:4 0/2:17:2 1/1:40:3
+20 1235237 . T . . . . GT 0/0 0|0 ./.
+X 10 rsTest AC A,ATG 10 PASS . GT 0 0/1 0|2
diff -r 000000000000 -r 34c3b846e0b6 tool_dependencies.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml Thu Mar 19 14:27:51 2015 -0400
@@ -0,0 +1,9 @@
+
+
+
+
+
+
+
+
+
diff -r 000000000000 -r 34c3b846e0b6 vcffilter.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/vcffilter.xml Thu Mar 19 14:27:51 2015 -0400
@@ -0,0 +1,78 @@
+
+ filter VCF data in a variety of attributes
+
+ macros.xml
+
+
+ tabix
+
+
+
+
+ ln -s "${input1}" input1.vcf &&
+ bgzip input1.vcf &&
+ tabix -p vcf input1.vcf.gz &&
+ vcffilter ${filterList} input1.vcf.gz > "${out_file1}"
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+You can specify the following options within the **Specify filtering expression** box in any combination::
+
+ -f, --info-filter specifies a filter to apply to the info fields of records, removes alleles which do not pass the filter
+ -g, --genotype-filter specifies a filter to apply to the genotype fields of records
+ -s, --filter-sites filter entire records, not just alleles
+ -t, --tag-pass tag vcf records as positively filtered with this tag, print all records
+ -F, --tag-fail tag vcf records as negatively filtered with this tag, print all records
+ -A, --append-filter append the existing filter tag, don't just replace it
+ -a, --allele-tag apply -t on a per-allele basis. adds or sets the corresponding INFO field tag
+ -v, --invert inverts the filter, e.g. grep -v
+ -o, --or use logical OR instead of AND to combine filters
+ -r, --region specify a region on which to target the filtering (must be used in conjunction with -f or -g)
+
+Filters are specified in the form {ID} {operator} {value}::
+
+ -f "DP > 10" # for info fields
+ -g "GT = 1|1" # for genotype fields
+ -f "CpG" # for 'flag' fields
+
+Any number of filters may be specified. They are combined via logical AND unless the --or option is specified. For convenience, you can specify "QUAL" to refer to the quality of the site, even though it does not appear in the INFO fields.
+
+Operators can be any of: =, !, <, >, pipe, &
+
+
+To restrict output to a specific location use the -r option (must be used in conjunction with -g or -f)::
+
+ -r chr20:14000-15000 # only output calls between positions 14,000 and 15,000 on chromosome 20
+ -r chrX # only output call on chromosome X
+
+-----
+
+Vcffilter @IS_PART_OF_VCFLIB@
+
+
+