Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
Somatic mutations of epithelial cells (HaplotypeCaller vs Mutect)
Hello, I'd like to extract somatic mutations from normal epithelial cells of cancer patients.
And it's confusing which tool to use. HaplotypeCaller or Mutect?
I know Mutect is for calling somatic mutations but the data is not about cancer cells.
Mutect is different from HaplotypeCaller because it is not based on diploid calculation.
Cancer cells are usually not diploid, so Mutect is useful in that case. (Is this correct?)
But I'm dealing with epithelial cells, so if reads don't show as diploid, then it's presumably sequencing error.
In the case above, it seems like triploid, so Mutect will detect the base pair as somatic variant while HaplotypeCaller might not call it as a variant.
Again, I'm using normal epithelial cells, not cancer cells. So it should be diploid.
Should I use HaplotypeCaller or Mutect?
Do you also think that removing all the variants which already exist in dbsnp database from results of HaplotypeCaller is a correct process for getting somatic mutations from normal epithelial cells?