VariantAnnotation as.vranges with asterisk in ALT field
1
0
Entering edit mode
@jeremy-leipzig-4924
Last seen 5.4 years ago

GATK has recently started using asterisks in the ALT field for certain variants (i believe a SNP and deletion at one site)

chr22    47461908    .    TAA    *    859.92    PASS    AC=2;AF=0.333;AN=6;BaseQRankSum=-3.105;ClippingRankSum=0.773;DP=119;FS=12.256;MLEAC=2;MLEAF=0.333;MQ=70;MQRankSum=0.042;QD=8.77;ReadPosRankSum=0.126;SOR=2.17;OLD_MULTIALLELIC=chr22:47461908:TAA/TA/*    GT:AD:DP:GQ:PL    0/1:16,14:30:99:350,0,527    0/.:20,0:39:99:193,254,652    0/1:17,15:32:99:347,0,483

This is causing the VRanges validator to crash

> variants <- as(vcf, "VRanges")

Error in validObject(.Object) : 

  invalid class “VRanges” object: if 'alt' is 'NA', then 'altDepth' should be 'NA'

 

using VariantAnnotation_1.16.4    and  GATK 3.4-46-gbc02625

variantannotation gatk • 1.6k views
ADD COMMENT
2
Entering edit mode
@michael-lawrence-3846
Last seen 2.4 years ago
United States

I removed the assertion in devel, 1.17.22.

ADD COMMENT

Login before adding your answer.

Traffic: 769 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6