User Tools

Site Tools


velvet

This is an old revision of the document!


Velvet

Velvet is a de novo genomic assembler specially designed for short read sequencing technologies, such as Solexa or 454, developed by Daniel Zerbino and Ewan Birney at the European Bioinformatics Institute (EMBL-EBI), near Cambridge, in the United Kingdom. You can find out more about Velvet here.

Environment Setup

To use Velvet first we will have to load the velvet module. You can do this by adding this line to your '.bashrc' file.

module load velvet

Before you run your code you'll want to create a directory to store the files as shown below.

razor-l1:jokinsey:~$ mkdir VELVET-TEST
razor-l1:jokinsey:~$ cd VELVET-TEST/
razor-l1:jokinsey:~/VELVET-TEST$ cp /share/apps/velvet/velvet_1.2.10/data/test_reads.fa .

Example Job

Now that we have our input file 'test_reads.fa', create a PBS script to run the job named 'velvet.pbs' that has the information below.

#!/bin/bash
#PBS -N velvet
#PBS -l nodes=1:ppn=12
#PBS -l walltime=00:05:00
#PBS -q tiny12core

cd $PBS_O_WORKDIR

velveth velvet1 29 test_reads.fa
velvetg velvet1

All that remains is to submit the job.

razor-l1:jokinsey:~/VELVET-TEST$ qsub velvet.pbs

The most important output files will be '/velvet1/configs.fa', '/velvet1/Graph', and '/velvet1/stats.txt'.

velvet.1507065248.txt.gz · Last modified: 2017/10/03 21:14 by jokinsey