Changes between Version 30 and Version 31 of SOPs/variant_calling


Ignore:
Timestamp:
03/21/19 14:56:10 (6 years ago)
Author:
thiruvil
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • SOPs/variant_calling

    v30 v31  
    109109== Annotate effect(s) of variants on genes ==
    110110
     111  * snpEff
     112
    111113Use snpEff (assuming snpEff has gene+protein annotations for your genome):
    112114{{{
     
    118120snpEff -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
    119121}}}
     122
     123
     124  * Ensembl's Variant Effect Predictor (VEP)
     125For 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
     128vep -i chr1_subset.vcf --database
     129}}}
     130
    120131
    121132See [wiki:SOPs/vcf Interpreting VCF files ] for more information.