Difference between revisions of "Velvet"

From UFRC
Jump to navigation Jump to search
 
(54 intermediate revisions by 5 users not shown)
Line 1: Line 1:
= Velvet<br>  =
+
__NOTOC__
 
+
[[Category:Software]][[Category:Biology]][[Category:Genomics]]
'''Sequence assembler for very short reads.''' 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.
+
{|<!--CONFIGURATION: REQUIRED-->
 
+
|{{#vardefine:app|velvet}}
<br>
+
|{{#vardefine:url|http://www.ebi.ac.uk/~zerbino/velvet/}}
 
+
<!--CONFIGURATION: OPTIONAL (1=ON)-->
Velvet currently takes in short read sequences, removes errors then produces high quality unique contigs. It then uses paired-end read and long read information, when available, to retrieve the repeated areas between contigs. <br>
+
|{{#vardefine:exe|1}} <!--RUNNING -->
 
+
|{{#vardefine:conf|}} <!--CONFIGS-->
<br>
+
|{{#vardefine:pbs|}} <!--PBS SCRIPTS-->
 
+
|{{#vardefine:policy|}} <!--POLICY-->
== Velvet has been installed on the cluster in several versions:<br>  ==
+
|{{#vardefine:testing|}} <!--TESTING/PROFILING-->
<blockquote><br></blockquote>  
+
|{{#vardefine:faq|}} <!--FAQ-->
{| cellspacing="1" cellpadding="1" border="1" style="width: 954px; height: 138px;"
+
|{{#vardefine:citation|}} <!--CITATION-->
|-
+
|{{#vardefine:installation|}} <!--INSTALLATION-->
| /apps/velvet/1.1.05/velveth<br>/apps/velvet/1.1.05/velvetg<br>
 
| These versions are the defailt compile options&nbsp;(max kmer=31, single threaded)<br>
 
|-
 
| /apps/velvet/1.1.05/velveth_OMP<br>/apps/velvet/1.1.05/velvetg_OMP<br>  
 
| These versions are compiled with the default kmer, but using OpenMP--see OpenMP notes below for details<br>
 
|-
 
| /apps/velvet/1.1.05/velveth_max99_OMP<br>/apps/velvet/1.1.05/velvetg_max99_OMP<br>  
 
| These versions are compiled with a max kmer of 99 and OpenMP--see OpenMP notes below for details<br>
 
 
|}
 
|}
<blockquote>Note that all versions were compiled with icc. <br> Please file a bugzilla request for additional compile options like color space, different kmers, LONGSEQUENCES, etc.<br> </blockquote>  
+
<!--BODY-->
<br>  
+
<!--Description-->
 
+
{{#if: {{#var: url}}|
== OpenMP Notes<br>  ==
+
{{App_Description|app={{#var:app}}|url={{#var:url}}|name={{#var:app}}}}|}}
  
<br>
+
Sequence assembler for very short reads.
  
If you use one of the OpenMP versions of velvet, you must set the enviromental variable OMP_THREAD_LIMIT to the same value as your processor request in your submission script. For example if you use #PBS -l nodes=1:ppn=8 in your script also include export OMP_THREAD_LIMIT=8 (for a bash script).<br>
+
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.  
  
<br>
+
Velvet currently takes in short read sequences, removes errors then produces high quality unique contigs. It then uses paired-end read and long read information, when available, to retrieve the repeated areas between contigs.
  
== Website<br>  ==
+
In the Velvet module version 1.2.06 [[MetaVelvet|MetaVelvet-1.1.01]] was made available. There are two versions compiled, meta-velvetg (using a max kmer of 31) and meta-velvetg_MAX99 (using a max kmer of 99). There are also multiple contributed scripts available in the same location as the velvet binaries.
  
[http://www.ebi.ac.uk/~zerbino/velvet/ Velvet home page].<br>
+
However, in the Velvet module version 1.2.07 only the contributed scripts are present. [[MetaVelvet|MetaVelvet-1.2.01]] is now provided by its own "metavelvet" module.
  
<br>  
+
Velvet memory requirements increase dramatically with each kmer length increase, so the default velveth and velvetg binaries only handle kmers up to 31 bp and we provide additional binaries to handle kmers up to 99, and 249 bp. The binaries are named with _maxNN suffixes e.g. velveth_max249. In addition _OMP suffixes are used to identify the binaries built for multi-threaded usage (see below).
 +
 +
<!--Modules-->
 +
==Environment Modules==
 +
Run <code>module spider {{#var:app}}</code> to find out what environment modules are available for this application.
 +
==System Variables==
 +
* HPC_{{uc:{{#var:app}}}}_DIR - installation directory
 +
* HPC_VELVET_BIN - executable directory
 +
<!--Additional-->
 +
{{#if: {{#var: exe}}|==Additional Information==
 +
# The user manual is available [http://www.google.com/url?sa=t&source=web&cd=2&sqi=2&ved=0CBkQFjAB&url=http%3A%2F%2Fhelix.nih.gov%2FApplications%2Fvelvet_manual.pdf&rct=j&q=velvet%20assembler%20manual&ei=Lr1WToqdIsGG0QGJyamyDA&usg=AFQjCNH2VB2-la3kWazMqdl3ms-xkCocVw here.]
  
== Example submission script<br>  ==
+
# If you use one of the OpenMP versions of velvet, you must set the enviroment variables '''OMP_NUM_THREADS''' and '''OMP_THREAD_LIMIT''' to the appropriate value based on the processor request in your submission script. For example if you use #PBS -l nodes=1:ppn=8 in your script, also include export OMP_NUM_THREADS=7; export OMP_THREAD_LIMIT=8 (for a bash script) in the script itself (see sample script below and/or the [http://www.google.com/url?sa=t&source=web&cd=2&sqi=2&ved=0CBkQFjAB&url=http%3A%2F%2Fhelix.nih.gov%2FApplications%2Fvelvet_manual.pdf&rct=j&q=velvet%20assembler%20manual&ei=Lr1WToqdIsGG0QGJyamyDA&usg=AFQjCNH2VB2-la3kWazMqdl3ms-xkCocVw user manual] for more information).
<blockquote>!/bin/bash<br>#PBS -N velvet<br>#PBS -o velvet.out<br>#PBS -e velvet.err<br>#PBS -M &lt;your e-mail addres&gt;<br>#PBS -m abe<br>#PBS -l walltime=12:00:00&nbsp;&nbsp;&nbsp;&nbsp; #&nbsp;HH:MM:SS<br>#PBS -l pmem=4gb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # 4GB&nbsp;of RAM&nbsp;per thread, or 4X8=32GB&nbsp;in this case<br>#PBS -l nodes=1:ppn=8 </blockquote><blockquote>cd $PBS_O_WORKDI</blockquote><blockquote>#Set OMP_THREAD_LIMIT--should be the same as ppn above
+
# Please file a [https://bugzilla.rc.ufl.edu/ bugzilla] request if you need an executable with different compile-time options such as color space, different kmers, LONGSEQUENCES, etc.
export OMP_THREAD_LIMIT=8  
+
}}
/apps/velvet/1.1.05/velveth_max99_OMP velvet_out/ 45 -fastq.gz -shortPaired my_paried_data.fastq.gz<br>/apps/velvet/1.1.05/velvetg_max99_OMP velvet_out/ -min_contig_lgth 250 -exp_cov 10 -ins_length 350<br></blockquote>
+
{{#if: {{#var: pbs}}|==PBS Script Examples==
<br>
+
See the [[{{PAGENAME}}_PBS]] page for {{#var: app}} PBS script examples.|}}

Latest revision as of 17:06, 14 December 2022

Description

velvet website  

Sequence assembler for very short reads.

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.

Velvet currently takes in short read sequences, removes errors then produces high quality unique contigs. It then uses paired-end read and long read information, when available, to retrieve the repeated areas between contigs.

In the Velvet module version 1.2.06 MetaVelvet-1.1.01 was made available. There are two versions compiled, meta-velvetg (using a max kmer of 31) and meta-velvetg_MAX99 (using a max kmer of 99). There are also multiple contributed scripts available in the same location as the velvet binaries.

However, in the Velvet module version 1.2.07 only the contributed scripts are present. MetaVelvet-1.2.01 is now provided by its own "metavelvet" module.

Velvet memory requirements increase dramatically with each kmer length increase, so the default velveth and velvetg binaries only handle kmers up to 31 bp and we provide additional binaries to handle kmers up to 99, and 249 bp. The binaries are named with _maxNN suffixes e.g. velveth_max249. In addition _OMP suffixes are used to identify the binaries built for multi-threaded usage (see below).

Environment Modules

Run module spider velvet to find out what environment modules are available for this application.

System Variables

  • HPC_VELVET_DIR - installation directory
  • HPC_VELVET_BIN - executable directory

Additional Information

  1. The user manual is available here.
  1. If you use one of the OpenMP versions of velvet, you must set the enviroment variables OMP_NUM_THREADS and OMP_THREAD_LIMIT to the appropriate value based on the processor request in your submission script. For example if you use #PBS -l nodes=1:ppn=8 in your script, also include export OMP_NUM_THREADS=7; export OMP_THREAD_LIMIT=8 (for a bash script) in the script itself (see sample script below and/or the user manual for more information).
  2. Please file a bugzilla request if you need an executable with different compile-time options such as color space, different kmers, LONGSEQUENCES, etc.