Skip to content

Commit

Permalink
Update README.md
Browse files Browse the repository at this point in the history
  • Loading branch information
VincentGardeux authored Apr 6, 2018
1 parent e2f0813 commit df752ca
Showing 1 changed file with 11 additions and 11 deletions.
22 changes: 11 additions & 11 deletions README.md
Original file line number Diff line number Diff line change
@@ -1,13 +1,13 @@
![](https://img.shields.io/badge/build-passing-green.svg)
![](https://img.shields.io/badge/version-1.1-blue.svg)
![](https://img.shields.io/badge/version-1.2-blue.svg)
![](https://img.shields.io/badge/picard-2.9.0-blue.svg)
![](https://img.shields.io/badge/java-1.8-red.svg)

# BRB-seq Tools 1.1
# BRB-seq Tools 1.2
A suite of tools for the pre-processing of BRB-seq data (bulk RNA-seq)

## Download software
BRB-seq command-line tools are provided as a [single executable jar file](../master/releases/BRBseqTools.1.1.jar?raw=true).
BRB-seq command-line tools are provided as a [single executable jar file](../master/releases/BRBseqTools.1.2.jar?raw=true).
The .jar file contains all required materials and can be run on any terminal.

## Dependencies
Expand Down Expand Up @@ -39,7 +39,7 @@ For further analyses (filtering, normalization, dimension reduction, clustering,
To check that BRB-seq Tools is working properly, run the following command:

```bash
java -jar BRBseqTools.1.1.jar
java -jar BRBseqTools.1.2.jar
```
As shown in the output of this command, BRB-seq tool suite allows the user to run 3 possible tools.

Expand Down Expand Up @@ -71,11 +71,11 @@ Options:

Example:
```bash
java -jar BRBseqTools.1.1.jar Demultiplex -r1 lib_example_R1.fastq.gz -r2 lib_example_R2.fastq.gz -c lib_example_barcodes.txt -p BU -UMI 14
java -jar BRBseqTools.1.2.jar Demultiplex -r1 lib_example_R1.fastq.gz -r2 lib_example_R2.fastq.gz -c lib_example_barcodes.txt -p BU -UMI 14
```
or, if no UMI:
```bash
java -jar BRBseqTools.1.1.jar Demultiplex -r1 lib_example_R1.fastq.gz -r2 lib_example_R2.fastq.gz -c lib_example_barcodes.txt -p B
java -jar BRBseqTools.1.2.jar Demultiplex -r1 lib_example_R1.fastq.gz -r2 lib_example_R2.fastq.gz -c lib_example_barcodes.txt -p B
```

> **Note:** When using this tool, UMIs are kept as indexes in all output .fastq files
Expand Down Expand Up @@ -106,7 +106,7 @@ Options:

Example:
```bash
java -jar BRBseqTools.1.1.jar CreateDGEMatrix -f lib_example_R1.fastq.gz -b lib_example_R2.bam -c lib_example_barcodes.txt -gtf Homo_sapiens.GRCh38.90.gtf.gz -p BU -UMI 14
java -jar BRBseqTools.1.2.jar CreateDGEMatrix -f lib_example_R1.fastq.gz -b lib_example_R2.bam -c lib_example_barcodes.txt -gtf Homo_sapiens.GRCh38.90.gtf.gz -p BU -UMI 14
```

> **Note:** The original BRB-seq protocol contains a UMI construct. But you can also generate a library without UMIs, or even not sequence the UMIs from R2 read. If UMIs are present, both UMI and read count matrices will be generated. If not, only the read count table will be generated.
Expand Down Expand Up @@ -136,7 +136,7 @@ Options:

Example:
```bash
java -jar BRBseqTools.1.1.jar AnnotateBAM -f lib_example_R1.fastq.gz -b lib_example_R2.bam -p BU -UMI 14
java -jar BRBseqTools.1.2.jar AnnotateBAM -f lib_example_R1.fastq.gz -b lib_example_R2.bam -p BU -UMI 14
```

> **Note:** Most of the options are optionals here. If not put, then the corresponding tag will simply not be added to the annotated BAM.
Expand All @@ -156,7 +156,7 @@ Options:

Example:
```bash
java -jar BRBseqTools.1.1.jar Trim -f lib_example_R2.fastq.gz
java -jar BRBseqTools.1.2.jar Trim -f lib_example_R2.fastq.gz
```

> **Note:** If you use STAR for alignment, this step is optional, as it will not change much the results of the alignment (our tests have shown that the improvement is real but very minor)
Expand Down Expand Up @@ -184,7 +184,7 @@ STAR --runMode genomeGenerate --genomeDir STAR_Index/ --genomeFastaFiles Homo_sa
When the index is built, you will never need to rebuild it. Then you can use only the following script:
```bash
# (Optional) Trim the read containing the sequence fragments (generates a 'lib_example_R2.trimmed.fastq.gz' file)
java -jar BRBseqTools.1.1.jar Trim -f lib_example_R2.fastq.gz
java -jar BRBseqTools.1.2.jar Trim -f lib_example_R2.fastq.gz
# Create output folder
mkdir BAM/
# Align only the R2 fastq file (using STAR, no sorting/indexing is needed)
Expand All @@ -194,7 +194,7 @@ STAR --runMode alignReads --genomeDir STAR_Index/ --outFilterMultimapNmax 1 --re
# (optional) Rename the output aligned BAM
mv BAM/Aligned.out.bam BAM/lib_example_R2.bam
# Demultiplex and generate output count/UMI matrix
java -jar BRBseqTools.1.0.jar CreateDGEMatrix -f lib_example_R1.fastq.gz -b BAM/lib_example_R2.bam -c lib_example_barcodes.txt -gtf Homo_sapiens.GRCh38.90.gtf -p BU -UMI 14
java -jar BRBseqTools.1.2.jar CreateDGEMatrix -f lib_example_R1.fastq.gz -b BAM/lib_example_R2.bam -c lib_example_barcodes.txt -gtf Homo_sapiens.GRCh38.90.gtf -p BU -UMI 14
# Note: This example suppose that R1 has barcode followed by 14bp UMI (see above for other cases)
# Note: 'lib_example_barcodes.txt' should be created by the user and should contain the mapping between the barcode and the sample name¹
```
Expand Down

0 comments on commit df752ca

Please sign in to comment.