Skip to content

Instantly share code, notes, and snippets.

View Phlya's full-sized avatar

Ilya Flyamer Phlya

  • FMI
  • Basel, Switzerland
  • 05:19 (UTC +02:00)
  • X @phlya
View GitHub Profile
@Phlya
Phlya / uneuploidy_project.yml
Last active July 11, 2023 14:23
uneuploidy_project.yml
#########################################
# THIS IS A TYPICAL project.yml TEMPLATE
# most of the settings present here
# are GO for mapping production data
# but nonetheless user must consider
# carefully every presented option
#########################################
#########################################
# When commmenting parameters out, make sure
(test-snakemake-wrappers) ilya@f346m-7eff66 snakemake-wrappers % pytest test.py -v -k cooltools
======================================================================= test session starts ========================================================================
platform darwin -- Python 3.10.5, pytest-7.1.2, pluggy-1.0.0 -- /Users/ilya/opt/anaconda3/envs/test-snakemake-wrappers/bin/python
cachedir: .pytest_cache
rootdir: /Users/ilya/Projects/snakemake-wrappers
collected 380 items / 379 deselected / 1 selected
test.py::test_cooltools_insulation FAILED [100%]
============================================================================= FAILURES =============================================================================
@Phlya
Phlya / Snakefile
Last active March 28, 2022 08:49
Pairtools dedup benchmarking
walkpolicies = ["all", "mask"]
backends = ["cython", "scipy", "sklearn"]
cores = [1, 2, 4, 8]
chunksizes = [1000, 10_000, 100_000, 1_000_000, 10_000_000]
mismatches = [0, 3]
carryovers = [10, 100, 1000]
extensions = ["nodups.pairs", "dups.pairs", "unmapped.pairs", "stats"]
cython = expand(
@Phlya
Phlya / expected_fetching.ipynb
Last active March 1, 2021 12:58
fetching arbitrary expected areas
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@Phlya
Phlya / project.yml
Created December 15, 2020 16:35
Distiller project for Hi-C processing of Phanstiel et al. 2017 data
# Fastqs can be provided as:
# -- a pairs of relative/absolute paths
# -- sra:<SRA_NUMBER>, optionally followed by the indices of the first and
# the last entry in the SRA in the form of "?start=<first>&end=<last>
# [to implement] -- as a path to a folder with fastqs '<base_folder>', with the structure
# <base_folder>/<library_name>/<run_name>/, with each folder containing only
# two fastq.gz files
# python bin/geo2yaml.py PRJNA385337 --title_column library_name --title_sub '\.[^.]*$' '' --title_sub '\.' '_' --group_sub '_[^_]+$' ''
input:
raw_reads_paths:
import sys
import argparse
import re
import numpy as np
import pandas as pd
import pysradb
# def unescaped_str(arg_str):
@Phlya
Phlya / quaich_config.yaml
Created September 11, 2020 15:28
Quaich config used for MCM paper bulk Hi-C data analsyis
genome: hg38
path_genomes: /exports/igmm/eddie/wendy-lab/ilia/genomes/
project_folder: /exports/igmm/eddie/wendy-lab/ilia/MCM/
filter: mapq_30
resolutions:
- 10000000
- 5000000
@Phlya
Phlya / Quaich_example.ipynb
Last active August 25, 2020 18:02
Example notebook for plotting quaich output
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@Phlya
Phlya / quaich_example_config.yaml
Last active August 26, 2020 16:15
Example config file for quaich
genome: hg38
path_genomes: /exports/igmm/eddie/wendy-lab/ilia/genomes/
project_folder: /exports/igmm/eddie/wendy-lab/ilia/quaich_example/
filter: mapq_30
resolutions:
- 10000000
- 5000000
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.