samtools view筛选cellranger比对结果 - 做专业的生信问答网站

登录

2: 关注

1952: 浏览

samtools view筛选cellranger比对结果

请问下，如果我想把cellranger产生的bam文件进行一个过滤，就是只保留包含barcode(tag值是CB:Z)和UMI(tag值是UB:Z)的比对结果，用samtools应该怎么操作呢？

好问题 0 评论收藏举报

查看全部 1 个回答

chatGPT机器人机器人用户来自于: IANA
2024-03-07 15:11

要筛选出只包含barcode和UMI信息的比对结果，可以使用samtools view命令结合-G参数进行操作。

具体操作如下：

```

samtools view -h input.bam | awk '{if($1 ~ /^@/ || ($1 ~ /^[^@]/ && (($14 ~ /CB:Z:/ && $16 ~ /UB:Z:/) || ($14 ~ /UB:Z:/ && $16 ~ /CB:Z:/)))) {print $0}}' | samtools view -b -o output.bam -

```

这条命令的意思是：

1. 使用samtools view读取输入的bam文件，并通过管道传递给awk进行过滤。

2. awk命令中的条件判断保留了包含barcode和UMI信息的比对结果。

3. 最后通过samtools view将过滤后的结果输出到output.bam文件。

执行这条命令后，output.bam文件中将只包含有barcode和UMI信息的比对结果。

赞同 0 0评论

关于作者

: bioinfo_2023 注册会员
生信新手

0: 回答

0: 文章

11: 问题

问题动态

发布时间: 2024-03-07 15:04

更新时间: 2024-03-07 15:11

关注人数: 2 人关注

相关问题

GATK和Samtools联合寻找变异: 2560 浏览 2 关注 2 回答 0 评论

lnc筛选: 1547 浏览 2 关注 1 回答 0 评论

分析CRISPR 高通量筛选数据: 2305 浏览 1 关注 2 回答 0 评论

Bulk-RNAseq多组数据差异表达基因的筛选: 2478 浏览 2 关注 2 回答 0 评论

当LIMMA计算计算结果中的调整P值都大于0.05时，可否使用p值<0.05+fc>2来筛选差异表达基因？: 2344 浏览 2 关注 1 回答 0 评论

问题：如何筛选WGCNA分析得出的consensus module中的hub基因？: 2197 浏览 2 关注 1 回答 1 评论

两个samtools命令之间用管道连接，最后的输出文件被吞了。: 1485 浏览 1 关注 1 回答 0 评论

VCONTACT2的结果文件genome_by_genome_overview如何进行统计分类: 2061 浏览 2 关注 1 回答 0 评论

用GWAS筛选受选择基因，样本数量不够，请问可以用选择消除分析吗: 1764 浏览 2 关注 1 回答 0 评论

samtools 有错误: 1802 浏览 2 关注 1 回答 0 评论

推荐内容

问单细胞RNA-seq，如何计算不同cluster之间的相关性: 2074 浏览 2 关注 1 回答 0 评论

问 bowtie2比对结果报告best，以及aligned concordantly: 2102 浏览 2 关注 2 回答 0 评论

问 seurat包，如何调节大热图左侧的基因字体大小: 3009 浏览 2 关注 1 回答 0 评论

问 pseudobulk分析: 3231 浏览 2 关注 3 回答 0 评论

问关于scrublet的使用: 2920 浏览 2 关注 1 回答 0 评论

问你好chatGPT,可以提供包含双细胞去除的单细胞分析完整流程吗？: 2394 浏览 2 关注 2 回答 0 评论

问单细胞RNA-seq分析流程: 2166 浏览 2 关注 2 回答 0 评论

问空转bin_size: 1913 浏览 2 关注 2 回答 0 评论

问 scRNA-seq有什么推荐的监督聚类方式吗？: 1989 浏览 2 关注 1 回答 0 评论

问去除批次效应选定靶标细胞群后如何进行亚群分析？: 2788 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025