bs_assign
Find groups of motifs around genomic landmarks of interest
Science Score: 54.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
✓CITATION.cff file
Found CITATION.cff file -
✓codemeta.json file
Found codemeta.json file -
✓.zenodo.json file
Found .zenodo.json file -
○DOI references
-
✓Academic publication links
Links to: zenodo.org -
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (9.7%) to scientific vocabulary
Keywords
Repository
Find groups of motifs around genomic landmarks of interest
Basic Info
Statistics
- Stars: 0
- Watchers: 1
- Forks: 0
- Open Issues: 0
- Releases: 1
Topics
Metadata Files
README.md
bs_assign
Author: aram.amalietti@gmail.com
Dependencies (these are the versions the script was developed with, newer versions should work, but if they don't, please use these versions):
python=3.7
pandas=1.2.3
numpy=1.19.2
pybedtools=0.8.1
iteration_utilities=0.11.0
Usage:
python3 <path_to_script> <motif1,motif2,...motifn> <xl_in> <prtxn> <fasta> <fai> <window> <len> <cores> <chunk_size> <output_dir> <consensus>
motif1,motif2,...motifn is the group of motifs that is searched for around the landmarks, for example AAAA,CCCC,GGGG;
xl_in is a BED file with landmarks around which the motifs are being searched for;
prtxn is the path to the file containing relevant positions where the motif needs to be in order to be detected, it is an output from another script, PEKA;
fasta is the path to the genome in fasta format;
fai is the path to the genome index file;
window is the flanking distance in bases around the landmarks (one base) within which the motifs are being searched for (30 is the usual value);
len is the length of the motifs (4-7 is the usual range);
cores is the number of threads used in the process;
chunk_size is the number of rows per thread (10000 is the usual value);
output_dir is the directory where the results will be saved (make sure it exists);
consensus is a meta motif representing the motif group and will be used in output file names;
Optionally this script can be ran without prtxn file and in that case it will use all the motifs found around crosslinks instead of only those on positions defined in prtxn file.
Outputs:
BED files representing the binding sites found within the given window around the given landmarks, binding sites are basically found motifs merged.
There are two merging distances and two corresponding files, one is merged0
Owner
- Name: Ulelab
- Login: ulelab
- Kind: organization
- Location: London
- Repositories: 10
- Profile: https://github.com/ulelab
Citation (CITATION.cff)
cff-version: 0.0.0 message: "If you use this software, please cite it as below." authors: - family-names: "Amalietti" given-names: "Aram Gustav" orcid: "https://orcid.org/0000-0002-9101-6918" title: "bs_assign" version: 0.0.0 doi: 10.5281/zenodo.8388764 date-released: 2021-06-30 url: "https://github.com/ulelab/bs_assign"