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.
i.e.,
ref: A
reads: AAACCCGGG
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?