Changes between Version 30 and Version 31 of SOPs/variant_calling
- Timestamp:
- 03/21/19 14:56:10 (6 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/variant_calling
v30 v31 109 109 == Annotate effect(s) of variants on genes == 110 110 111 * snpEff 112 111 113 Use snpEff (assuming snpEff has gene+protein annotations for your genome): 112 114 {{{ … … 118 120 snpEff -c /usr/local/share/snpEff/snpEff.config -no-downstream -no-intergenic -no-intron -no-upstream -no-utr -s A_snpEff.html SacCer_Apr2011.18 A_reads.bt2.sorted_unique.filtered.vcf > A_reads.bt2.sorted_unique.filtered.snpEff.vcf 119 121 }}} 122 123 124 * Ensembl's Variant Effect Predictor (VEP) 125 For large number of variants, download the annotation and use --cache or --offline instead of --database, see[[http://www.ensembl.org/info/docs/tools/vep/script/vep_options.html | documentation]] 126 {{{ 127 #default human, set by --species 128 vep -i chr1_subset.vcf --database 129 }}} 130 120 131 121 132 See [wiki:SOPs/vcf Interpreting VCF files ] for more information.