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!
Thanks for any help!
-
-
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. -
Yes, exactly. I was wondering what those programs do to determine if a locus is homozygous or heterozygous.Comment
-
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
-
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).Originally posted by andy11 View Post
The field "AF" (allele frequency) in the vcf-format file gives you this information.
Hope that helps.
Cheers,
SophiaComment
-
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
Comment