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

Next revision
Previous revision
homework_7_functional_annotation [2015/12/02 15:43]
scott Created page with "==== Notes ==== ==== Code for Homework ==== <syntaxhighlight lang="bash"> # PSYC 7102 -- Statistical Genetics # Homework #7: Functional Annotation # Due: December 10th # Ove..."
homework_7_functional_annotation [2015/12/03 14:22] (current)
scott /* Code for Homework */
Line 20: 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 39: 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 65: 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.1449096199.txt.gz · Last modified: 2015/12/02 15:43 by scott