TopDownCrawl is a tool for aligning quantitative binding data for k-mers from experiments such as SELEX-seq and SMiLE-seq
Project description
TopDownCrawl
Top-Down-Crawl is able to rapidly align quantitative binding data from experiments such as SELEX-seq and SMiLE-seq. Methods like these are able to provide highly reprodible k-mer level measurements of binding affinity, but understanding where an enriched k-mer falls along the bidning site depends on alignment. k-mer level data is more informative than simplistic PWM-based models which assume independent contributions between positions of the binding site.
The python script (TDC.py) accepts a table (.txt, .tsv, .csv, .xlsx, .xls) as input including the sequences in column 1, and the associated binding measurements in column 2. Values for reverse complements will be averaged if provided, or assumed to be equivalent if not. The first row is expected to contain header names of your choosing. An example input file is included (example_MEF2B_R1_E_k7.tsv).
STEP 1: Install Package with Dependencies
pip install TopDownCrawl
STEP 2: Run Top-Down Crawl
# Accepted formats (.txt, .tsv, .csv, .xlsx, .xls)
TopDownCrawl <input>
ie: TopDownCrawl MEF2B_R1_k7_enrichment.tsv (can be found in example folder)
OUTPUT
# Aligned and reverse-complement-averaged input (tab-delimited)
<input>_aligned.tsv
# Show how many sequences align to each shift
<input>_summary.tsv
# Position weight matrix (PWM) of the aligned input, weighting each sequence by its associated binding metric, removing the 2 most extreme shifts
<input>_PWM.png
# Position weight matrix (PWM) of the aligned input, weighting each sequence by its associated binding metric, removing the 2 most extreme shifts
<input>_PWM_rc.png
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distributions
Built Distribution
Hashes for TopDownCrawl-1.0.1-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 86f134825d4fc0fd99f1080c16678baed6d4ff257705031b5011e6c017e3e2c0 |
|
MD5 | ea9ab2b29311f242865b34a95668ad35 |
|
BLAKE2b-256 | bd0fb3f05669b7ad56ef7ae78e20f1974d335cb38ce16016e5f4b10a1a49ffef |