软件介绍
Flye 是一种用于单分子测序读数（如 PacBio 和牛津纳米孔技术公司生产的读数）的从头组装器。它适用于各种数据集，从小型细菌项目到大型哺乳动物级组装。该软件包代表了一个完整的流水线：它将原始的 PacBio / ONT 读取数据作为输入，并输出精加工的等位基因。
软件的安装
1.直接conda安装
conda install flye
2.编译安装
在github上获得最新的软件包，进行编译
git clone https://github.com/fenderglass/Flye
cd Flye
make
3.python直接编译安装，命令直接进入环境变量
git clone https://github.com/fenderglass/Flye
cd Flye
python setup.py install
软件的使用
flye将原始的 PacBio / ONT 数据作为输入
1.先查看软件的帮助信息
$ flye -h
usage: flye (--pacbio-raw | --pacbio-corr | --pacbio-hifi | --nano-raw |
--nano-corr | --nano-hq ) file1 [file_2 ...]
--out-dir PATH
[--genome-size SIZE] [--threads int] [--iterations int]
[--meta] [--polish-target] [--min-overlap SIZE]
[--keep-haplotypes] [--debug] [--version] [--help]
[--scaffold] [--resume] [--resume-from] [--stop-after]
[--read-error float] [--extra-params]
[--deterministic]
Assembly of long reads with repeat graphs
options:
-h, --help show this help message and exit
--pacbio-raw path [path ...]
PacBio regular CLR reads (<20% error)
--pacbio-corr path [path ...]
PacBio reads that were corrected with other methods (<3% error)
--pacbio-hifi path [path ...]
PacBio HiFi reads (<1% error)
--nano-raw path [path ...]
ONT regular reads, pre-Guppy5 (<20% error)
--nano-corr path [path ...]
ONT reads that were corrected with other methods (<3% error)
--nano-hq path [path ...]
ONT high-quality reads: Guppy5+ SUP or Q20 (<5% error)
--subassemblies path [path ...]
[deprecated] high-quality contigs input
-g size, --genome-size size
estimated genome size (for example, 5m or 2.6g)
-o path, --out-dir path
Output directory
-t int, --threads int
number of parallel threads [1]
-i int, --iterations int
number of polishing iterations [1]
-m int, --min-overlap int
minimum overlap between reads [auto]
--asm-coverage int reduced coverage for initial disjointig assembly [not set]
--hifi-error float [deprecated] same as --read-error
--read-error float adjust parameters for given read error rate (as fraction e.g. 0.03)
--extra-params extra_params
extra configuration parameters list (comma-separated)
--plasmids unused (retained for backward compatibility)
--meta metagenome / uneven coverage mode
--keep-haplotypes do not collapse alternative haplotypes
--no-alt-contigs do not output contigs representing alternative haplotypes
--scaffold enable scaffolding using graph [disabled by default]
--trestle [deprecated] enable Trestle [disabled by default]
--polish-target path run polisher on the target sequence
--resume resume from the last completed stage
--resume-from stage_name
resume from a custom stage
--stop-after stage_name
stop after the specified stage completed
--debug enable debug output
-v, --version show program's version number and exit
--deterministic perform disjointig assembly single-threaded
Input reads can be in FASTA or FASTQ format, uncompressed
or compressed with gz. Currently, PacBio (CLR, HiFi, corrected)
and ONT reads (regular, HQ, corrected) are supported. Expected error rates are
<15% for PB CLR/regular ONT; <5% for ONT HQ, <3% for corrected, and <1% for HiFi. Note that Flye
was primarily developed to run on uncorrected reads. You may specify multiple
files with reads (separated by spaces). Mixing different read
types is not yet supported. The --meta option enables the mode
for metagenome/uneven coverage assembly.
To reduce memory consumption for large genome assemblies,
you can use a subset of the longest reads for initial disjointig
assembly by specifying --asm-coverage and --genome-size options. Typically,
40x coverage is enough to produce good disjointigs.
You can run Flye polisher as a standalone tool using
2.通过查看参数，对应我们的原始数据，选择正确的参数进行运行
--pacbio-raw | --pacbio-corr | --pacbio-hifi | --nano-raw |--nano-corr | --nano-hq
3.看看实例
flye --pacbio-hifi hifi_reads.fastq.gz --out-dir ./03.flye --iterations 3 --genome-size 1.4g --threads 40
#参数介绍
--pacbio-hifi #我们的是PacBio HiFi数据，所以选择这个参数
--out-dir #输出文件的路劲
--iterations 3 #进行3次抛光（polish）
--genome-size 1.4g #预估的基因组大小（单位m,g）
--threads 40 #调取的线程数
