[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