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

Empty ALT genotypes? #81

Open
rajwanir opened this issue Apr 16, 2024 · 1 comment
Open

Empty ALT genotypes? #81

rajwanir opened this issue Apr 16, 2024 · 1 comment

Comments

@rajwanir
Copy link

Hello @jzieve

I see some of the genotypes have empty ALT allele. This could be simply missing calls. But instead of using . or N, the field is simply empty. In other records I do see the ALT being populated with an N. I can fix it by adding . or N in place of empty ALT field. But wanted to check with you if there is any reason that it might be intentionally set as empty by the GTCtoVCF.

Example records:

chr1 47851 cnvi0146654 C . PASS . GT:GQ ./.:0
chr1 50251 cnvi0146656 T . PASS . GT:GQ ./.:0
chr1 51938 cnvi0151530 T N,A . PASS . GT:GQ ./.:0
chr1 52651 cnvi0146655 T . PASS . GT:GQ ./.:0
chr1 55338 cnvi0159124 A N . PASS . GT:GQ ./.:0
chr1 64251 cnvi0146663 A . PASS . GT:GQ ./.:0
chr1 65451 cnvi0147451 A . PASS . GT:GQ ./.:0
chr1 80386 rs3878915 C A . PASS . GT:GQ ./.:0
chr1 82154 rs4477212 A T,C . PASS . GT:GQ 1/1:3
chr1 82620 cnvi0052563 A N . PASS . GT:GQ ./.:0

Thanks in advance for your help.

@jzieve
Copy link
Collaborator

jzieve commented Apr 16, 2024

I don't know of any reason why it should be empty off the top of my head. What product/manifest is this? Some of these seem to be CNV related targets and maybe could be filtered out of the VCF via the --filter-loci option. Ref/Alt is based off the manifest fields such as SNP, SourceSeq, RefStrand and the bases when looked up on the reference genome. So something must be different about these loci.

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