cortes-ciriano-lab / SComatic

A tool for detecting somatic variants in single cell data
Other
145 stars 20 forks source link

How does SingleCellGenotype.py work and what does it mean if substantial number of cell have different expected and observed cell types? #52

Open prakashraaz opened 5 months ago

prakashraaz commented 5 months ago

Thank you for making this wonderful tool. I have used this pipeline on six samples and was surprised to see very few somatic variants after all the filtering. I noticed that in a substantial number of cases, the reason the filtering failed was because the expected cell type and observed cell type were different. Even though you have tried to answer similar questions elsewhere, do you think you could explain how computing the genotypes for each cell at variant sites works? Also, what does it mean for all the cells of a particular cell type have mismatches in the expected cell type and observed cell type columns? Does this point towards there being some issue with the cell type annotation quality?