Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

alt allels error #8953

Open
gaochenx opened this issue Aug 14, 2024 · 1 comment
Open

alt allels error #8953

gaochenx opened this issue Aug 14, 2024 · 1 comment

Comments

@gaochenx
Copy link

Instructions

gatk --java-options "-Xmx256g -XX:ParallelGCThreads=10" GenotypeGVCFs -R ./Zm_Mo17.fa -stand-call-conf 0 -ploidy 1 -V gendb:///share/home/gaochenxu/01.prj/50.pan-maize/2.sv/pan/03.vcf/gatk10 -O ./gatk10.vcf.gz --cloud-prefetch-buffer 10000 --cloud-index-prefetch-buffer 10000 --genomicsdb-max-alternate-alleles 110 --max-alternate-alleles 100 --gcs-max-retries 1000

Bug Report

Hello,
when I was using GATK GenotypeGVCFs function call variant, I found that the result could be more than 3 alleles in alt. However, none of the samples in the generated vcf had genotypes other than 012

Expected behavior

alt genotype, there should be corresponding sample support

Actual behavior

But in fact, all the samples have only genotype 0/0 1/1 2/2, and no genotype 3/3 or above

So I was wondering, is there a problem with where I set it up

@gokalpcelik
Copy link
Contributor

When ploidy is set to 1 none of your samples should have GTs with 0/0 1/1 2/2. Are you sure that your original VCF calls are made with ploidy 1 as well?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants