Skip to content

VEP on non-toplevel coords

Marek Szuba requested to merge github/fork/willmclaren/master into master

Created by: willmclaren

Adding in some code and supporting stuff to allow the VEP endpoints to deal with input on non-toplevel coords; specifically this is targetted at input on alt haplotypes, e.g.

/vep/human/region/HSCHR6_MHC_MANN_CTG1:3848158-3848158:1/G?content-type=application/json

On previous versions this just returns an intergenic result, which is misleading as a) there's no error generated since the slice fetch test passes and b) if you transform to toplevel you do get consequences:

[ { assembly_name: "GRCh38", end: 32582565, seq_region_name: "CHR_HSCHR6_MHC_MANN_CTG1", strand: 1, transcript_consequences: [ { variant_allele: "G", hgnc_id: "HGNC:4948", strand: -1, gene_symbol: "HLA-DRB1", transcript_id: "ENST00000411959", gene_id: "ENSG00000227826", biotype: "protein_coding", gene_symbol_source: "Uniprot_gn", consequence_terms: [ "intron_variant" ], impact: "MODIFIER" }, etc etc

There is a small issue in that the returned coordinates are then different to the input, but there's no real way around this.

Merge request reports