diff options
Diffstat (limited to 'doc/blog/using-covid-19-pubseq-part1.html')
-rw-r--r-- | doc/blog/using-covid-19-pubseq-part1.html | 82 |
1 files changed, 50 insertions, 32 deletions
diff --git a/doc/blog/using-covid-19-pubseq-part1.html b/doc/blog/using-covid-19-pubseq-part1.html index 0e6136c..5fd86d1 100644 --- a/doc/blog/using-covid-19-pubseq-part1.html +++ b/doc/blog/using-covid-19-pubseq-part1.html @@ -3,7 +3,7 @@ "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"> <html xmlns="http://www.w3.org/1999/xhtml" lang="en" xml:lang="en"> <head> -<!-- 2020-07-17 Fri 05:05 --> +<!-- 2020-07-19 Sun 02:32 --> <meta http-equiv="Content-Type" content="text/html;charset=utf-8" /> <meta name="viewport" content="width=device-width, initial-scale=1" /> <title>COVID-19 PubSeq (part 1)</title> @@ -248,20 +248,20 @@ for the JavaScript code in this tag. <h2>Table of Contents</h2> <div id="text-table-of-contents"> <ul> -<li><a href="#org0db5db0">1. What does this mean?</a></li> -<li><a href="#orge5267fd">2. Fetch sequence data</a></li> -<li><a href="#orgfbd3adc">3. Predicates</a></li> -<li><a href="#org08e70e1">4. Fetch submitter info and other metadata</a></li> -<li><a href="#org9194557">5. Fetch all sequences from Washington state</a></li> -<li><a href="#org76317ad">6. Discussion</a></li> -<li><a href="#orgeb871a1">7. Acknowledgements</a></li> +<li><a href="#orgb852bf7">1. What does this mean?</a></li> +<li><a href="#orge6db105">2. Fetch sequence data</a></li> +<li><a href="#orgf3b8001">3. Predicates</a></li> +<li><a href="#org11097b0">4. Fetch submitter info and other metadata</a></li> +<li><a href="#org4f8467e">5. Fetch all sequences from Washington state</a></li> +<li><a href="#orge9b18e2">6. Discussion</a></li> +<li><a href="#orga0badf8">7. Acknowledgements</a></li> </ul> </div> </div> -<div id="outline-container-org0db5db0" class="outline-2"> -<h2 id="org0db5db0"><span class="section-number-2">1</span> What does this mean?</h2> +<div id="outline-container-orgb852bf7" class="outline-2"> +<h2 id="orgb852bf7"><span class="section-number-2">1</span> What does this mean?</h2> <div class="outline-text-2" id="text-1"> <p> This means that when someone uploads a SARS-CoV-2 sequence using one @@ -313,9 +313,8 @@ initiative! </div> </div> - -<div id="outline-container-orge5267fd" class="outline-2"> -<h2 id="orge5267fd"><span class="section-number-2">2</span> Fetch sequence data</h2> +<div id="outline-container-orge6db105" class="outline-2"> +<h2 id="orge6db105"><span class="section-number-2">2</span> Fetch sequence data</h2> <div class="outline-text-2" id="text-2"> <p> The latest run of the pipeline can be viewed <a href="https://workbench.lugli.arvadosapi.com/collections/lugli-4zz18-z513nlpqm03hpca">here</a>. Each of these @@ -339,8 +338,8 @@ these identifiers throughout. </div> </div> -<div id="outline-container-orgfbd3adc" class="outline-2"> -<h2 id="orgfbd3adc"><span class="section-number-2">3</span> Predicates</h2> +<div id="outline-container-orgf3b8001" class="outline-2"> +<h2 id="orgf3b8001"><span class="section-number-2">3</span> Predicates</h2> <div class="outline-text-2" id="text-3"> <p> To explore an RDF dataset, the first query we can do is open and gets @@ -446,15 +445,18 @@ select (COUNT(distinct ?dataset) as ?num) } </pre> </div> + +<p> +Run this <a href="http://sparql.genenetwork.org/sparql/?default-graph-uri=&query=PREFIX+pubseq%3A+%3Chttp%3A%2F%2Fbiohackathon.org%2Fbh20-seq-schema%23MainSchema%2F%3E%0D%0Aselect+%28COUNT%28distinct+%3Fdataset%29+as+%3Fnum%29%0D%0A%7B%0D%0A+++%3Fdataset+pubseq%3Asubmitter+%3Fid+.%0D%0A+++%3Fid+%3Fp+%3Fsubmitter%0D%0A%7D&format=text%2Fhtml&timeout=0&debug=on&run=+Run+Query+">query</a>. +</p> </div> </div> - -<div id="outline-container-org08e70e1" class="outline-2"> -<h2 id="org08e70e1"><span class="section-number-2">4</span> Fetch submitter info and other metadata</h2> +<div id="outline-container-org11097b0" class="outline-2"> +<h2 id="org11097b0"><span class="section-number-2">4</span> Fetch submitter info and other metadata</h2> <div class="outline-text-2" id="text-4"> <p> -To get dataests with submitters we can do the above +To get datasets with submitters we can do the above </p> <div class="org-src-container"> @@ -468,6 +470,10 @@ select distinct ?dataset ?p ?submitter </div> <p> +Run this <a href="http://sparql.genenetwork.org/sparql/?default-graph-uri=&query=PREFIX+pubseq%3A+%3Chttp%3A%2F%2Fbiohackathon.org%2Fbh20-seq-schema%23MainSchema%2F%3E%0D%0Aselect+distinct+%3Fdataset+%3Fp+%3Fsubmitter%0D%0A%7B%0D%0A+++%3Fdataset+pubseq%3Asubmitter+%3Fid+.%0D%0A+++%3Fid+%3Fp+%3Fsubmitter%0D%0A%7D&format=text%2Fhtml&timeout=0&debug=on&run=+Run+Query+">query</a>. +</p> + +<p> Tells you one submitter is "Roychoudhury,P.;Greninger,A.;Jerome,K." with a URL <a href="http://purl.obolibrary.org/obo/NCIT_C42781">predicate</a> (<a href="http://purl.obolibrary.org/obo/NCIT_C42781">http://purl.obolibrary.org/obo/NCIT_C42781</a>) explaining "The individual who is responsible for the content of a @@ -526,6 +532,10 @@ select distinct ?sid ?sample ?p1 ?dataset ?submitter </div> <p> +Run <a href="http://sparql.genenetwork.org/sparql/?default-graph-uri=&query=%0D%0APREFIX+pubseq%3A+%3Chttp%3A%2F%2Fbiohackathon.org%2Fbh20-seq-schema%23MainSchema%2F%3E%0D%0Aselect+distinct+%3Fsid+%3Fsample+%3Fp1+%3Fdataset+%3Fsubmitter%0D%0A%7B%0D%0A+++%3Fdataset+pubseq%3Asubmitter+%3Fid+.%0D%0A+++%3Fid+%3Fp+%3Fsubmitter+.%0D%0A+++FILTER%28CONTAINS%28%3Fsubmitter%2C%22Roychoudhury%22%29%29+.%0D%0A+++%3Fdataset+pubseq%3Asample+%3Fsid+.%0D%0A+++%3Fsid+%3Fp1+%3Fsample%0D%0A%7D%0D%0A&format=text%2Fhtml&timeout=0&debug=on&run=+Run+Query+">query</a>. +</p> + +<p> which shows pretty much <a href="http://sparql.genenetwork.org/sparql/?default-graph-uri=&query=PREFIX+pubseq%3A+%3Chttp%3A%2F%2Fbiohackathon.org%2Fbh20-seq-schema%23MainSchema%2F%3E%0D%0Aselect+distinct+%3Fsid+%3Fsample+%3Fp1+%3Fdataset+%3Fsubmitter%0D%0A%7B%0D%0A+++%3Fdataset+pubseq%3Asubmitter+%3Fid+.%0D%0A+++%3Fid+%3Fp+%3Fsubmitter+.%0D%0A+++FILTER%28CONTAINS%28%3Fsubmitter%2C%22Roychoudhury%22%29%29+.%0D%0A+++%3Fdataset+pubseq%3Asample+%3Fsid+.%0D%0A+++%3Fsid+%3Fp1+%3Fsample%0D%0A%7D&format=text%2Fhtml&timeout=0&debug=on&run=+Run+Query+">everything known</a> about their submissions in this database. Let's focus on one sample "MT326090.1" with predicate <a href="http://semanticscience.org/resource/SIO_000115">http://semanticscience.org/resource/SIO_000115</a>. @@ -543,21 +553,26 @@ select distinct ?sample ?p ?o </div> <p> -This <a href="http://sparql.genenetwork.org/sparql/?default-graph-uri=&query=PREFIX+pubseq%3A+%3Chttp%3A%2F%2Fbiohackathon.org%2Fbh20-seq-schema%23MainSchema%2F%3E%0D%0APREFIX+sio%3A+%3Chttp%3A%2F%2Fsemanticscience.org%2Fresource%2F%3E%0D%0Aselect+distinct+%3Fsample+%3Fp+%3Fo%0D%0A%7B%0D%0A+++%3Fsample+sio%3ASIO_000115+%22MT326090.1%22+.%0D%0A+++%3Fsample+%3Fp+%3Fo+.%0D%0A%7D&format=text%2Fhtml&timeout=0&debug=on&run=+Run+Query+">query</a> tells us the sample was submitted "2020-03-21" and +Run <a href="http://sparql.genenetwork.org/sparql/?default-graph-uri=&query=%0D%0APREFIX+pubseq%3A+%3Chttp%3A%2F%2Fbiohackathon.org%2Fbh20-seq-schema%23MainSchema%2F%3E%0D%0APREFIX+sio%3A+%3Chttp%3A%2F%2Fsemanticscience.org%2Fresource%2F%3E%0D%0Aselect+distinct+%3Fsample+%3Fp+%3Fo%0D%0A%7B%0D%0A+++%3Fsample+sio%3ASIO_000115+%22MT326090.1%22+.%0D%0A+++%3Fsample+%3Fp+%3Fo+.%0D%0A%7D&format=text%2Fhtml&timeout=0&debug=on&run=+Run+Query+">query</a>. +</p> + +<p> +This query tells us the sample was submitted "2020-03-21" and originates from <a href="http://www.wikidata.org/entity/Q30">http://www.wikidata.org/entity/Q30</a>, i.e., the USA and is a biospecimen collected from the back of the throat by swabbing. -We can track it back to the original GenBank <a href="http://identifiers.org/insdc/MT326090.1#sequence">submission</a>. +We can track it back to the original GenBank <a href="http://identifiers.org/insdc/MT326090.1#sequence">submission</a> using the +<a href="http://identifiers.org/insdc/MT326090.1">http://identifiers.org/insdc/MT326090.1</a> link. </p> <p> We have also added country and label data to make it a bit easier -to view/query the database. +to view/query the database and place the sequence on the <a href="http://covid19.genenetwork.org/">map</a>. </p> </div> </div> -<div id="outline-container-org9194557" class="outline-2"> -<h2 id="org9194557"><span class="section-number-2">5</span> Fetch all sequences from Washington state</h2> +<div id="outline-container-org4f8467e" class="outline-2"> +<h2 id="org4f8467e"><span class="section-number-2">5</span> Fetch all sequences from Washington state</h2> <div class="outline-text-2" id="text-5"> <p> Now we know how to get at the origin we can do it the other way round @@ -574,8 +589,8 @@ and fetch all sequences referring to Washington state </div> <p> -which lists 300 sequences originating from Washington state! Which is almost -half of the set coming out of GenBank. +which lists 300 sequences originating from Washington state! Which in +April was almost half of the set coming out of GenBank. </p> <p> @@ -591,12 +606,15 @@ entity is <a href="https://www.wikidata.org/wiki/Q43">Q43</a>: } </pre> </div> + +<p> +Run <a href="http://sparql.genenetwork.org/sparql/?default-graph-uri=&query=%0D%0Aselect+%3Fseq+%3Fsample%0D%0A%7B%0D%0A++++%3Fseq+%3Chttp%3A%2F%2Fbiohackathon.org%2Fbh20-seq-schema%23MainSchema%2Fsample%3E+%3Fsample+.%0D%0A++++%3Fsample+%3Chttp%3A%2F%2Fpurl.obolibrary.org%2Fobo%2FGAZ_00000448%3E+%3Chttp%3A%2F%2Fwww.wikidata.org%2Fentity%2FQ43%3E%0D%0A%7D&format=text%2Fhtml&timeout=0&debug=on&run=+Run+Query+">query</a>. +</p> </div> </div> - -<div id="outline-container-org76317ad" class="outline-2"> -<h2 id="org76317ad"><span class="section-number-2">6</span> Discussion</h2> +<div id="outline-container-orge9b18e2" class="outline-2"> +<h2 id="orge9b18e2"><span class="section-number-2">6</span> Discussion</h2> <div class="outline-text-2" id="text-6"> <p> The public sequence uploader collects sequences, raw data and @@ -607,8 +625,8 @@ referenced in publications and origins are citeable. </div> </div> -<div id="outline-container-orgeb871a1" class="outline-2"> -<h2 id="orgeb871a1"><span class="section-number-2">7</span> Acknowledgements</h2> +<div id="outline-container-orga0badf8" class="outline-2"> +<h2 id="orga0badf8"><span class="section-number-2">7</span> Acknowledgements</h2> <div class="outline-text-2" id="text-7"> <p> The overall effort was due to magnificent freely donated input by a @@ -623,7 +641,7 @@ Garrison this initiative would not have existed! </div> </div> <div id="postamble" class="status"> -<hr><small>Created by <a href="http://thebird.nl/">Pjotr Prins</a> (pjotr.public768 at thebird 'dot' nl) using Emacs org-mode and a healthy dose of Lisp!<br />Modified 2020-07-17 Fri 05:02</small>. +<hr><small>Created by <a href="http://thebird.nl/">Pjotr Prins</a> (pjotr.public768 at thebird 'dot' nl) using Emacs org-mode and a healthy dose of Lisp!<br />Modified 2020-07-19 Sun 02:32</small>. </div> </body> </html> |