Following on from the first blog post about htsget-rs and Crypt4GH, this post goes into further details about how htsget works and illustrates more complex use cases.

Let's start by querying an example file. The deployed GA4GH htsget instance has access to example files from the htsget-rs repository. Recall from the first htsget-rs blog post that the reads endpoint can serve BAM files, for example:

curl "https://htsget.ga4gh-demo.org/reads/htsnexus_test_NA12878"

This will return a set of URL "tickets" inside the "urls" field of the JSON response. These "tickets" contain URLs that should be fetched and concatenated to produce the response. Additionally, there is a "headers" field that contains HTTP headers that should included when requesting the url in the ticket. Take a look at the htsget spec for more details.

To simplify fetching and concatenating URL tickets, use a htsget client, such as the GA4GH client.

As a simple example, query the header of the file example file by passing class=Header:

htsget "https://htsget.ga4gh-demo.org/reads/htsnexus_test_NA12878?class=header" > out.bam

Internally, this yields a JSON with a URL that can be fetched along with a "Range" header:

{
  "htsget": {
    "format": "BAM",
    "urls": [
      {
        "url": "...",
        "headers": {
          "Range": "bytes=0-4667"
        },
        "class": "header"
      }
    ]
  }
}

The client takes care of fetching the URLs and concatenating bytes.

A strength of the htsget protocol is that the output represents a small part of the full file, allowing the user to query specific regions of a file without needing to obtain the entire file.

In this case, the output represents the BAM header of the file:

samtools view -H out

Note that "..." inside the JSON example responses represents some data or a URL. This will be different when executing the query.

A more interesting query would involve selecting a specific region, for example chr11. This can be accomplished by using the referenceName parameter. Viewing the output will show data for that specific region:

htsget "https://htsget.ga4gh-demo.org/reads/htsnexus_test_NA12878?referenceName=11" | samtools view

Similarly, the query can be refined further by specifying specific start and end ranges, so that only those regions are returned:

htsget "https://htsget.ga4gh-demo.org/reads/htsnexus_test_NA12878?referenceName=11&start=500000&end=5001000" | samtools view

Internally, the output from htsget-rs will contain multiple URL tickets that represent the specific data queried:

{
  "htsget": {
    "format": "BAM",
    "urls": [
      {
        "url": "...",
        "headers": {
          "Range": "bytes=0-273085"
        }
      },
      {
        "url": "...",
        "headers": {
          "Range": "bytes=499249-574358"
        }
      },
      {
        "url": "...",
        "headers": {
          "Range": "bytes=627987-647345"
        }
      },
      {
        "url": "...",
        "headers": {
          "Range": "bytes=824361-842100"
        }
      },
      {
        "url": "...",
        "headers": {
          "Range": "bytes=977196-996014"
        }
      },
      {
        "url": "...",
        "headers": {
          "Range": "bytes=2596771-2596798"
        }
      }
    ]
  }
}

Moving on to a more complex example, we will now incorporate querying Crypt4GH encrypted files from htsget-rs. To decrypt Crypt4GH files, install the Crypt4GH CLI and get the keys from the htsget-rs repository.

Then, query like before, except add the encryptionScheme parameter:

curl "https://htsget.ga4gh-demo.org/reads/htsnexus_test_NA12878?class=header&encryptionScheme=C4GH"

This is an experimental parameter that htsget-rs supports but is not yet part of the htsget spec.

This will return a JSON that contains encrypted data when concatenated. Here, there are additional URLs that are base64 encoded. These URLs represent inline data to the JSON ticket, and just need to be decoded to obtain the bytes. They follow the same semantics as the other URLs and should be concatenated after decoding.

{
  "htsget": {
    "format": "BAM",
    "urls": [
      {
        "url": "data:;base64,..."
      },
      {
        "url": "...",
        "headers": {
          "Range": "bytes=16-123"
        }
      },
      {
        "url": "data:;base64,..."
      },
      {
        "url": "...",
        "headers": {
          "Range": "bytes=124-65687"
        }
      }
    ]
  }
}

Putting it all together, and using the keys from the htsget-rs repo, the data can be accessed by running:

htsget "https://htsget.ga4gh-demo.org/reads/htsnexus_test_NA12878?class=header&encryptionScheme=C4GH" | crypt4gh decrypt --sk bob.sec | samtools view -H

As an extra section, the htsget protocol has a compliance suite that can be run on htsget-rs. This contains tests that ensure that htsget-rs runs as expected.

In order to run the compliance tests, follow the installation instructions in the htsget-compliance repository and then run the following on the deployed htsget-rs instance:

htsget-compliance https://htsget.ga4gh-demo.org | jq '.["summary"]'