Question: How do I determine the correct filtering threshold in VCFfilter after running FreeBayes
0
gravatar for jamieburke1020
24 months ago by
jamieburke1020 • 0 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 COMMENT • link • modified 23 months ago by chen.randy • 110 • written 24 months ago by jamieburke1020 • 0

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 REPLY • link written 23 months ago by Jennifer Hillman Jackson ♦ 25k

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 REPLY • link written 23 months ago by jamieburke1020 • 0
1
gravatar for bon1068
23 months ago by
bon1068 • 30
bon1068 • 30 wrote:

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

ADD COMMENT • link written 23 months ago by bon1068 • 30

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 REPLY • link written 23 months ago by Jennifer Hillman Jackson ♦ 25k
1
gravatar for chen.randy
23 months ago by
chen.randy • 110
chen.randy • 110 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 COMMENT • link written 23 months ago by chen.randy • 110
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