Haplotype Calling Using GATK HaplotypeCaller to Get Exact 35718025 Positions/Loci
Hi,
I have a VCF file that has 35718025 positions/loci. I am trying to add more accessions (let's say accession1, accession2, and accession3) into this VCF using bcftools merge. When I perform haplotype calling using the GATK HaplotypeCaller, how can I output exact 35718025 positions/loci for the accession1, accession2, and accession3?
I have tried using GATK HaplotypeCaller --force-call-filtered-alleles --alleles <the VCF file that has 35718025 loci>, but the output VCF files have more than 35718025 positions/loci. Does anybody have any suggestions for me to solve this problem?
• 1,425 views
•
link
written
by
yenon118 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
GATK SelectVariants behavior with gvcfs
written by Gianfilippo •
-
HaplotypeCaller Not Detecting Variants
written by akshit •
-
Generate VCF version 4.1 using GATK version 4
written by taniamahmood38 •
-
GATK based variant calling yields excessive Indels
written by evoclive •
-
Calling variants on reads with MAPQ=0 on HaplotypeCaller or bcftools mpileup
written by Anupriya Dalmia •
-
GATK variant calling error
written by endretoth •
-
Bcftools Mpileup bias for alternative alleles?
written by cetin.m
-
How to force GATK haplotypecaller to output all genomic position (even non variant positions)
written by Nicolas Rosewick
-
SNPs not filtered from a GATK vcf file
written by evelyn
-
incorporate indels to gtf/gff file
written by grant.hovhannisyan
what is an "accession" ?
Accession and sample are the same things. The Bcftools use the word sample in their manual. You can use
bcftools query -l <VCF file>to print out all the accessions/samples in your VCF file.https://en.wikipedia.org/wiki/Accession_number_(bioinformatics)