How to output all records in INFO fild ?

Bioconductor / VariantAnnotation

Annotation of Genetic Variants

https://bioconductor.org/packages/VariantAnnotation

23 stars 20 forks source link

How to output all records in INFO fild ? #23

Closed qindan2008 closed 5 years ago

qindan2008 commented 5 years ago

Hello, I have 64539 records in a vcf file. The annotation information is in [INFO] fild, and [INFO] isn't the same in each record. For example, wechatimg4 Not all the records have the "ExAC_Het_AMR" annotation. How to output @info@listData$ExAC_Het_AMR@unlistData the same length as number of vcd records? And the missing data represented by "NA" or "." ?

Thank you very much !

lawremi commented 5 years ago

You could call drop(info(vcf)$ExAC_Het_AMR). Btw, please avoid accessing the internal slots of the object. It's more trouble than it's worth. Also, you might want tot ask future questions like this on support.bioconductor.org.

qindan2008 commented 5 years ago

OK, Thanks.

vobencha commented 5 years ago

Looks like this one is resolved. Closing the issue.