Velvet is a hash-based deBrujin assembler. I occasionally use it to build contigs from short read RNAseq data, although it can be used for long reads, metagenomics, etc. To install, just download, extract and make. There are two parts to the program: 1. velveth – prepare the sequences 2. velvetg – build the graph and...

GMAP is an old (circa 2006) software for long read alignment. Its use case is for mapping RNA reads back to a genome. It has found new life in the world of long read RNAseq such as from Pacbio reads. Perhaps because of its age and architecture, it has some quirks and dependencies that seem...

Its handy to have a webserver where you can put up BAM files for sharing with collaborators. But if you have tons of files or available webserver, you can use Amazon S3 to serve up your files. While S3 is pretty simple, setting it up for serving BAM (and index) files has a few tricky...

Packaging your data in R The goal of this is to make distribution of your data easy and consistent. The steps are not hard if you are using RStudio. The New Project wizard will do some of the work for you. Create a directory for the R project (myproject) Create a subdirectory (myproject/data-raw) In the...