This is an old revision of the document!
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.
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 .
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'.