Variant annotation tutorial: Difference between revisions
Jump to navigation
Jump to search
No edit summary |
No edit summary |
||
| Line 2: | Line 2: | ||
== Slicing and dicing of VCF files == | == Slicing and dicing of VCF files == | ||
tabix -h allbt.vcf.gz 18 >BT18.vcf | |||
tabix -p vcf BT18.vcf.gz | |||
Revision as of 12:46, 6 March 2015
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