Skip to content

Instantly share code, notes, and snippets.

View rpetit3's full-sized avatar
🦬
Yo!

Robert A. Petit III rpetit3

🦬
Yo!
View GitHub Profile
@rpetit3
rpetit3 / teton.config
Created June 16, 2022 17:24
nextflow config for teton
/*
This file includes default values for Teton.
*/
process {
executor = 'slurm'
queue = 'teton,moran'
scratch = params.slurm_use_scratch
time = 60.m
clusterOptions = '--account healthdatasci'
@rpetit3
rpetit3 / .bashrc
Created June 7, 2022 18:01
bash-full-of-colors with conda support
#!/usr/bin/env bash
# ~/.bashrc: executed by bash(1) for non-login shells.
# If not running interactively, don't do anything
[ -z "$PS1" ] && return
# don't put duplicate lines in the history. See bash(1) for more options
# ... or force ignoredups and ignorespace
HISTCONTROL=ignoredups:ignorespace
@rpetit3
rpetit3 / submit-cgc.py
Last active February 17, 2022 22:01
Example script for submitting and downloading results through CGC pipeline
#! /usr/bin/env python3
if __name__ == '__main__':
import os
import sys
import time
import argparse as ap
import sevenbridges as sbg
from sevenbridges.errors import SbgError
@rpetit3
rpetit3 / bam-to-cov.sh
Created January 7, 2022 16:28
per-base coverage from bam file
#! /bin/bash
GS_PATH=$1
# Copy Local
gsutil -q cp ${GS_PATH} ./
LOCAL_FILE=$(basename ${GS_PATH})
OUTDIR="./"
if echo "${LOCAL_FILE}" | grep ".rg.sorted.bam"; then
OUTDIR="./nanopore"
else
@rpetit3
rpetit3 / test-scrubber.sh
Created November 10, 2021 16:28
Quick script to compare two versions of sra-human-scrubber
#! /bin/bash
SCRUBBER_DB=${SCRUBBER_SHARE}/data/human_filter.db
# Current scrubber
echo "Run test on current Scrubber"
which scrub.sh
scrub.sh test
echo "Run test on PR Scrubber"
@rpetit3
rpetit3 / compare-coverages.sh
Last active July 23, 2021 18:48
Comparing Clearlabs reported coverage to observed coverage
#! /usr/bin/env bash
# https://www.ncbi.nlm.nih.gov/nuccore/1798174254
WUHAN=29903
while IFS=$'\t' read -r POSITION SAMPLE_ID ANALYSIS SEQUENCING_COVERAGE ASSEMBLY_COVERAGE; do
if [ "${POSITION}" == "POSITION" ]; then
printf "%s\t%s\t%s\t%s\t%s\n" "sample_id" "clearlabs_cov" "observed_cov" "mapped_coverage" "unmapped_coverage"
else
FASTQ=$(ls fastq/ | grep ${SAMPLE_ID})
if [ -f "fastq/${FASTQ}" ]; then
@rpetit3
rpetit3 / BaseSpaceRunDownloader_v2.py
Created June 9, 2021 15:55
Script to download projects from BaseSpace
from urllib2 import Request, urlopen, URLError
import json
import math
import sys
import os
import socket
import optparse
def arg_parser():
cwd_dir = os.getcwd()
@rpetit3
rpetit3 / HCMVL54.fasta
Last active March 2, 2021 14:59
Convert FASTA gene seqeunce to codon with position info
>X17403.1:c80631-76903 Human cytomegalovirus strain AD169 complete genome
ATGTTTTTCAACCCGTATCTGAGCGGCGGCGTGACCGGCGGTGCGGTCGCGGGTGGCCGGCGTCAGCGTT
CGCAGCCCGGCTCCGCGCAGGGCTCGGGCAAGCGGCCGCCACAGAAACAGTTTTTGCAGATCGTGCCGCG
AGGTGTCATGTTCGACGGTCAGACGGGGTTGATCAAGCATAAGACGGGACGGCTGCCTCTCATGTTCTAT
CGAGAGATTAAACATTTGTTGAGTCATGACATGGTTTGGCCGTGTCCTTGGCGCGAGACCCTGGTGGGTC
GCGTGGTGGGACCTATTCGTTTTCACACCTACGATCAGACGGACGCCGTGCTCTTCTTCGACTCGCCCGA
AAACGTGTCGCCGCGCTATCGTCAGCATCTGGTGCCTTCGGGGAACGTGTTGCGTTTCTTCGGGGCCACA
GAACACGGCTACAGTATCTGCGTCAACGTTTTCGGGCAGCGCAGCTACTTTTACTGTGAGTACAGCGACA
CCGATAGGCTGCGTGAGGTCATTGCCAGCGTGGGCGAACTAGTGCCCGAACCGCGGACGCCATACGCCGT
GTCTGTCACGCCGGCCACCAAGACCTCCATCTATGGGTACGGGACGCGACCCGTGCCCGATTTGCAGTGT
@rpetit3
rpetit3 / readme.md
Last active December 31, 2020 19:19
Setup NCOV

Started here: https://nextstrain.org/docs/getting-started/introduction

Installation

Installed on Chlamy (behind Emory VPN)

mkdir nextstrain
cd nextstrain
curl http://data.nextstrain.org/nextstrain.yml --compressed -o nextstrain.yml
conda env create -f nextstrain.yml
conda activate nextstrain
@rpetit3
rpetit3 / main.nf
Created November 15, 2020 06:24
bactopia-wrapper-nextflow
#! /usr/bin/env nextflow
import groovy.json.JsonSlurper
import groovy.text.SimpleTemplateEngine
import groovy.util.FileNameByRegexFinder
import java.nio.file.Path
import java.nio.file.Paths
import nextflow.util.SysHelper
PROGRAM_NAME = workflow.manifest.name
VERSION = workflow.manifest.version