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.

UnifiedGenotyper failure

I have twice run UnifiedGenotyper and the resultant .vcf file contains only part of chromosome 20. I do not see what I am doing wrong. Neither do the other two people in the lab who have extensive experience with GATKcat

Best Answer

Answers

  • Nick_BeesonNick_Beeson Member

    Human error. I had a -L 20 argument in the command that I did not want, but somehow got in (cut and paste?). That is so human and makes me feel so stupid. But it is just the way we are built. We humans see what we expect.. This is exactly why GATK exists, it never overlooks anything it is set to detect. So sorry to bother you and thanks for the response.

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie admin

    Hah, don't beat yourself up, it happens all the time :)
    I'm glad your problem is resolved.

Sign In or Register to comment.