Skip to content

Commit

Permalink
Updated commands to match modifications in setup.py
Browse files Browse the repository at this point in the history
  • Loading branch information
Dana Elizabeth Wyman committed Sep 17, 2019
1 parent f98b82a commit 76b62a4
Showing 1 changed file with 10 additions and 7 deletions.
17 changes: 10 additions & 7 deletions example/README.md
Original file line number Diff line number Diff line change
@@ -1,21 +1,24 @@
# TALON example

## User note
This example was initially constructed for TALON v4.2. In the packaged versions since, the commands have changed since TALON does not need to be run from a specific path. If you are using a TALON version below 4.3, please see the archived documentation here for the correct commands: https://github.com/dewyman/TALON/wiki/Archived-TALON-Example-Instructions-(v4.2).

## Tutorial
Use the provided files to try TALON out on chromosome 21 PacBio SAM reads from two biological replicates of human cell line GM12878.

First, initialize the TALON database from the provided GTF annotation of chromosome 21:

```
python ../initialize_talon_database.py \
talon_initialize_database \
--f gencode.v29.chr21.gtf \
--a gencode_v29 \
--g hg38 \
--o example_talon
```



Now, use the provided config file and the newly initialized database to annotate and quantify the reads. The database is modified in place.
```
python ../talon.py \
talon \
--f config.csv \
--db example_talon.db \
--build hg38 \
Expand All @@ -25,15 +28,15 @@ The file 'example_talon_QC.log' contains a record of each input transcript along

Now, we can explore the results of the run. To summarize how many of each transxcript were found (prior to any filtering), run:
```
python ../post-TALON_tools/summarize_datasets.py \
talon_summarize \
--db example_talon.db \
--v \
--o example
```

To create an abundance matrix without filtering (for use computing gene expression), we run the following:
```
python ../post-TALON_tools/create_abundance_file_from_database.py \
talon_abundance \
--db example_talon.db \
-a gencode_v29 \
--build hg38 \
Expand All @@ -42,7 +45,7 @@ python ../post-TALON_tools/create_abundance_file_from_database.py \

We can also generate a filtered abundance matrix where we require transcripts to be either a) known, or b) detected in both replicates (for transcript-level expression):
```
python ../post-TALON_tools/create_abundance_file_from_database.py \
talon_abundance \
--db example_talon.db \
-a gencode_v29 \
--filter \
Expand Down

0 comments on commit 76b62a4

Please sign in to comment.