Calculate TSS enrichment for ATAC-seq data
Project description
tssenrich
Calculate TSS enrichment for ATAC-seq data. (see also the ENCODE standards)
Installation
pip3 install tssenrich
or
pip3 install --user tssenrich
Usage
usage: tssenrich [-h] [--genome {hg38,hg19,mm10}] [--memory <float>]
[--processes <int>] [--mapping-quality <int>]
[--samtools-path <path/to/samtools>]
[--log <path/to/log.txt>] [--tmp-dir <temp/file/dir/>]
[--flank-distance <int>] [--flank-size <int>]
<path/to/file.bam> [<path/to/file.bam> ...]
calculate TSS enrichment for ATAC-seq data
positional arguments:
<path/to/file.bam> Path to input BAM file
optional arguments:
-h, --help show this help message and exit
--genome {hg38,hg19,mm10}
genome build [hg38]
--memory <float> memory limit in GB [2]
--processes <int> number of processes/threads to use [1]
--mapping-quality <int>
ignore reads with mapping quality below the given value [0]
--samtools-path <path/to/samtools>
path to an alternate samtools executable [/usr/local/bin/samtools]
--log <path/to/log.txt>
path to log file
--tmp-dir <temp/file/dir/>
directory to use for temporary files
--flank-distance <int>
distance from tss of outer ends of flanks [1000]
--flank-size <int> size of flanks (for determining average depth) [100]
ENCODE standards:
| Genome | Concerning | Acceptable | Ideal |
| ------ | ---------- | ---------- | ----- |
| hg19 | < 6 | 6 - 10 | > 10 |
| hg38 | < 5 | 5 - 7 | > 7 |
| mm10 | < 10 | 10 - 15 | > 15 |
Examples
tssenrich example.bam > score.txt
tssenrich --genome hg19 --log log.txt --memory 2 --processes 2 example1.bam example2.bam > scores.txt
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
tssenrich-1.2.2.tar.gz
(17.9 MB
view hashes)
Built Distribution
tssenrich-1.2.2-py3-none-any.whl
(17.9 MB
view hashes)
Close
Hashes for tssenrich-1.2.2-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 92754083d331955bf5f5e396f5ee31ab0aa6f57762138cac40a63781d48d0922 |
|
MD5 | e060bf0d5f3758621e5c982316af1cc1 |
|
BLAKE2b-256 | f0ba4dd8a2249b94ed6d838095c5e944a7add42a7b8ee894ccc5298f8167d4c2 |