Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

Wrong genotype correction for hemizygote after PhaseByTransmission

Hello all,

I have a WEX data set for a trio and I followed the best practices for variant calling. PhaseByTransmission updated some genotypes, including one we already validated by sanger.

See attachments for the variants and the PED. First variant in the file is after phasing. The affected male is in the last column, his parents are the preceding 2. The variant gets updated from 1/1 to 1|0 but since it is chrX and the affected individual is male I would expect this never to happen.

The command I used:

java -Xmx4g -jar GenomeAnalysisTK-2.8-1-g932cd3a/GenomeAnalysisTK.jar -R hg19.fa -T PhaseByTransmission \
-V UG8.filtered.SNP.InDel.vcf \
-ped /home/corona/Desktop/pedfiles/VUSerie1.ped \
-o UG8.phased.filtered.SNP.InDel.vcf \
-pedValidationType SILENT

I had to use -pedValidationType SILENT because I had a family in my VCF that was not in my ped. I extracted the samples from the vcf and tried again with STRICT validation and got the same results.

Cheers,
Martin

Best Answer

Answers

  • ZaagZaag Member ✭✭
Sign In or Register to comment.