2
关注
1031
浏览

Picard Markduplicate

查看全部 2 个回答

chatGPT机器人 机器人 用户来自于: IANA
2023-07-31 11:26

回答:

1. Picard MarkDuplicates是根据以下原理和条件进行标记去重复的: - Picard MarkDuplicates根据测序读取的起始位置、方向、长度和序列等信息来确定是否为重复读取。 - 它使用了一种基于比对位置和序列相似性的算法来鉴定重复读取。 - MarkDuplicates会将第一个出现的读取标记为原始,而将其他相同位置的重复读取标记为重复。 - 重复读取的定义是:在同一样本中,起始位置和方向完全相同,且序列相似度高于设定的阈值。 2. 统计deduplicate reads的数值可以使用Picard工具中的CollectDuplicateMetrics来完成。具体步骤如下: - 首先运行Picard MarkDuplicates命令,得到去重后的bam文件。 - 然后使用Picard CollectDuplicateMetrics命令,对去重后的bam文件进行统计。 - CollectDuplicateMetrics会生成一个metrics文件,其中包含了各种统计信息,包括去重前后的reads数目、重复reads的数目等。 使用以下命令可以进行统计: ``` java -jar picard.jar CollectDuplicateMetrics \ I=input.bam \ O=metrics.txt \ ASSUME_SORTED=true ``` 其中,input.bam为去重后的bam文件,metrics.txt为生成的统计文件。 通过查看metrics文件中的相关统计信息,可以得到deduplicate reads的数值。

关于作者

li-nwafu 注册会员

这家伙很懒,还没有设置简介

问题动态

发布时间
2023-07-31 11:24
更新时间
2023-08-15 17:14
关注人数
2 人关注

相关问题

Picard fixmate是否有必要

推荐内容

生信行业目前薪资怎样?谈谈吧
关于RSEM和RPKM
Bowtie2使用中遇到的内存不足问题
测了实验组与对照组的mRNA与miRNA的转录组,miRNA-mRNA互作分析网络怎么处理呢?
二代测序得到测序数据组装完成后如何进行丰度计算
GWAS和WGS的选择消除分析有什么区别呀?
生存分析KM-plot交叉问题
人工智能技术目前在生信分析中有哪些应用?
ATAC-seq与RNA-seq的联合分析
All Rights Reserved Powered BY WeCenter V4.1.0 © 2024