LoRDEC 利用二代数据纠错PacBio 数据（ lordec-src

LoRDEC是一个利用二代测序短序列reads 矫正三代测序长reads 的工具，利用准确度高的NGS数据构建de Bruijn Graph（DBG）对三代长Reads进行纠错，适用于PacBio和Oxford 平台。

下载安装方法：
最新版本v0.9
（1）支持coda 安装，conda install -c atgc-montpellier lordec

(2) 可下载Binary package https://gite.lirmm.fr/lordec/lordec-releases/uploads/710113d83c210b6989ccfbdbafa89234/lordec-bin_0.9_linux64.tar.bz2
(3）可选择source 自己编译

wget https://gite.lirmm.fr/lordec/lordec-releases/uploads/800a96d81b3348e368a0ff3a260a88e1/lordec-src_0.9.tar.bz2
tar -jxvf lordec-src_0.9.tar.bz2
cd lordec-src_0.9
make clean all -j8

“all”—— 这个伪目标是所有目标的目标，其功能一般是编译所有的目标。

“clean” —— 这个伪目标功能是删除所有被make创建的文件。

软件需要 libc6 >=2.23 以上版本，
若低于此版本可升级新版本glibc

wget http://ftp.gnu.org/gnu/glibc/glibc-2.27.tar.gz
tar -xf glibc-2.27.tar.gz
cd glibc-2.27
mkdir buildcd build
../configure --prefix=/usr --disable-profile --enable-add-ons --with-headers=/usr/include --with-binutils=/usr/bin
make -j 10
make  install

软件使用说明

LoRDEC v0.9
using GATB v1.4.1
website : http://www.atgc-montpellier.fr/lordec/
FAQ : https://www.lirmm.fr/~rivals/lordec/FAQ/
./lordec-correct
-i|--long_reads <long read FASTA/Q file>
-2|--short_reads <short read FASTA/Q file(s)>
-k|--kmer_len <k-mer size>
-o|--corrected_read_file <output reads file>
-s|solid_threshold <solid k-mer abundance threshold>
[-t|--trials <number of paths to try from a k-mer>]
[-b|--branch <maximum number of branches to explore>]
[-e|--errorrate <maximum error rate>]
[-T|--threads <number of threads>]
[-S|--stat_file <out statistics file>]
[-c|--complete_search]
[-a|--abundance-max <abundance max threshold for k-mers>]
[-O|--out-tmp <GATB graph creation temporary files directory>]
[-p|--progress]
[-g|--graph_named_like_output]

../lordec-correct  -T 12  -2 ill-test-5K-1.fa,ill-test-5K-2.fa -k 19 -s 3  -a 10000 -i pacbio-reads.fa  -o corrected-pacbio-reads.fa &> mylog.log

-2 输入NGS 数据，支持fastq和fasta 格式，以及.gz的压缩格式
双端reads 以逗号隔开。如果输入多个fastq文件,则矫正时候会以所有的reads构建graph,然后矫正每条pacbio reads。
注意：如过是pair-end 序列，LoRDEC并不考虑配对信息，两个单端的library 无本质区别，且pair-end 的strand 方向对矫正无影响。
输入文件最好不要指定类似路径~/ ../ ，可能会出现以下错误

HDF5-DIAG: Error detected in HDF5 (1.8.18) thread 0
major: File accessibilty
minor: Unable to open file

另外，如果多个数据用到同一组NGS数据矫正，可以保存graph 文件为.h5格式，用于以后的矫正。

-a 矫正使用的最高coverage 可之指定高的数值，例如10000
-s 矫正使用的最低coverage 一般指定为3
-k kmer 数值越大，花费时间越多，可以测试多个数值，如19， 31, 41

LoRDEC 利用二代数据纠错PacBio 数据（ lordec-src_0.9.tar.bz2）相关推荐

三代测序数据纠错的方法、装置和计算机可读存储介质与流程
三代测序数据纠错的方法.装置和计算机可读存储介质与流程文档序号:15616049发布日期:2018-10-09 21:24 导航: X技术> 最新专利>计算;推算;计数设备的制造及其应用 ...
SPAdes混合组装二代、三代测序数据
导读 SPAdes是2012年发表在Journal of Computational Biology上的一篇文章提出的二代测序组装软件,是目前引用量已经达到6200+,在宏基因组组装软件中引用量最高[ ...
原核生物基因组三代数据（pacbio/nanopore）组装
动物随随便便几个G,植物动不动十多个G,而细菌基因组一般就几M,基因组结构简单,没有重复序列杂合度这些组装障碍.我现在很多文章细菌基因组开始用三代测序组装了,还有些文章是二代结合三代测序???纳税人的 ...
Pacbio 数据相关的几个重要概念
Pacbio 数据相关的几个重要概念 pacbio 基本概念 Polymerase Read, Subreads,Circular Consensus Sequence 的概念和相互关系 Pacbio ...
车载应用--AUXIN 利用 surfaceView 预览 Camera 数据
前文介绍了倒车后视,主要是利用 WindowManager 来加载布局,在布局中增加 SurfaceView 来预览 camera 数据.车载倒车后视–利用 WindowManager 预览 came ...
Case Study: 利用PHP获取关系型数据库中多张数据表的数据
一.目标该笔记的目的是引导读者借助WampServer平台和MySQL数据库,利用HTML/CSS/JS/PHP设计一个多数据表关联的网页.在上一个案例(Case Study: 利用JS实现数据库网 ...
Case Study: 利用JS实现数据库网页的数据分页、数据选择、数据详细信息查看功能
一.目标该笔记的目的是引导读者借助WampServer平台和MySQL数据库,利用HTML/CSS/JS/PHP设计一个能够进行实现数据分页显示.数据选择.数据详细信息查看功能的数据库网页.该数据库 ...
利用SQl对数据库实行数据拆分与组合
利用SQl对数据库实行数据拆分与组合实现提供以下几种方案: 方法一: WITH CTE AS ( SELECT A.Id,A.[Uid],UserName FROM ( SELECT A.[id], ...
【深度学习】利用一些API进行图像数据增广
[深度学习]利用一些API进行图像数据增广文章目录 [深度学习]利用一些API进行图像数据增广 1 先送上一份最强的翻转代码(基于PIL) 2 Keras中的数据增强API种类概述 3 特征标准化 ...

LoRDEC 利用二代数据纠错PacBio 数据（ lordec-src_0.9.tar.bz2）

LoRDEC 利用二代数据纠错PacBio 数据（ lordec-src_0.9.tar.bz2）相关推荐

最新文章

热门文章