<div dir="ltr">Hi all,<br><br>I've been trying to get information of overlapping protein domains for one variant using VEP, but it looks as if the REST API returns more domains than the command line tool. Domains here means the output of the command line switch `--domains`, which, as far as I can tell, is the same as `domains=1` with the `GET vep/:species/id/:id` API request.<br><br>For example, for this single variant I'm using for testing, EGFR T790M, with the GET method above `<a href="https://rest.ensembl.org/vep/human/id/rs121434569?domains=1&content-type=application/json`">https://rest.ensembl.org/vep/human/id/rs121434569?domains=1&content-type=application/json`</a> the `domains` list of the `transcript_consequences` object, lists several ENSP_mappings and also information from CDD, Pfam, PROSITE_profiles, and others. I'm more interested in the Pfam information, which in this case corresponds to a protein tyrosine and serine/threonine kinase, PF07714.<br><br>However, when I run this same variant in the command line (using a VCF file with this single variant as input), I can only obtain information from the ENSP_mappings, but all other databases appear to be missing. The command used was the following:<br><br>```<br>vep --domains --dir_cache /opt/bioResources/vep_106/ --fasta /opt/bioResources/vep_106/homo_sapiens_refseq/106_GRCh38/Homo_sapiens.GRCh38.dna.toplevel.fa.gz --input_file T790M.vcf --output_file T790M.vep.json --cache --offline --json --force_overwrite<br>```<br><br>Does anyone know if this is expected, or how to get the same output of the REST API (regarding the list of protein domains) when using the command line tool? Are custom annotations needed for these cases?<br><br>Many thanks,<br>Pedro<br clear="all"><div><br></div></div>