Skip to content

Instantly share code, notes, and snippets.

View jts's full-sized avatar

Jared Simpson jts

View GitHub Profile
Slice 1/1, container offset 549, file offset 2589
Slice content type MAPPED_SLICE
Slice ref seq -1
Slice ref start 1
Slice ref span 1
Slice MD5 00000000000000000000000000000000
Rec counter 0
No. records 4311
No. blocks 26
Blk IDS: {11, 12, 25, 30, 32, 5779802, 6371138, 6377306, 6384474, 6382182, 6369894, 6376538, 6370921, 6369641, 6381926, 6370409, 6374746, 6370665, 6377050, 6386010, 6371418, 5779546, 6375258, 6370138, 6386778}
@jts
jts / Makefile
Last active September 25, 2018 19:16
nanopolish test makefile
#
# Sub directories containing source code, except for the main programs
SUBDIRS := src src/hmm src/thirdparty src/thirdparty/scrappie src/common src/alignment src/pore_model
#
# Set libraries, paths, flags and options
#
#Basic flags every build needs
@jts
jts / cnv_example
Last active August 29, 2015 14:02
cnv example
##fileformat=VCFv4.1
##fileDate=20140620
##reference=hs37d5
##INFO=<ID=CIEND,Number=2,Type=Integer,Description="Confidence interval around END for imprecise variants">
##INFO=<ID=CIPOS,Number=2,Type=Integer,Description="Confidence interval around POS for imprecise variants">
##INFO=<ID=END,Number=1,Type=Integer,Description="End position of the variant described in this record">
##INFO=<ID=IMPRECISE,Number=0,Type=Flag,Description="Imprecise structural variation">
##INFO=<ID=SVLEN,Number=.,Type=Integer,Description="Difference in length between REF and ALT alleles">
##INFO=<ID=SVTYPE,Number=1,Type=String,Description="Type of structural variant">
##ALT=<ID=DEL,Description="Deletion">
a45ff304-451c-4027-83fa-e37ba9c01bc1 133 19 2998502 0 * = 2998502 0 TGACATCATGCCTGGTGCACACCCACTCTGGTTCTGATACCGCCAGCTTCCAGAAAGGGCAGCACTTCAGTTCCATCTCTCTTGGTATTGAGGGAATGGAG
5437e56c-9c19-4b53-9568-b870629bbbd3 69 19 2998550 0 * = 2998550 0 GTTCTGATACCGCCAGCTTCCAGAAAGGGCAGCACTTCAGTTCCATCTCTCTTGGTATTGAGGGAATGGAGGCCCTGTTCTGCCCATTGGACTTAGCTCCA
8f570a1f-426c-48e1-afac-a866caa12613 69 19 2998551 0 * = 2998551 0 CATCATGCCTGGTGCACACCCACTCTGGTTCTGATACCGCCAGCTTCCAGAAAGGGCAGCACTTCAGTTCCATCTCTCTTGGTATTGAGGGAATGGAGGCC
18226cbc-505d-4eef-89a2-48575fd09b5e 133 19 2998698 0 * = 2998698 0 ACATCATGCCTGGTGCACACCCACTCTGGTTCTGATACCGCCAGCTTCCAGAAAGGGCAGCACTTCAGTTCCATCTCTCTTGGTATTGAGGGAATGGAGGC
19503915-e648-4036-98d1-2701784bc1d3 101 19 2998934 0 * = 2998934 0 CTGGTGCACACCCACTCTGGTTCTGATACCGCCAGCTTCCAGAAAGGGCAGCACTTCAGTTCCATCTCTCTTGGTATTGAGGGAATGGAGGCCCTGTTCTG
9548554d-cf08-4d29-883b-80d05f48c067 1
@jts
jts / jsbsub.pl
Created December 16, 2011 10:13
LSF submission wrapper
#! /usr/bin/perl
#
# Submission wrapper for LSF
# Recognizes parameters for the number
# of cpus to use, the amount of memory to reserve
# and the amount of disk space required.
#
use strict;
use Getopt::Long;
Getopt::Long::Configure('pass_through');
@jts
jts / gist:1119852
Created August 2, 2011 09:00
wavelet tree profiling
// Test case: ~600 Mbp of 150bp E. coli sequence reads
//
// SGA FM-index
//
Symbols: 605264456 Bytes: 121737401 (0.20 bytes per symbol)
Calculated 100000000 occurrence values in 21.28 seconds (0.212800 us/call)
//