This shows you the differences between two versions of the page.
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 / | -r / | ||
- | -g / | + | -g / |
-p / | -p / | ||
-c / | -c / | ||
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 ' | + | zgrep ' |
### | ### |