源代碼編譯安裝velvets

apt-get默認(rèn)安裝的velvet最大km值只有31；源代碼編譯安裝的話可以修改這個(gè)值到127

# 源碼編譯安裝velvet
# cd ~/Biosofts/
# wget https://www.ebi.ac.uk/~zerbino/velvet/velvet_1.2.10.tgz
tar zvxf /disk/shares/velvet_1.2.10.tgz -C ~/Biosofts/
cd ~/Biosofts/velvet_1.2.10/
make 'CATEGORIES=10' 'MAXKMERLENGTH=127' 'LONGSEQUENCES=1'  'OPENMP=1' 'BUNDLEDZLIB=1'
 ./velveth -h

最后一步./velveth -h返回截圖如上

比較不同kmer值對(duì)velvet組裝結(jié)果的影響

nano kmerselection.sh
#輸入如下代碼
################################
#!/bin/bash

Data1=/disk/shares/Seqs/test_7942raw_1.fq.gz
Data2=/disk/shares/Seqs/test_7942raw_2.fq.gz
Velvet_dir=~/Biosofts/velvet_1.2.10

for k in `seq 23 8 127`;do
    mkdir velvet$k
    $Velvet_dir/velveth velvet$k $k -shortPaired -fastq -separate $Data1 $Data2
    $Velvet_dir/velvetg velvet$k -exp_cov auto -cov_cutoff auto -very_clean yes

 done
############################
./kmerselection.sh

最后一步運(yùn)行結(jié)果

Quast安裝

#預(yù)編譯安裝
tar zvxf /disk/shares/quast-5.0.0.tar.gz -C ~/Biosofts/
#進(jìn)入軟件目錄
cd Biosofts/quast-5.0.0/
#測試是否可以運(yùn)行
./quast.py
#添加環(huán)境變量
echo 'export PATH=~/Biosofts/quast-5.0.0:$PATH' >> ~/.bashrc
source ~/.bashrc

比較不同km值拼接效果

quast.py -o 
compare_velvet_kmer 
velvet23/contigs.fa 
velvet31/contigs.fa 
velvet39/contigs.fa 
velvet47/contigs.fa
 velvet55/contigs.fa
 velvet63/contigs.fa
 velvet71/contigs.fa 
velvet79/contigs.fa 
velvet87/contigs.fa
 velvet95/contigs.fa 
velvet103/contigs.fa
 velvet111/contigs.fa 
velvet119/contigs.fa 
velvet127/contigs.fa

Quast截圖

Quast結(jié)果

Quast結(jié)果