I am using the GATK UnifiedGenotyper to generate a gVCF file based using the Emit_All_Sites option.
The problem I have is that if I have a single chromosome BAM such as BAM for only say chromosome 3 and I want gVCF data for only chromosome 3, is it possible to do that? So I want a gVCF file with variant and non variant sites only for chromosome 3
Right now when I run GATK then it starts processing all the chromosomes. Is there a way to only generate gVCF data for one chromosome even if genome truly means all the chromosomes?
1 answer
You want the "intervals" argument 'L', as described in the docs.
From the detailed description, linked from the UnifiedGenotyper page
You can use samtools-style intervals either explicitly on the command line (e.g. -L chr1 or -L chr1:100-200) or by loading in a file containing a list of intervals (e.g. -L myFile.intervals).
Log in to answer this question.
-
CombineGVCFs with GATK4.0
written by Yong •
-
How to convert gvcf file to normal vcf file ?
written by DareDevil
-
How to force GATK haplotypecaller to output all genomic position (even non variant positions)
written by Nicolas Rosewick
-
GATK UnifiedGenotyper Stuck at "Preparing for traversal over 1 BAM files|"
written by sevvalakturk96 •
-
How to combine 2000 samples in multi-sample SNP calling
written by yuzhe891 •
-
GATK joint genotyping
written by win
-
VariantRecalibration step issue
written by win
-
Block break from gVCF files
written by win
-
using GATK UnifiedGenotyper on a single chromosome
written by stefano.iantorno •
-
Query Gvcf Against My Variants
written by win
