Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

ASK: How can I filter controversal imputed result? #248

Open
duaghk opened this issue Dec 18, 2024 · 1 comment
Open

ASK: How can I filter controversal imputed result? #248

duaghk opened this issue Dec 18, 2024 · 1 comment

Comments

@duaghk
Copy link

duaghk commented Dec 18, 2024

Hello,

I am currently using GLIMPSE2 for rare variant detection and came across a confusing result.

image

In the attached image, you can see that there are multiple imputed results (3 rows) for the same genomic position (chr1:10397). Each result appears to differ in the REF and ALT alleles, while the genotype calls for all three are 0|0.

Here are my questions:

How should I interpret these results? Why are there multiple imputation outcomes for the same position?
How should I filter or prioritize these results? For downstream analysis, which entry should I use or trust in this case?
If you have any suggestions on how to handle or clean such data, I’d greatly appreciate it.

I used GLIMPSE v2.0.0, and input vcf was generated by: LP-WGS -> DeepVariant call -> Filter PASS & DP >=5

Thank you in advance!

@srubinacci
Copy link
Collaborator

Hi,

Your reference panel contain a complex indel (STR maybe?) and the position is represented as 4 different bi-allelic indels and a SNP. This is one of the cases where you typically want to filter out those INDEL variants from the reference panel, as very likely quite unreliable at the position (and probably not well phased too, for the same reason).

Nothing wrong with GLIMPSE here. If unsure, always restrain your reference panel to INDELS to avoid difficult situations like this.
Simone

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants