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.
SNP calling in a "case-control" setting
I wish to look for SNPs that are unique to my case samples, and do not appear in my control samples. To be more precise, their enrichment in the case samples should be significant. Something like a fisher's test.
I have many bam files of control samples, and many of cases. How can I run this?
This is bacterial DNA, so as far as we know, there is a single copy of each chromosome.