Does anyone have any experience with SNP calling using GATK and exome sequencing data (illumina paired end)? We tried to run it through normally, like it was whole genome sequence and we did not get any SNPs called (the reads were aligned using bwa). Do we need to do anything with exome positions vs just the whole genome using an option like -L? Any help would be appreciated.
-
-
I get plenty of SNPs in an Illumina sequenced, bwa aligned exome with GATK, and I don't use intervals to restrict calling to the bait/target positions.
Care to share your workflow? Hard to diagnose otherwise.Last edited by Bukowski; 07-10-2012, 11:46 AM. -
We did alignment using bwa run as default. We then used samtools to convert to bam and sort. Then we ran the file through picard AddOrReplaceReadGroups to fix the read groups so that GATK would like it. Then we ran the file through GATK using the UnifiedGenotyper using default options. Our GATK isn't the newest version. Is there an issue with older versions of GATK?Comment
-
There are many issues with old versions of GATK, but having used most of them 'not calling genotypes' isn't one of them. Which version? Can you please post your UnifiedGenotyper command line?Originally posted by ercfrtz View Post
To reiterate it is impossible to diagnose issues where you do not post what you're doing.
The first thing to do though is to update GATK, if you attempt to go for support on https://getsatisfaction.com/gsa it is the first thing they will tell you to do.Comment
-
The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.
This convergence of genetics, immunology, and computation...09-01-2026, 05:41 AM
Comment