From 487db9abc324453497353b954ee710fb3668404e Mon Sep 17 00:00:00 2001
From: Shifu Chen <chen@haplox.com>
Date: Tue, 16 Feb 2016 22:50:53 +0800
Subject: [PATCH] Add gencode example

---
 README.md | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/README.md b/README.md
index 9991a2e..493db62 100644
--- a/README.md
+++ b/README.md
@@ -94,3 +94,18 @@ while (row = gtf_read_row(stream)) != false
     # do something with row ...
 end
 ```
+
+***locate the gene/exon/intron***
+```julia
+using OpenGene, OpenGene.Reference
+
+# load the gencode dataset, it will download a file from gencode website if it not downloaded before
+# once it's loaded, it will be cached so future loads will be fast
+index = gencode_load("GRCh37")
+
+# locate which gene chr:pos is in
+gencode_locate(index, "chr5", 149526621)
+# it will return
+# 1-element Array{Any,1}:
+#  Dict{ASCIIString,Any}("gene"=>"PDGFRB","number"=>1,"transcript"=>"ENST00000261799.4","type"=>"intron")
+```