Would I want to call variants on hs37d5 for any reason?
Of course, one should use hs37d5.fa for read mapping.
However, should I bother to call variants on the decoy contig hs37d5? GATK Haplotype Caller is spending ages on the contig and I will never use the variant calls.
Would this contig be of any use downstream? Would it make sense to use the GATK variant recalibration on the variant calls of this contig or could they be safely ignored?
Thanks!
1 answer
Hi,
we included hs37d5 for the mapping step to get rid of mismapped reads and to increase the SNV-calling sensitivity. See also GATK forum decoy discussion
For downstream analysis we ignore scaffolds and contigs, since most databases and tools are unaware of these.
Log in to answer this question.
More posts like this
-
GATK - Output differences in Calls when Split by chromosome intervals and without splitting
written by Maverick •
-
Estimating Nucleotide diversity (Pi) from VCF file in polyploid organism (tetraploid)
written by Michael Dondrup
-
Haplotype Caller GATK Query
written by ramshahaya •
-
Query regarding somatic variant calling
written by iti.gupta •
-
What is the different if I call germline SNP/Indel with mutect and HaplotypeCaller?
written by Makplus T
-
Running time in different callers
written by zizigolu
-
GATK; Variant calling RNA-Seq best practices: bam processing steps
written by user230613
-
Multiple RNAseq Library to SNP Call using GATK
written by Bioinfonext
-
Filter/Annotate VCF with dbSNP during variant calling workflow
written by alons
-
Using Gatk On Inbred Species - What To Do With The Heterozygous Calls?
written by William
