Skip to content

Commit ca2401e

Browse files
authored
Update README.md
1 parent c40d6b7 commit ca2401e

1 file changed

Lines changed: 4 additions & 2 deletions

File tree

README.md

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -26,7 +26,7 @@ pip install -e .
2626

2727
## Usage
2828
### Peak-based applications (ChIP-seq, ATAC-seq, DNase-seq, etc)
29-
#### Pre-processing
29+
#### Pre-processing and alignment
3030
Using Allo requires a few pre-processing steps. In most ChIP-seq, ATAC-seq, and DNase-seq pipelines, the default behavior of aligners is to assign multi-mapped reads to random locations within their mappings without retaining information on the other locations. Both Bowtie1/2 and BWA can be used for single-end. Unfortunately, BWA cannot be used for paired-end reads prior to Allo due to constraints in how it outputs multi-mapped reads. The following arguments should be used:
3131

3232
*Bowtie1*
@@ -75,12 +75,14 @@ allo testRunPE.sam -seq pe
7575
#### Additional tips
7676
It is recommended to run Allo on both the control and target sequencing files in order to balance out background in the samples. We recommend running Allo using the --random argument on the control file. This generally results in higher confidence peaks.
7777

78-
### Pre-processing for RNA-seq
78+
### RNA-seq application
79+
#### Pre-processing and alignment
7980
Allo is compatible with STAR alignments. We recommend using the "--outFilterType BySJout" argument if you choose to use the "--splice" function in Allo in order to only consider high quality junctions. An example of a paired-end STAR alignment keeping up to 25 locations per read is shown below:
8081
```
8182
STAR --genomeDir GENOMEDIR --readFilesIn fASTQ_1 FASTQ_2 --outSAMtype BAM Unsorted --outSAMmultNmax 25 --outFilterType BySJout --outFileNamePrefix ALIGNEROUTPUT
8283
```
8384

85+
#### Running Allo
8486
To use Allo, first sort your file:
8587
```
8688
samtools collate -o ALIGNEROUTPUT_SORT.BAM ALIGNEROUTPUT_FILTER.BAM

0 commit comments

Comments
 (0)