Is there a command line option for a single region (not a region file)?
Does anyone have suggestions for rapidly running GenotypeGVCF? The last step (merging) took over three weeks for ~ 150 human genomes.
Thanks!
Is there a command line option for a single region (not a region file)?
Does anyone have suggestions for rapidly running GenotypeGVCF? The last step (merging) took over three weeks for ~ 150 human genomes.
Thanks!
--intervals / -L
One or more genomic intervals over which to operate
Use this option to perform the analysis over only part of the genome. This argument can be specified multiple times. 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).
Engine parameters should work with every tool.
After genotyping the gVCFs over small intervals, you can use bcftools concat to produce a single VCF for all individuals over all regions.
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.