Meaning of "%" in expressions for bcftools filter
The current bcftools HOWTO gives several examples of variant filtering, such as:
bcftools filter -i'%QUAL>20' calls.vcf.gzWhat does the "%" mean? Why not simply using -i'QUAL>20' ?
I couldn't find details on this in the EXPRESSIONS section of the manual ("man bcftools").
• 3,818 views
•
link
updated
by
lstbl •
•
written
by
Timothée Flutre •
•
written
by
Timothée Flutre •
1 answer
This is because when you include a % inside of quotes, it will act as a stdin for the program.
See this post: http://serverfault.com/questions/274475/escaping-double-quotes-and-percent-signs-in-cron
Log in to answer this question.
More posts like this
-
Blog: NGS one-liner to call variants
written by barslmn •
-
Weird behavior of `--max-depth` argument from `bcftools mpileup` command when its value is too low
written by sbstevenlee •
-
How To Make Samtools Mpileup Produce Genotype Information At Every Base?
written by komal.rathi
-
Vcf Sort According To Order Of The Reference File
written by Tonyzeng •
-
Why Is The Variant File Produced By Samtools So Large? It Is About ~34G.
written by Jordan
-
Tool: Gemini: Integrative Exploration Of Genetic Variation And Genome Annotations
written by Aaronquinlan
-
Filtering Variant Calls Based On Genotype Quality
written by Kssr
-
Help With Custom Vcf-Annotate Feature In Vcftools
written by 714 •
-
Comparison Of Two Versions Of Gatk Filtering For Exome Sequencing Data.
written by michealsmith
-
How Do I Convert 454 Ace To A Regular Ace?
written by Lee Katz
I'm curious about this as well. I know you have to put it in the single quotes so the shell doesn't think that you are want to to use the modulus function but I tried it both ways and got the same results. I saw earlier that the header format was %CHROM, etc. so maybe its more of a format/best practice kind of thing.