Skip to content

Instantly share code, notes, and snippets.

View marcelm's full-sized avatar

Marcel Martin marcelm

  • Stockholm
  • 22:07 (UTC +02:00)
View GitHub Profile
@marcelm
marcelm / snakemake-pure-python.py
Last active November 29, 2023 00:45
pure Python module that uses snakemake to construct and run a workflow
#!/usr/bin/env python3
"""
Running this script is (intended to be) equivalent to running the following Snakefile:
include: "pipeline.conf" # Should be an empty file
shell.prefix("set -euo pipefail;")
rule all:
input:
@marcelm
marcelm / info-file-with-qualities.patch
Created April 16, 2015 14:58
hack to print out qualities in cutadapt’s info file (will crash when trimming FASTA files)
diff --git i/cutadapt/scripts/cutadapt.py w/cutadapt/scripts/cutadapt.py
index 855721d..2eaf435 100755
--- i/cutadapt/scripts/cutadapt.py
+++ w/cutadapt/scripts/cutadapt.py
@@ -155,6 +155,7 @@ class AdapterCutter(object):
# TODO write only one line, even for multiple matches
for match in matches:
seq = match.read.sequence
+ qualities = match.read.qualities
if match is None:
@marcelm
marcelm / star-segfault.fasta
Created September 23, 2014 16:54
trying to map this to Danio rerio Zv9 with STAR leads to a segfault
>segfault
AACTGGAGCAGGTCCTGGTATTCTAACATTAATGTCTCCAGATGAACGCCCATAAACTTGGCTTTCACCTCAAACACTCC
CACAGTCTCTGATGGCGAGATCTCGAAAATTGCATTCTTATACTGATTTGTCTGCAGATCATCGATGTCAATAAGCACTC
CCTTCTCATGCAGACGGGCTGCTGTGTACTTCTGTGATACCTGCTTGCTCTTCTTTGCCTGTTTATCTCCAGGTTTCTTA
GAAACTTTTCCCTTGCTCCCTAAATTGTCCATGCAGGTCTTGATGTACTGATTGTAGTAGTCAATCTGCTCGTTGTAGAA
GGTTGCGGATGTCCTTGGCGATGTCATTGATTAGGTCCTGGTATTTCTTCTCAGGGTGTACTTTGCCCAACTCTCCCAGT
CTTTGCAGGTTGCTCTTGATCTTGTCCTTCTTGCCTTGCAAAGTAAGACTGTCATCTACAACTGGTTTGGCTTGCTTCAT
CTTCTCTGGGGTCTTTGCATCGCGAATGGCTCGCCGCTGCATAGCTCGCTGGTACTCTGCTTCCTGCTCAGGTGAAGCTA
CAGTCTCCAGAATCTCAGTAAGTGTGTCTCCTGGCTGAAATCTGATCACATCTACTATGAGCCTCTTGGTGTTAAGGAGA
AGAGTTTTGGCATCCATTTCAGCATTGGCCTCATCGGGAACATCAAATTTGTTGGTGAGGGTGAGGGAGACTTCTGTCTT