User Tools

Site Tools


homework_7_functional_annotation

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
homework_7_functional_annotation [2015/12/02 15:44]
scott
homework_7_functional_annotation [2015/12/03 14:22] (current)
scott /* Code for Homework */
Line 2: Line 2:
 ==== Notes ==== ==== Notes ====
  
-  * BioFrontiers system admins have not installed ANNO just yet. Should be ready by class time tomorrow (12/3).* 
  
  
Line 21: Line 20:
 # Some modules you'll need # Some modules you'll need
 module load tabix_0.2.6 module load tabix_0.2.6
-module load anno +module load anno_1.0
-module load plink_latest+
  
 # I encourage you to annotate your own genome, but for purposes of this  # I encourage you to annotate your own genome, but for purposes of this 
Line 40: Line 38:
      -o chrALL.filtered.PASS.beagled.HG00096.rsID.anno.vcf.gz \      -o chrALL.filtered.PASS.beagled.HG00096.rsID.anno.vcf.gz \
      -r /Users/scvr9332/reference_data/gotcloud.ref/hs37d5.fa \      -r /Users/scvr9332/reference_data/gotcloud.ref/hs37d5.fa \
-     -g /Users/scvr9332/reference_data/anno/refFlat.gencode.v11.gz \+     -g /Users/scvr9332/reference_data/anno/refFlat_hg19.txt.gz \
      -p /Users/scvr9332/reference_data/anno/priority.txt \      -p /Users/scvr9332/reference_data/anno/priority.txt \
      -c /Users/scvr9332/reference_data/anno/codon.txt      -c /Users/scvr9332/reference_data/anno/codon.txt
Line 66: Line 64:
 ### The following command tells me at how many sites this individual ### The following command tells me at how many sites this individual
 ### has a heterozygous stop gain allele: ### has a heterozygous stop gain allele:
-zgrep '\s0|1:\|\s1|0:' stop_gains.vcf | wc -l+zgrep '\s0|1:\|\s1|0:' stop_gains.vcf.gz | wc -l
  
 ###----------------------------------### ###----------------------------------###
homework_7_functional_annotation.1449096251.txt.gz · Last modified: 2015/12/02 15:44 by scott