Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • andy11
    Junior Member
    • Dec 2010
    • 5

    #1

    Calling heterozygosity in NGS data

    In genotyping from NGS data, what is the best way to account for bias (sequencing, mapping, etc) when identifying a locus as homozygous vs heterozygous?

    Thanks for any help!
  • Heisman
    Senior Member
    • Dec 2010
    • 534

    #2
    Not quite sure what you mean by locus. Do you mean specific SNP calls? If so there are many programs out there to call SNPs and will attempt to correctly determine if they are homozygous or heterozygous.

    As for bias, it's good to make sure to remove any PCR duplicates.

    Comment

    • andy11
      Junior Member
      • Dec 2010
      • 5

      #3
      Yes, exactly. I was wondering what those programs do to determine if a locus is homozygous or heterozygous.

      Comment

      • Heisman
        Senior Member
        • Dec 2010
        • 534

        #4
        SAMtools mpileup is pretty good. It's the only common one I'm really familiar with so I don't want to try to compare it to others in any way, but I know it's used fairly commonly.

        Comment

        • sdvie
          Member
          • Jul 2010
          • 68

          #5
          Originally posted by andy11 View Post
          Yes, exactly. I was wondering what those programs do to determine if a locus is homozygous or heterozygous.
          The programs that I know (samtools mpileup + bcf tools and GATK unified genotyper), basically count the number of reads in your experiment supporting one or another allele at a SNP position, and then they estimate whether the ratio of 1st allele to 2nd allele is 50 : 50 or 100 : 0 (usually one of the alleles corresponds to the reference sequence, and the other to the alternative sequence).

          The field "AF" (allele frequency) in the vcf-format file gives you this information.

          Hope that helps.

          Cheers,
          Sophia

          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, Yesterday, 11:13 AM
          0 responses
          16 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, Yesterday, 10:48 AM
          0 responses
          17 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-15-2026, 10:52 AM
          0 responses
          35 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...