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

Issues with group-expr in Mutect2 VCF #142

Open
fridells51 opened this issue Aug 22, 2022 · 0 comments
Open

Issues with group-expr in Mutect2 VCF #142

fridells51 opened this issue Aug 22, 2022 · 0 comments

Comments

@fridells51
Copy link

I am having difficulty getting slivar expressions to show up when using a VCF generated by GATK Mutect2. The command I am using:

./slivar expr --vcf git_sample.txt \
      --alias test.txt \
      -o test.vcf \
      --group-expr 'tumor_normal:(tumors.het && normals.hom_ref)' \
      --group-expr 'tumor_het:(tumors.alts == 1)'

The only time I have ever gotten an expression to work is when i do something like "normals.alts != 0" or "tumors.alts != 0"
The variants were decomposed at multiallelic sites using bcftools. Here is the console output:
Screen Shot 2022-08-22 at 4 45 11 PM

The group file and sample vcf are attached:

test.txt

git_sample.txt

Thank you for your time.

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

1 participant