bam文件去除pcr重复和去除线粒体信息

最新推荐文章于 2025-04-09 14:33:24 发布

小兔叽乖乖

最新推荐文章于 2025-04-09 14:33:24 发布

阅读量1k

点赞数

文章标签： linux

本文链接：https://blog.youkuaiyun.com/m0_70611557/article/details/130380757

版权

sambamba markdup -r input.bam output.bam

samtools view -h input.bam |grep -v "chrM" |samtools view -bS -o output.bam

确定要放弃本次机会？

福利倒计时

: :

立减 ¥

普通VIP年卡可用

立即使用

小兔叽乖乖

关注关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
1
评论
分享

复制链接

分享到 QQ

分享到新浪微博

扫一扫
举报

举报

生信技能26 - BWA比对、去除PCR重复和统计最佳实践

LittleComputerRobot的博客

08-14

568

BWA MEM 比对、去除PCR重复和统计最佳实践

生信技能56 - 去除重复BAM文件的窗口reads计数方法

LittleComputerRobot的博客

08-14

332

去除重复BAM文件的窗口reads计数方法

1 条评论您还未登录，请先登录后发表或查看评论

生信软件23 - Samtools和GATK去除PCR重复方法汇总

LittleComputerRobot的博客

06-21

1968

Samtools和GATK去除PCR重复方法汇总

归纳:bam去重

qq_39306047的博客

03-10

2984

####rmdup一般已经过时，一般不要用这个方法！ picard Picard也是一款去重的软件，只不过要用Java来启动它，所以你的服务器上要安装Java 下载地址：https://github.com/broadinstitute/picard/releases/tag/2.21.6 选择第一个，直接下载就可以用了但是利用picard去重所用的sam文件或者bam文件是需要加header的，如果你的文件没有，我们可以用samtools加个表头： #加header samtools view -Sb

frp 进行内网穿透Sambamba 去除重复

qq_39306047的博客

04-14

411

通过 frp 进行内网穿透并实现Rstudio server远程访问https://www.jianshu.com/p/52956240685b Sambamba 去除重复工具https://www.jianshu.com/p/18fc84d4fbba

零基础小白笔记7 | 对bam文件进行去重和过滤

tmrtmr___的博客

01-01

6068

1.samtools index是 Samtools 工具包中的一个命令，用于为 BAM 文件创建索引，以便于快速的随机访问。通过 Picard 工具，可以识别和标记 PCR 重复，这对于测序数据的后续分析是很重要的。（2.1）markduplicates：是 Picard 工具集中的一个命令行工具，主要用于识别和标记测序数据中的 PCR 重复。1.samtools view 是 Samtools 工具包中的一个命令，用于查看和转换 SAM/BAM 格式的文件，也可用于过滤掉质量比较低的序列；

ATAC-seq数据分析流程

weixin_57975238的博客

12-31

1326

ATAC-seq分析流程

二代和三代测序数据的生殖系变异分析最佳实践综述

最新发布

LittleComputerRobot的博客

04-09

1038

二代和三代测序数据的生殖系变异分析最佳实践综述

16/10/2019 一步步学会分析ATAC-seq

uqyuan0515的博客

10-29

2959

https://zhuanlan.zhihu.com/p/20702684 关于illumina 测序的原理 so, 我的测序结果是paired-end sequencing 并且是在不同的lanes上测序的。原始名称为： larvae1_S1_L001_R1_001.fastq.gz, larvae1_S1_L001_R2_001.fastq.gz, larvae1_S1_L002...

强烈推荐！Encode官方的ATAC数据分析流程

庐州月光的博客

01-19

2631

欢迎关注”生信修炼手册”!Encode不仅共享了大量的组学数据，还开源了自己的数据分析pipeline, ATAC的pipeline网址如下https://github.com/ENCO...

samtools and bwa

03-11

samtools and bwa参数例子

bamtools拆分bam文件

weixin_46530443的博客

12-07

1521

bamtools拆分bam文件 sRNA使用shortstack比对到基因组上的比对文件是合并的bam文件，需要按照样品将其拆分，使用bamtools进行拆分： bamtools split -in merged_alignments.bam -tag RG 运行时间较久，有更好的方式再来更新。 ...

使用sambamba去重遇到的问题

探序基因的博客

10-24

778

2>>log/picard_MarkDuplicates_log.txt & # 将标准错误输出追加到日志文件，后台执行。M=metrics.txt \ # 重复读数的统计信息文件。注意java版本，如果版本不满足，会报错，直接下载满足条件的java环境就行。MAX_FILE_HANDLES=800 \ # 最大文件句柄数。I=sample_sorted.bam \ # 输入BAM文件。O=sample_rmDup.bam \ # 输出BAM文件。探序基因肿瘤研究院整理。

生信c语言面试题,生信（十一）htslib处理bam或sam文件的简单示例

weixin_35174422的博客

05-19

703

原创：hxj7本文给出了一个示例，介绍如何使用htslib编写c程序来处理bam/sam文件。(本文写于2020年初，随着将来htslib和samtools库的更新，本文部分内容可能会不适用，请读者注意官网的更新动态。)我们通常会使用samtools软件来处理bam/sam文件。但有时候我们也需要对bam/sam文件进行一些个性化的处理，这个时候就需要自己编写程序。从samtools的github...

如何用c语言读bam文件,【一起学生信】bam文件统计覆盖深度、靶向捕获效率

weixin_30618305的博客

05-23

668

ItemAnnotation[Total] Raw Reads (All reads)All reads in the bam file(s).[Total] QC Fail readsReads number failed QC, this flag is marked by other software,like bwa. See flag in the bam structure.[Tota...

samtools sort -@ 4 -t CB -O BAM -o $new $idsamtools排序原理过程

生信小博士的博客

03-26

5094

读取每个小文件的第一个read，将它们按照位置大小排序，选出最小的一个read，将该read写入输出文件中，并从对应的小文件中读取下一个read，并继续进行排序，重复此过程直到所有read都被排序。因此，在进行单细胞的scvelo分析时，通常需要对BAM文件进行排序，以确保后续的数据分析结果的准确性和可靠性。读取每个小文件的第一个read，将它们按照位置大小排序，选出最小的一个read，将该read写入输出文件中，并从对应的小文件中读取下一个read，并继续进行排序，重复此过程直到所有read都被排序。

sambamba — samtools 的高效平替工具

long_1998的博客

02-23

2037

sambamba 不仅可以对 BAM 文件进行排序、索引、过滤、统计、标记重复等常见的操作，还可以进行一些特殊的功能，比如区域过滤，标记重复序列，检测结构变异等。flagstat: 统计 BAM 文件中的 reads 的标志位，可以显示 QC 通过和失败的 reads 数目，以及各种配对情况和重复情况。view: 查看 SAM/BAM 文件的内容，可以指定输出格式和过滤条件，也可以查看参考序列的信息。提供一个或多个覆盖度阈值，对于每个阈值，会添加一个额外的列，显示区域中覆盖度超过此值的碱基的百分比。

samtools用法详解

热门推荐

sunchengquan的博客

12-21

1万+

文章目录下载安装测试数据命令详解dictfaidxindexreheaderrmdupcatmergempileup查看参数mpileup生成的结果有参考序列的pileup使用生成一个简单的vcf文件sort查看用法主要参数释义：splitfastqfastabedcovdepthflagstatidxstatsstatsflagstviewviewbam文件转换为sam文件sam文件转换为bam...

NGS测序中PCR重复序列的判定方法

庐州月光的博客

01-04

6229

欢迎关注”生信修炼手册”!在NGS的数据分析中，去除PCR重复序列是一个常见的分析步骤，无论是WES/WGS的snp calling，还是chip_seq, ATAC_seq，都需要对原...

如何去除bam 文件中的pcr 重复

03-31

PCRs (聚合酶链式反应) 重复可以通过多种方法去除，以下是一些常用的方法： 1. Picard 工具 Picard 是一个 Java 工具包，可用于处理 BAM 文件。其中，MarkDuplicates 工具可以检测和去除 PCR 重复。该工具可以标记 PCR 重复的序列，也可以将它们从 BAM 文件中删除。安装 Picard 工具： ``` conda install -c bioconda picard ``` 使用 MarkDuplicates 工具： ``` picard MarkDuplicates \ I=input.bam \ O=output.bam \ M=metrics.txt \ REMOVE_DUPLICATES=true ``` 2. SAMtools 工具 SAMtools 是一个流行的工具包，用于处理 SAM 和 BAM 文件。其中，rmdup 命令可以检测和去除 PCR 重复。安装 SAMtools 工具： ``` conda install -c bioconda samtools ``` 使用 rmdup 命令： ``` samtools rmdup input.bam output.bam ``` 3. GATK 工具 GATK (基因组分析工具包) 是一个广泛使用的工具包，用于变异检测和基因组分析。其中，MarkDuplicates 工具可以检测和去除 PCR 重复。安装 GATK 工具： ``` conda install -c bioconda gatk4 ``` 使用 MarkDuplicates 工具： ``` gatk MarkDuplicatesSpark \ -I input.bam \ -O output.bam \ --remove-sequencing-duplicates \ --metrics-file metrics.txt ``` 以上是一些常用的方法，可以根据实际需求选择适合的工具和方法。