Skip to content

Instantly share code, notes, and snippets.

seqid Name
5 Transcript_10 EDO31303
9 Transcript_100 EDO42116
19 Transcript_101 EDO42087
24 Transcript_102 EDO43514
39 Transcript_103 tr|A7SUP6|A7SUP6_NEMVE Predicted protein OS=Nematostella vectensis OX=45351 GN=v1g217764 PE=4 SV=1
58 Transcript_105 EDO33432
60 Transcript_107 tr|A7S0C0|A7S0C0_NEMVE Predicted protein OS=Nematostella vectensis OX=45351 GN=v1g184098 PE=4 SV=1
61 Transcript_109 EDO49229
71 Transcript_11 tr|A7RQ19|A7RQ19_NEMVE Predicted protein OS=Nematostella vectensis OX=45351 GN=v1g239897 PE=4 SV=1
>hairpin_1_fwd
GATCGGAAGAGCACACGTCTGAACTCCAGTC
>hairpin_1_revcomp
GACTGGAGTTCAGACGTGTGCTCTTCCGATC
>hairpin_2_fwd
ACACTCTTTCCCTACACGACGCTCTTCCGATC
>hairpin_2_revcomp
GATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT
>Illumina_top
ACACTCTTTCCCTACACGACGCTCTTCCGATC
##########################################
## Intersect and Venn Diagram Functions ##
##########################################
## Author: Thomas Girke
## Last update: August 15, 2009
## Utilities:
## (1) Venn Intersects
## Computation of Venn intersects among 2-20 or more sample sets using the typical
## 'only in' intersect logic of Venn comparisons, such as: objects present only in
## set A, objects present only in the intersect of A & B, etc. Due to this restrictive
# custom script by Igor Dolgolev at NYUMC
plotPCAWithSampleNames = function(x, intgroup="condition", ntop=500)
{
library("genefilter")
library("lattice")
rv = rowVars(assay(x))
select = order(rv, decreasing=TRUE)[seq_len(min(ntop, length(rv)))]
###################################################################
"""Function: Take in a dammit gff3, dammit fasta, dammit namemap.
Output gff3 and fasta with trinity names.
by Dr. Tessa Pierce
"""
###################################################################
import sys
import os
import argparse
#outputs top 20 barcodes in file, if there are any
zcat lane1_sample_L001_R1_001.fastq.gz | head -400000 | grep '^@HWI' | cut -d : -f 10 | sort | uniq -c | sort -nr | head -20
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
>DNA_CS
GCCATCAGATTGTGTTTGTTAGTCGCTTTTTTTTTTTGGAATTTTTTTTTTGGAATTTTTTTTTTGCGCTAACAACCTCCTGCCGTTTTGCCCGTGCATATCGGTCACGAACAAATCTGATTACTAAACACAGTAGCCTGGATTTGTTCTATCAGTAATCGACCTTATTCCTAATTAAATAGAGCAAATCCCCTTATTGGGGGTAAGACATGAAGATGCCAGAAAAACATGACCTGTTGGCCGCCATTCTCGCGGCAAAGGAACAAGGCATCGGGGCAATCCTTGCGTTTGCAATGGCGTACCTTCGCGGCAGATATAATGGCGGTGCGTTTACAAAAACAGTAATCGACGCAACGATGTGCGCCATTATCGCCTAGTTCATTCGTGACCTTCTCGACTTCGCCGGACTAAGTAGCAATCTCGCTTATATAACGAGCGTGTTTATCGGCTACATCGGTACTGACTCGATTGGTTCGCTTATCAAACGCTTCGCTGCTAAAAAAGCCGGAGTAGAAGATGGTAGAAATCAATAATCAACGTAAGGCGTTCCTCGATATGCTGGCGTGGTCGGAGGGAACTGATAACGGACGTCAGAAAACCAGAAATCATGGTTATGACGTCATTGTAGGCGGAGAGCTATTTACTGATTACTCCGATCACCCTCGCAAACTTGTCACGCTAAACCCAAAACTCAAATCAACAGGCGCCGGACGCTACCAGCTTCTTTCCCGTTGGTGGGATGCCTACCGCAAGCAGCTTGGCCTGAAAGACTTCTCTCCGAAAAGTCAGGACGCTGTGGCATTGCAGCAGATTAAGGAGCGTGGCGCTTTACCTATGATTGATCGTGGTGATATCCGTCAGGCAATCGACCGTTGCAGCAATATCTGGGCTTCACTGCCGGGCGCTGGTTATGGTCAGTTCGAGCATAAGGCTGACAGCCTGATTGCAAAATTCAAAGAAGCGGGCGGAACGGTCAGAGAGATTGATGTATGAGCAGAGTCACCGCGATTATCTCC
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.