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
            16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
            by SEQadmin2





            NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
            10-30-2026, 06:46 AM
          • 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

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, Today, 07:59 AM
          0 responses
          13 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 11-05-2026, 09:25 AM
          0 responses
          32 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 11-03-2026, 11:54 AM
          0 responses
          45 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-28-2026, 09:57 AM
          0 responses
          169 views
          0 reactions
          Last Post SEQadmin2  
          Working...