Skip to content

Instantly share code, notes, and snippets.

@inutano
Created February 23, 2017 05:24
Show Gist options
  • Save inutano/843de747a735bd5011405ef511cb3f3f to your computer and use it in GitHub Desktop.
Save inutano/843de747a735bd5011405ef511cb3f3f to your computer and use it in GitHub Desktop.
# Calculate average throughput for long read sequencing data
PREFIX sos: <http://purl.jp/bio/01/quanto/ontology/sos#>
PREFIX quanto: <http://purl.jp/bio/01/quanto/resource/>
PREFIX rdfs: <http://www.w3.org/2000/01/rdf-schema#>
PREFIX dct: <http://purl.org/dc/terms/>
PREFIX pav: <http://purl.org/pav/>
SELECT (AVG(?total_sequences AS ?average_throughput))
FROM <http://quanto.dbcls.jp>
WHERE {
?quanto dct:identifier ?id;
sos:totalSequences ?x1;
sos:overallMedianBaseCallQuality ?x2;
sos:medianSequenceLength ?x3 .
?x1 rdf:value ?total_sequences .
?x2 rdf:value ?median_bc_quality .
?x3 rdf:value ?median_seq_length .
FILTER(?median_seq_length > 500).
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment