2
关注
1473
浏览

samtools view筛选cellranger比对结果

查看全部 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信息的比对结果。

问题动态

发布时间
2024-03-07 15:04
更新时间
2024-03-07 15:11
关注人数
2 人关注

相关问题

samtools 有错误
lnc筛选
当LIMMA计算计算结果中的调整P值都大于0.05时,可否使用p值<0.05+fc>2来筛选差异表达基因?
GATK和Samtools联合寻找变异
cellranger运行结果分析
分析CRISPR 高通量筛选数据
WGCNA筛选到感兴趣的模块后,如何找这个模块的hub基因呢?
cellranger使用问题
两个samtools命令之间用管道连接,最后的输出文件被吞了。
Bulk-RNAseq多组数据差异表达基因的筛选

推荐内容

使用seurat包,导出特定cluster的细胞-基因counts矩阵。初始数据命名为pbmc
​单细胞RNA-seq,如何计算不同cluster之间的相关性
空转bin_size
单细胞——小鼠T细胞注释的markers
你好chatGPT,可以提供包含双细胞去除的单细胞分析完整流程吗?
reads数越多,滑窗内GC含量降低
单细胞RNA-seq分析流程
singularity查到不到指定输入文件位置
monocle3中plot_cells更改其画出来的细胞群的颜色
celseq2转换单细胞原始数据
All Rights Reserved Powered BY WeCenter V4.1.0 © 2025