二代组装

Gapcloser

2019-03-28  本文已影响0人  nitrostarch

下载和安装
https://sourceforge.net/projects/soapdenovo2/files/GapCloser/

使用

GapCloser -a scaffolds.fasta -b input_library.txt -o gapcloser_scaffolds.fasta -l 149 -t 4

需要编辑input文件

#maximal read length
max_rd_len=149
[LIB]
#average insert size
avg_ins=408
#if sequence needs to be reversed
reverse_seq=0
#in which part(s) the reads are used
asm_flags=4
#use only first 100 bps of each read
rd_len_cutoff=100
#in which order the reads are used while scaffolding
rank=1
# cutoff of pair number for a reliable connection (at least 3 for short insert size)
pair_num_cutoff=3
#minimum aligned length to contigs for a reliable read location (at least 32 for short insert size)
map_len=32
#a pair of fastq file, read 1 file should always be followed by read 2 file
q1=/mnt/e/linux/experiment/data/fastp/corrected_cutf_uni_10m_1.fastq
q2=/mnt/e/linux/experiment/data/fastp/corrected_cutf_uni_10m_2.fastq
上一篇下一篇

猜你喜欢

热点阅读