GATK4 -known-sites input Recalibrate Base Quality Scores
0
0
Entering edit mode
3.1 years ago
missTique • 0

I am working with a non-reference plant species that I want to call variants from after aligning to a closely related species reference genome. I am following the GATK4 best practices pipeline, but I would like to know how I should proceed with the BaseRecalibrator step, given that it requires "-known-sites" input.Since the data I am working with is non-reference, I don't have any known sites. What could I use as an input for this step? Or am I forced to move on without it?

Plant GATK NGS Variant-Calling • 1.7k views
ADD COMMENT
0
Entering edit mode

There is some guidance on this on the official site. https://gatk.broadinstitute.org/hc/en-us/articles/360035890531-Base-Quality-Score-Recalibration-BQSR-

Essentially you will have to create your own database of variants, and then use this as your --known-sites.

ADD REPLY

Login before adding your answer.

Traffic: 2008 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6