Clinical-Genomics / genmod

Annotate models of genetic inheritance patterns in variant files (vcf files)
http://moonso.github.io/genmod/
MIT License
74 stars 18 forks source link

GenMod error: 'list' object has no attribute 'split' #91

Closed fgvieira closed 5 years ago

fgvieira commented 5 years ago

While running genmod (v3.7.3) I got an error:

$ zcat ../joint_variant_calls_annot.vcf.gz | genmod models -f family.ped --vep - 
[2019-07-18 11:47:46,398] WARNING : genmod.commands.annotate_models: 'list' object has no attribute 'split'
Aborted!

Any idea of what it might be? It seems to be on this line:

1 69270 . A G 1061.95 . AC=8;AF=1;AN=8;DP=32;ExcessHet=3.0103;FS=0;MLEAC=8;MLEAF=1;MQ=30.54;QD=33.19;SOR=6.994;CSQ=G|synonymous_variant|LOW|OR4F5|ENSG00000186092|Transcript|ENST00000335137|protein_coding|1/1||ENST00000335137.3:c.180A>G|ENSP00000334393.3:p.Ser60=|180|180|60|S|tcA/tcG|rs201219564||1||SNV|HGNC|14825|YES|||CCDS30547.1|ENSP00000334393|Q8NH21||UPI0000041BC1|||||hmmpanther:PTHR26451&hmmpanther:PTHR26451:SF72&Gene3D:1.20.1070.10&Superfamily_domains:SSF81321&Low_complexity_(Seg):seg&Transmembrane_helices:TMhelix&Prints_domain:PR00237&PROSITE_profiles:PS50262|||||||||||0.838|0.3591|0.7932|0.8482|0.9984|0.8821|0.9146|0.868|0.9005|0.9984|gnomAD_EAS||||||||||rs201219564|6.72432e-01||||||||,G|regulatory_region_variant|MODIFIER|||RegulatoryFeature|ENSR00000918288|TF_binding_site||||||||||rs201219564||||SNV||||||||||||||||||||||||||0.838|0.3591|0.7932|0.8482|0.9984|0.8821|0.9146|0.868|0.9005|0.9984|gnomAD_EAS||||||||||rs201219564|6.72432e-01|||||||| GT:AD:DP:GQ:PL 1/1:0,10:10:30:337,30,0 1/1:0,8:8:24:268,24,0 1/1:0,5:5:15:149,15,0 1/1:0,9:9:27:316,27,0

PS - it seems that this error also affects genmod compound

moonso commented 5 years ago

There is an annotation in the VEP field that breaks the VCF spec.: |ENSP00000334393.3:p.Ser60=|