This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/usr/bin/env perl | |
use strict; | |
$/ = ">"; | |
<STDIN>; # Discard "junk", if any, at beginning of the file. | |
my @sequencelengths; | |
my $gccontent = 0; | |
my $atcontent = 0; | |
my $combinedlength = 0; |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
# -q: output in Fastq format | |
# -Q: ignore BAM quality flags | |
# -P: paired-end data | |
bam2fastx -qQP -o clean.fq <(bwa mem contaminants.fasta reads.1.fq reads.2.fq | \ | |
perl -ne '@v = split(/\t/); print if(m/^@/ or ($v[1] & 4 and $v[1] & 8))' | \ | |
samtools view -bhS -) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
ProbandSex | Chrom | MotherDist | FatherDist | ProbandDist | |
---|---|---|---|---|---|
Male | X | mu | mu/2 | mu/2 | |
Male | Y | ??? | mu/2 | mu/2 | |
Female | X | mu | mu/2 | mu |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Proband | Mother | Father | ValidAutosomal | ValidSonChrX | ValidSonChrY | ValidDaughterChrX | |
---|---|---|---|---|---|---|---|
0/0 | 0/0 | 0/0 | + | + | + | + | |
0/0 | 0/0 | 0/1 | + | + | + | + | |
0/0 | 0/0 | 1/1 | - | - | - | - | |
0/0 | 0/1 | 0/0 | + | + | - | + | |
0/0 | 0/1 | 0/1 | + | + | - | + | |
0/0 | 0/1 | 1/1 | - | - | - | - | |
0/0 | 1/1 | 0/0 | - | - | - | - | |
0/0 | 1/1 | 0/1 | - | - | - | - | |
0/0 | 1/1 | 1/1 | - | - | - | - |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
alias list='ls -lhp' | |
alias vsc='open -a /Applications/Visual\ Studio\ Code.app/' | |
alias msu='ssh hpcc.msu.edu' | |
alias cab='ssh cabernet.genomecenter.ucdavis.edu' | |
alias rum='ssh rum.genomecenter.ucdavis.edu' | |
alias gremlin='ssh gremlin2.soic.indiana.edu' | |
export CLICOLOR=1 | |
export HISTSIZE=1000000 | |
export HISTFILESIZE=1000000 | |
export EDITOR=nano |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/usr/bin/env python | |
import argparse | |
import kevlar | |
import random | |
import sys | |
def resv_samp(objstream, n=100, stopafter=None, filterfunc=None): | |
sample = list() | |
for counter, obj in enumerate(objstream): | |
if filterfunc: |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
INDIVIDUALS = ('proband', 'mother', 'father') | |
KSIZES = ('27', '31', '35', '39', '43', '47', '51') | |
rule all: | |
input: | |
expand('{coverage}x_k{ksize}_kevlar_calls_like.vcf', | |
coverage=('30'), ksize=KSIZES) | |
rule novel_reads: |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/usr/bin/env perl | |
# Copyright (c) 2012-2013, Daniel S. Standage <daniel.standage@gmail.com> | |
# | |
# Permission to use, copy, modify, and/or distribute this software for any | |
# purpose with or without fee is hereby granted, provided that the above | |
# copyright notice and this permission notice appear in all copies. | |
# | |
# THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES | |
# WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/usr/bin/env python | |
import argparse | |
from collections import defaultdict | |
import khmer | |
from math import log | |
import re | |
import scipy.stats | |
import sys |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#include <iostream> | |
#include <assert.h> | |
#include <cmath> | |
#include "filter.hpp" | |
template<typename ElementType, typename CounterType, size_t maxcount> | |
filter<ElementType, CounterType, maxcount>::filter(std::vector<size_t> array_sizes) | |
: _cells_occupied(array_sizes.size(), 0), | |
_arrays(array_sizes.size(), std::vector<CounterType>()) | |
{ |
NewerOlder