Skip to content

Instantly share code, notes, and snippets.

Greg Caporaso gregcaporaso

Block or report user

Report or block gregcaporaso

Hide content and notifications from this user.

Learn more about blocking users

Contact Support about this user’s behavior.

Learn more about reporting abuse

Report abuse
View GitHub Profile
@gregcaporaso
gregcaporaso / README.md
Created Nov 17, 2012
quick and dirty script to create a barcode read fastq file from a sequence read fastq file with barcodes in the headers
View README.md

USAGE: extract_fastq_barcodes_from_header.py input_reads.fastq barcode_reads.fastq

@gregcaporaso
gregcaporaso / notes.md
Last active Mar 28, 2019
vsearch-based sequence dereplication through generation of a biom table
View notes.md

This depends on biom version >= 2.1.5, < 2.2.0 and vsearch >= 1.7.0.

Please note that all of this is highly experimental. I'm keeping these notes as I work with this approach. For published work, I still recommend using standard pipelines, such as those in QIIME 1.9.1.

$ biom --version
biom, version 2.1.5

$ vsearch --version
vsearch v1.7.0_osx_x86_64, 16.0GB RAM, 8 cores
@gregcaporaso
gregcaporaso / neighbor-joining-experiments.ipynb
Created Jun 11, 2014
initial experiments with implementing neighbor joining for scikit-bio
View neighbor-joining-experiments.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@gregcaporaso
gregcaporaso / qiime2-issue190.ipynb
Created Dec 1, 2016
QIIME 2 issue #190 example
View qiime2-issue190.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View jgcaporaso-qiime2-scipy2016.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@gregcaporaso
gregcaporaso / qiime-exercise.ipynb
Created Apr 5, 2016
Final assignment for BIO/CS 499/599 (Spring 2016)
View qiime-exercise.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View q2d3-test.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View artifact-demo.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@gregcaporaso
gregcaporaso / check_illumina_barcodes.py
Created Dec 3, 2013
Script for performing some basic testing on Illumina amplicon sequencing primers as described in: http://www.nature.com/ismej/journal/v6/n8/full/ismej20128a.html
View check_illumina_barcodes.py
#!/usr/bin/env python
# File created on 01 Dec 2011
from __future__ import division
__author__ = "Greg Caporaso"
__copyright__ = "Copyright 2011, The QIIME project"
__credits__ = ["Greg Caporaso"]
__license__ = "GPL"
__version__ = "1.3.0-dev"
__maintainer__ = "Greg Caporaso"
View fq_qual.py
#!/usr/bin/env python
from sys import argv
from random import random
from cogent.parse.fastq import MinimalFastqParser
from cogent.draw.distribution_plots import generate_box_plots
from qiime.quality import ascii_to_phred33
from qiime.util import qiime_open
def fastq_quality_plots(fastq_records,
You can’t perform that action at this time.