Skip to content

Commit

Permalink
Update README.md
Browse files Browse the repository at this point in the history
  • Loading branch information
ArtPoon authored Sep 10, 2019
1 parent 1819eb2 commit 40b57b1
Showing 1 changed file with 4 additions and 1 deletion.
5 changes: 4 additions & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -101,7 +101,7 @@ Successfully installed sierralocal-0.2.1
To demonstrate the use of *sierra-local* to predict drug resistance levels
from next-generation sequence data processed by *MiCall-Lite*, we'll use a published
NGS data set from a
[study of HIV-1 drug resistance in Uganda](https://www.liebertpub.com/doi/abs/10.1089/AID.2017.0205).
[study of integrase-inhibitor treatment failure in Uganda](https://www.liebertpub.com/doi/abs/10.1089/AID.2017.0205).
(To obtain this data set, you'll need to install NCBI's
[sratools](https://ncbi.github.io/sra-tools/), which provides the *fasterq-dump*
file transfer program).
Expand Down Expand Up @@ -147,6 +147,9 @@ In this output, the `header` column indicates how nucleotide polymorphisms in th
`0.010` indicates that any nucleotide polymorphism at a frequency above 1% should be incorporated into the consensus sequences, resulting in a mixture or ambiguous base call.
Lowering this frequency threshold increases the number of mutations that are available for resistance scoring, as demonstrated by the longer list of mutations in the second row in this example.

Note that the Sierra pipeline automatically detects which *pol* gene products are covered in the NGS data.
In this example, the reads only span HIV-1 integrase and so only resistance scores associated with integrase inhibitors are reported.


## How it works

Expand Down

0 comments on commit 40b57b1

Please sign in to comment.