Variant annotation tutorial

From HPCwiki
Revision as of 12:46, 6 March 2015 by Megen002 (talk | contribs)
Jump to navigation Jump to search


Slicing and dicing of VCF files

 tabix -h allbt.vcf.gz 18 >BT18.vcf
 tabix -p vcf BT18.vcf.gz


Annotating VCF with rs-numbers

tabix -h BT18.vcf.gz 18:100000-101000 | vcf-annotate -a BT_incl_cons.18.vcf.gz -c CHROM,FROM,ID

something

Variant Effect Predictor

 gunzip -c /home/formacion/COMUNES/IAMZ/data/CIHEAM/MULTISAMPLE_VCF/all.fb.vcf.gz | perl /home/formacion/COMUNES/IAMZ/soft/ensembl-tools-release-78/scripts/variant_effect_predictor/variant_effect_predictor.pl --dir /home/formacion/COMUNES/IAMZ/data/CIHEAM/ReferenceGenome/VEP/ --species bos_taurus -o test2.vep --fork 4 --canonical --sift b --coding_only --no_intergenic --offline --force_overwrite --vcf

snpEff

Polyphen

Provean