Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • ercfrtz
    Member
    • Aug 2010
    • 23

    #1

    Exome SNP Calling

    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.
  • Bukowski
    Senior Member
    • Jan 2010
    • 388

    #2
    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.

    Comment

    • ercfrtz
      Member
      • Aug 2010
      • 23

      #3
      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

      • Bukowski
        Senior Member
        • Jan 2010
        • 388

        #4
        Originally posted by ercfrtz View Post
        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?
        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?

        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

        Latest Articles

        Collapse

        • SEQadmin2
          Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
          by SEQadmin2



          Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

          Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
          ...
          10-14-2026, 09:55 AM
        • SEQadmin2
          New Genomics Technologies Take Aim at Long-Standing Limits
          by SEQadmin2


          Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

          We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
          ...
          09-28-2026, 10:25 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, 10-22-2026, 11:13 AM
        0 responses
        23 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-22-2026, 10:48 AM
        0 responses
        20 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-15-2026, 10:52 AM
        0 responses
        36 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-14-2026, 05:39 AM
        0 responses
        42 views
        0 reactions
        Last Post SEQadmin2  
        Working...