[ensembl-dev] VEP text output sort order

cj5 at sanger.ac.uk cj5 at sanger.ac.uk
Mon Oct 17 16:18:39 BST 2011


Hi,
I've noticed that the VEP 2.2. text output order is occasionally not in
genomic order for a small percentage of records. I'm using the hs cache
for this like so:

variant_effect_predictor.pl -v -i <my 120K record vcf> --sift b --polyphen
b --condel b --gene --hgnc --format vcf -o <my output file>
--force_overwrite --cache --dir <my vep_cache>

Output:
...
6_4998993_C/T   6:4998993   T   ENSG00000124787 ENST00000380051 Transcript
 SYNONYMOUS_CODING   561 516 172 K   aaG/aaA -   HGNC=RPP40
10_93502_C/G    10:93502    G   ENSG00000173876 ENST00000482075 Transcript
 DOWNSTREAM
...
[several hundred chr10 records then jump back to chr 6]

Is this a know issue?

Thanks
Chris






More information about the Dev mailing list