why I have a 1/1 genotype on chromosome X for males
AnsweredHi GTAK team:
I just found the male has the 1/1 genotype on the chromosome X. This should not happen theoretically.
a) GATK version used: 4.1.8
b) Exact command used: standard WDL scripts provided by GATK team
If not an error, choose a category for your question(REQUIRED):
b) What does it mean? there are some error on the haplotypecaller step. But I didn't receive any error info related this.
Much thanks
-
Unless you specify ploidy for chromosome X you will get hemizygous segments as homozygous alt. Though PAR's will give you both homozygous and heterozygous calls.
-
Hey,
Thanks for you answer. But I didn't find the argument to specify ploidy for chromosome X, and only find the --sample-ploidy argument to specify the number of the chromosome.
Could you please send me a link fot that.
Much thanks
-
You need to call X chromosome Non-PAR region seperately using -L argument and by defining --sample-ploidy 1 in the command line GATK will call Non-PAR variants with monoploidy active. There is no other way of defining it. Make sure that you do not call PAR with --sample-ploidy 1 otherwise you will be missing some important variants.
-
Thanks a lot.
-
Thank you for providing your input SkyWarrior!
-
GATK should issue a warning when allosomes are detected without using appropriate ploidy options
Post is closed for comments.
6 comments