velvet怎么安裝使用

本篇內(nèi)容主要講解“velvet怎么安裝使用”,感興趣的朋友不妨來看看。本文介紹的方法操作簡單快捷,實用性強(qiáng)。下面就讓小編來帶大家學(xué)習(xí)“velvet怎么安裝使用”吧!

成都創(chuàng)新互聯(lián)于2013年成立,先為廣西等服務(wù)建站,廣西等地企業(yè),進(jìn)行企業(yè)商務(wù)咨詢服務(wù)。為廣西企業(yè)網(wǎng)站制作PC+手機(jī)+微官網(wǎng)三網(wǎng)同步一站式服務(wù)解決您的所有建站問題。

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

安裝過程如下

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'

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

  1. velveth

  2. velvetg

軟件的運行過程對應(yīng)的也分成兩步

1. velveth

用法如下

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

第一個參數(shù)Assem, 代表輸出結(jié)果的目錄;shortPaired指定測序類型,fasta指定輸入的序列格式。

對于二代測序平臺的數(shù)據(jù),常用的測序類型包括以下兩種情況

  1. short

  2. shortPaired

short用于單獨數(shù)據(jù),shortPaired用于雙端數(shù)據(jù)。輸入的序列文件支持以下格式:

  1. fasta/fasta.gz

  2. fastq/fastq.gz

  3. sam/bam

通過不同的參數(shù)指定輸入文件的格式,-fasta對應(yīng)fasta格式;-fastq對應(yīng)fastq格式,-fastq.gz對應(yīng)fastq.gz格式,-fasta.gz對應(yīng)fasta.gz格式,-sam對應(yīng)sam格式,-bam對應(yīng)bam格式。

對于雙端數(shù)據(jù),有以下兩種格式

  1. interleaved

  2. separate

R1和R2端序列保存在兩個文件中,就是separate格式;interleaved是雙端序列的一種格式,R1端和R2端的序列保存在一個文件當(dāng)中,每一條序列的R1端之后緊跟著就是R2端序列;對于雙端測序,默認(rèn)是interleaved格式,如果是separate格式,要顯示的聲明。

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

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

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

運行結(jié)束后,會在輸出目錄生成以下文件

  1. Sequences

  2. Roadmaps

2. velvetg

基本用法如下

velvetg Assembly/ -min_contig_lgth 100

第一個參數(shù)為上一步的輸出目錄,min_contig_lgth代表contig的最小長度,小于該長度的contig會被刪除,不會出現(xiàn)在最終的結(jié)果中。

運行結(jié)束后,輸出目錄下的contigs.fa就是最終的組裝結(jié)果。

到此,相信大家對“velvet怎么安裝使用”有了更深的了解,不妨來實際操作一番吧!這里是創(chuàng)新互聯(lián)網(wǎng)站,更多相關(guān)內(nèi)容可以進(jìn)入相關(guān)頻道進(jìn)行查詢,關(guān)注我們,繼續(xù)學(xué)習(xí)!

網(wǎng)站標(biāo)題:velvet怎么安裝使用
文章URL:http://muchs.cn/article16/gdipgg.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供品牌網(wǎng)站建設(shè)、外貿(mào)建站、網(wǎng)站導(dǎo)航、網(wǎng)站維護(hù)、云服務(wù)器Google

廣告

聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請盡快告知,我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時需注明來源: 創(chuàng)新互聯(lián)

h5響應(yīng)式網(wǎng)站建設(shè)