Skip to content

Instantly share code, notes, and snippets.

@jaudoux
jaudoux / Snakefile
Created April 20, 2018 09:46
Generate k-mer matrice from FASTQ files
RAW_DIR = "raw"
SAMPLES = ['SRR2966453', 'SRR2966456','SRR2966471','SRR2966474', 'SRR2966454', 'SRR2966457', 'SRR2966472', 'SRR2966475', 'SRR2966455', 'SRR2966458', 'SRR2966473', 'SRR2966476']
rule all:
input: "raw_counts.tsv.gz.gbi"
rule jellyfish_counts:
input:
r1 = RAW_DIR + "/{sample}_1.fastq.gz",
@jaudoux
jaudoux / sampleTosampleJaccardIndex.pl
Created April 9, 2018 07:37
Comptute jaccard index from a k-mer abundance matrix
#! /usr/bin/perl
#
use strict;
use warnings;
my %jaccard_h;
while(<>) {
chomp;
my($k,@v) = split "\t", $_;
for(my $i = 0; $i < scalar @v; $i++) {
@jaudoux
jaudoux / gitio-test
Created June 30, 2014 10:26
Testing gist.io
# Testing gist.io
is-it cool?