加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

使用SPAdes软件对基因组进行拼接

(2015-08-04 09:56:21)
标签:

it

分类: 生物信息学
1:网址:http://bioinf.spbau.ru/spades#benchmark

2:可以拼接来自多个平台的测序数据,如果只有一个库,可以利用来自不同的软件的组装结果(参见软件参数:--untrusted-contigs或--trusted-contigs  ),除了不兼容Illumina与IonTorrent文库的混合文库拼接外,其他的平台数据不同文库都可以实现混合拼接。

3:该软件包含了前期对测序数据的错误纠正(

Quake\

hammer)以及基于kmer的错配纠正( mismatch correction

4:参数:(http://spades.bioinf.spbau.ru/release3.5.0/manual.html#sec3.4
-1 forward reads
-2 reverse reads

--pacbio     
File with PacBio CLR reads.

--trusted-contigs 
Reliable contigs of the same genome, which are likely to have no misassemblies and small rate of other errors (e.g. mismatches and indels). This option is not intended for contigs of the related species.

--untrusted-contigs
 Contigs of the same genome, quality of which is average or unknown. Contigs of poor quality can be used but may introduce errors in the assembly. This option is also not intended for contigs of the related species.

To correct and assemble the reads:

spades.py -k 21,33,55,77,99,127 --careful  -o spades_output

实际例子
 spades.py --untrusted-contigs assembly_out_velvet/contigs.fa -k 21,23,25,27,29,31,33,47,51,55,77,81,83,85,87,97 --cov-cutoff auto --careful -1 correct.1_1.fastq.paired.fq -2 correct.1_2.fastq.paired.fq -o spades_erro

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有