dpryan79

rawFastqs = Channel.fromFilePairs(['data/*_{1,2}.fastq.gz', 'data/*_R{1,2}.fastq.gz'])
databases = Channel.of(['label1', '/some/path/database'], ['label2', '/some/other/path/database'], ...)

/+
The general structure of what I'd like to achieve, just without manually putting in all of the someProcess() calls
That won't work for a variety of reasons, among them being that someProcess() would need an alias each time
*/
workflow {
  someProcess(rawFastqs, databases.first())
  someProcess(someProcess.out.fastq, databases[2])

dpryan79

#!/usr/bin/env python
import pyBigWig
import argparse
import os.path
import sys

parser = argparse.ArgumentParser(description="Subset a single chromosome from a bigWig file.")
parser.add_argument("--headerToo", action="store_true", help="Subset the header too, which can expedite some downstream programs but possibly cause problems for others.")
parser.add_argument("input", help="Input bigWig file.")
parser.add_argument("output", help="Output bigWig file.")

dpryan79

#!/usr/bin/env python
# This can take up to ~20 minutes and use up to ~2GB RAM for mammals
import argparse
import py2bit
from deeptoolsintervals import GTF, tree
from deeptoolsintervals.parse import openPossiblyCompressed, parseExonBounds, findRandomLabel

parser = argparse.ArgumentParser(description="Generate a blacklist file of polyX stretches of a given minimum length not within a specified distance of a TES")
parser.add_argument("--output", "-o", help="Output file", required=True)
parser.add_argument("--tb", help="2bit file", required=True)

dpryan79

This demonstrates the issue discussed in https://bitbucket.org/snakemake/snakemake/issues/916/integrated-r-scripts-and-directory

To reproduce that and test any fixes, simply clone this github respository and chmod -R a-w it. Then go to a different directory and snakemake -s path/to/the/snakefile/here. With an unpatched copy of snakemake you'll get a permissions error. Note that the Snakefile has examples using R and Rmd files, so you'll need to ensure that R is in your PATH.

dpryan79

#!/usr/bin/env python
import argparse
import pyBigWig

parser = argparse.ArgumentParser(description="Convert the chromosome names of a bigWig file.")
parser.add_argument("conv", metavar="conversion.txt", help="Text file with two columns, the first a chromosome name and the second the converted chromosome name.")
parser.add_argument("input", metavar="input.bigWig", help="Input bigWig file")
parser.add_argument("output", metavar="output.bigWig", help="Output bigWig file name")
args = parser.parse_args()

dpryan79

#!/usr/bin/env python
import argparse
import numpy as np
import sys

from deeptools import parserCommon, bamHandler, utilities
from deeptools.mapReduce import mapReduce
from deeptools._version import __version__
import deeptools.config as cfg

dpryan79

#!/usr/bin/env python
import deeptools.utilities
from deeptools import bamHandler
from deeptools import mapReduce
from deeptools import writeBedGraph
from deeptools import utilities
from deeptools.utilities import getCommonChrNames, toBytes

import sys
import shutil

dpryan79

#!/usr/bin/env python3
import urllib.request
import json
from subprocess import Popen
from time import sleep


def checkUp(path, timeout):
    p = Popen(['ls', path])
    try:

dpryan79

#!/usr/bin/env python
e2g = []
e2u = []
g2e = []
g2u = []
u2e = []
u2g = []
ifile = open("GCA_000001405.22_GRCh38.p7_assembly_report.txt", "r")
for line in ifile:
    if line.startswith("#"):

dpryan79

#!/usr/bin/env python
# This is just a quick hack and hasn't been tested, but it'll likely work
import pysam
import argparse
import gzip

def writeAlignment(of, al):
    name = al.query_name
    seq = al.query_sequence
    qual = al.query_qualities