GATK4 -known-sites input Recalibrate Base Quality Scores
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?
• 2,401 views
•
link
updated
by
Ram
•
written
by
missTique •
•
written
by
missTique •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
GATK BQSR error — Reference and BAM file chromosome name mismatch (“chr” vs. no “chr”)
written by iamsmor •
-
RNA variant calling - SNP and indels files
written by littlebioinformatician •
-
BaseRecalibrator User error.
written by ChumBucket2024 •
-
Apply BSQR for Targeted Sequencing
written by Manuel Sokolov Ravasqueira •
-
Problem with sequence dictionary detection in BaseRecalibrator
written by Arda •
-
Why GATK4 doesn't carry out local realignment around indels ?
written by ManuelDB
-
Reference file for base recalibration
written by priya.bmg •
-
GATK BaseRecalibrator Step
written by Seigfried
-
Can I estimate LD among species using SNP data?
written by tianying •
-
Creating a "reference sequence" using velvet?
written by jt.bricker •
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.