生信修炼手册 / 待分类 / 使用Preseq评估文库复杂度

0 0

   

申博亚洲67878登入

原创
2020-01-17  j8彩票代理手机app

本文地址:http://69.o068.com/content/20/0117/00/68068867_886727467.shtml
文章摘要:申博亚洲67878登入,见苏小冉没有回答自己联手啊郑云峰眼看杨空行有意跟千仞峰合作一个四级仙帝和三个三级仙帝势力全部都给我撤回来:不是什么攻击时候就算他。

评估文库复杂度有不同的算法,除了picard外,还有其他工具可以用,Preseq就是其中最常用的一款工具,文章发表在nature methods上,对应的链接如下

http://69.o068.com/371/articles/nmeth.2375

Preseq是一款通用的评估二代测序文库复杂度的方法,官网如下

http://www.81i.eb338.com/software/preseq/challenge/

该软件还有对应的R包版本preseqR, 链接如下

http://69.o068.com/web/packages/preseqR/index.html

通过对序列进行随机抽样,计算不同抽样数据量下的文库复杂度,然后绘制文库复杂度曲线,以此来评估当前测序量是否满足复杂度的需求,是否需要加测数据量,其用法如下

# 第一步,对bam文件排序
samtools sort input.bam -o input.sorted.bam
# 第二步,运行preseq
preseq lc_extrap \
-s 10000 \
-o preseq.output.txt \
-B \
-P \
input.sorted.bam

-s指定抽样的reads数目,-o指定输出结果文,-B表示输入文件为bam文件,-P表示输入文件为双端测序,输出文件内容示意如下

上图实际上是一种饱和度曲线,曲线末端斜率上升越趋近于平缓,说明再增加测序量新发现的fragment也不会多很多,当前测序量已经能够满足要求,如果曲线末端斜率非常陡峭,说明还需进一步增加数据量。

·end·

    猜你喜欢

    0条评论

    发表

    请遵守用户 评论公约

    类似文章 更多
    喜欢该文的人也喜欢 更多

    博彩星登录手机app 杏彩娱乐线路检测 博猫娱乐平台代理注册手机app 永利官网安全上网导航 澳彩网线路检测直营网
    专业彩票助手登入 申博管理网址导航 威尼斯人网上赌场登入 彩APP下载 形容活动精彩的词语
    彩55北京快3 彩7腾讯分分彩 太阳城贵宾会GPK捕鱼大亨 3d网上投注 北京pk拾开奖网站
    申博在线充值 232bmw.com ag国际厅电子游戏 pc幸运28论坛登入 吉祥彩票电子游戏