velvet怎么安裝使用

本篇內容主要講解“velvet怎么安裝使用”，感興趣的朋友不妨來看看。本文介紹的方法操作簡單快捷，實用性強。下面就讓小編來帶大家學習“velvet怎么安裝使用”吧!

velvet是由EMBL-EBI開發的一款基因組組裝工具

安裝過程如下

wget https://www.ebi.ac.uk/~zerbino/velvet/velvet_1.2.10.tgz
tar xzvf velvet_1.2.10.tgz
cd velvet_1.2.10/
make 'MAXKMERLENGTH=127'

默認情況下，velvet支持的kmer最大長度為31bp, 如果想要只會更大長度的kmer,在編譯時需要設置MAXKMERLENGTH的值。編譯完成后，會生成如下兩個可執行文件

1. velveth

2. velvetg

   
   

軟件的運行過程對應的也分成兩步

1. velveth

用法如下

velveth Assem 31 -shortPaired -fasta -separate left.fa right.fa

第一個參數Assem, 代表輸出結果的目錄；shortPaired指定測序類型，fasta指定輸入的序列格式。

對于二代測序平臺的數據，常用的測序類型包括以下兩種情況

1. short

2. shortPaired

   
   

short用于單獨數據，shortPaired用于雙端數據。輸入的序列文件支持以下格式：

1. fasta/fasta.gz

2. fastq/fastq.gz

3. sam/bam

   
   

通過不同的參數指定輸入文件的格式，-fasta對應fasta格式；-fastq對應fastq格式，-fastq.gz對應fastq.gz格式，-fasta.gz對應fasta.gz格式，-sam對應sam格式，-bam對應bam格式。

對于雙端數據，有以下兩種格式

1. interleaved

2. separate

   
   

R1和R2端序列保存在兩個文件中，就是separate格式；interleaved是雙端序列的一種格式，R1端和R2端的序列保存在一個文件當中，每一條序列的R1端之后緊跟著就是R2端序列；對于雙端測序，默認是interleaved格式，如果是separate格式，要顯示的聲明。

還需要注意的一個用法就是kmer長度，在實際分析時，通常會采用一系列的kmer長度分別組裝，然后挑選一個最佳的結果。velvet 的kmer參數可以設置為一個梯度，示例如下

velveth Assem 31,37,2 -shortPaired -fasta -separate left.fa right.fa

上述用法中的31,37,2表示從kmer=31開始組裝，然后進行遞增，步長為2，依次進行33, 35， 37共4個kmer長度的組裝。這樣的參數設計非常的貼合實際需求。

運行結束后，會在輸出目錄生成以下文件

1. Sequences

2. Roadmaps

2. velvetg

基本用法如下

velvetg Assembly/ -min_contig_lgth 100

第一個參數為上一步的輸出目錄，min_contig_lgth代表contig的最小長度，小于該長度的contig會被刪除，不會出現在最終的結果中。

運行結束后，輸出目錄下的contigs.fa就是最終的組裝結果。

到此，相信大家對“velvet怎么安裝使用”有了更深的了解，不妨來實際操作一番吧！這里是億速云網站，更多相關內容可以進入相關頻道進行查詢，關注我們，繼續學習！