When to remove releated samples in GATK
Hello Guys,
I am generating my VCF files with the GATK pipeline including joint calling GVCFs. The joint calling process calculates cross-sample information like QD, MQ... based on all samples in the GVCF. However afterwards when receiving the VCF files I (might) remove samples according to ancestry or relatedness. Thus the filters might be partially based on samples not present anymore in the data.
How do you solve this? Do you just apply 2 pipelines? One to remove samples and one to process the remaining ones?
Best,
Andreas
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Missing ML/PL FORMAT fields in GVCFtyper joint calling output (Ultima Genomics + DNAscope)
written by SeoG •
-
DRAGEN gVCF files and joint calling GATK-GenomicsDB
written by brunomiwa •
-
Liftover using gVCF
written by williamsbrian5064
-
Preparing data for kinship calculations
written by Filago •
-
> 90% samples, DP > 10
written by Filago •
-
Considerations when merging a gVCF with a VCF?
written by jon.klonowski
-
Strange Allele Frequency in vcf file
written by GP_PG •
-
Platypus based joint variant calling
written by agoel •
-
Filtering out the variants annotated as "FILTER" from vcf file using GATK
written by lakhujanivijay
-
QD value from GATK in joint VCF
written by deepmehra.atgc •
