How to handle allosomes and chrM in GATK?
I'm using whole-genome sequencing (human) data for variant calling using GATK HaplotypeCaller. By default, the HaplotypeCaller returns diploid calls on all the chromosomes, but for male haploid call should be made for chrY, non-pseudoautosomal chrX and chrM/chrMT for both male and female. So, when dealing with these calls should I use these regions as intervals in HaplotypeCaller and run using ploidy=1?
The 1000 Genome project (Phase 3) method section(supplementary info page 78) they have mentioned about following a different protocol for chrY but the details are not very clear to me.
• 1,887 views
•
link
written
by
Arup Ghosh
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
whatshap ploidy error on female sample
written by eebloom •
-
low sequencing read depth on Y chromosome
written by karen2 •
-
Identifying/filtering pairwise IBD regions in diploid genomes
written by shpak.max •
-
Genotyping low read depth sequences using GATK
written by shpak.max •
-
gender determination and chrX CN calls
written by erikt •
-
bcftools for plant SNPs
written by Ric
-
Diploid calls on chrY?
written by dec986
-
Ploidy in copy number analysis
written by rahul.nahar •
-
Genotype representation for male chromosome X and Y
written by bioinforesearchquestions
-
Getting pileup information from Sam files
written by anjali.gopal91 •
Doesn't answer your question for using HaplotypeCaller., but as an alternative have a look at
bcftools. There you can define the ploidy based on region and sample sex. See the description of--ploidy-fileand--samples-filehere.