# Change what's inside the quotes into your MBL STAMPS server username, or something else unique
# reminder! hit shift-ENTER to execute cell & move onto next cell
username="CHANGEME"

cd /mnt

mkdir -p $username

cd $username

# get a small E. coli data set
!curl -O http://public.ged.msu.edu.s3.amazonaws.com/ecoli-500k.fq.gz

!velveth k.31 31 -short -fastq.gz ecoli-500k.fq.gz

!velvetg k.31 -exp_cov auto

!head k.31/contigs.fa

import screed
contig_lengths = [ len(record.sequence) for record in screed.open('k.31/contigs.fa') ]

hist(contig_lengths, bins=20);