Question: How do I determine the correct filtering threshold in VCFfilter after running FreeBayes
0
gravatar for jamieburke1020
24 months ago by
jamieburke10200 wrote:

I would like to perform a VCFfilter on my FreeBayes data but I'm not sure how to choose the correct parameters for filtering. I visualized my variants with vcf.iobio. Where do I find the correct number from my data for the following fields: -f "DP > 10" # for info fields -g "GT = 1|1" # for genotype fields -f "CpG" # for 'flag' fields

snp vcffiter • 2.0k views
ADD COMMENTlink modified 23 months ago by chen.randy110 • written 24 months ago by jamieburke10200

Hello - I am not quite sure what the question is. Have you tried these filters and examined the results? Or are you unsure about how to enter them on the form? Thanks for the extra info, Jen, Galaxy team

ADD REPLYlink written 23 months ago by Jennifer Hillman Jackson25k

My question should be: How do I use the VCFfilter to only include sites where the chance of a false positive call is 1 in 10,000 or better according to the VCF qual field? This is how the question was stated in one of steps for a project within Coursera.

ADD REPLYlink written 23 months ago by jamieburke10200
1
gravatar for bon1068
23 months ago by
bon106830
bon106830 wrote:

-f "QUAL > 40" (that is 1 out of 10,000)

ADD COMMENTlink written 23 months ago by bon106830

Correct! The full example with a workflow and history is on this page for review: https://usegalaxy.org/u/coursera/p/genomic-data-science-with-galaxyidentify-polymorphic-sites

ADD REPLYlink written 23 months ago by Jennifer Hillman Jackson25k
1
gravatar for chen.randy
23 months ago by
chen.randy110
chen.randy110 wrote:

“-f ‘QUAL > 40’” parameter will need to be used. A phred scale score of 40 equates to a 1 in 10,000 probability of a false positive call.

ADD COMMENTlink written 23 months ago by chen.randy110
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 172 users visited in the last hour