请问在不考虑测序质量的情况下，如何根据bam文件获得比对到指定position的碱基？

2: 关注

1279: 浏览

请问在不考虑测序质量的情况下，如何根据bam文件获得比对到指定position的碱基？

好问题 0 评论收藏举报

usertq 初级会员用户来自于: 上海市杨浦区
2023-08-10 17:04

chatgpt写的有些繁琐，直接
```
samtools mpileup demo.bam > demo.mpileup

```
demo.mpileup，每列分别为：chr，pos，ref（不用-f选项会显示N），深度(Number of reads covering this position)，碱基，碱基对应的质量。
具体堆积文件的格式见 samtools-mpileup(1) manual page (htslib.org)

这种堆积文件一般作为后面bcftool call或者varscan call的输入文件，或者是用来统计bam平均深度之类。

如有错误，还请指教。

赞同 2 1评论

关于作者

: m_c101dc02f446cf5872a522aeb243f264 注册会员
这家伙很懒，还没有设置简介

0: 回答

0: 文章

3: 问题

问题动态

发布时间: 2023-08-08 17:34

更新时间: 2023-08-10 17:04

关注人数: 2 人关注

相关问题

请问芯片数据分析的一般流程及涉及的常用算法？: 1043 浏览 2 关注 3 回答 0 评论

x <- x[keep.exprs,, keep.lib.sizes=FALSE] 请问一下这条命令该如何解读呢？: 537 浏览 1 关注 0 回答 2 评论

关于蛋白质文件的疑问: 798 浏览 2 关注 1 回答 0 评论

用ls *sra |while read id; do fastq-dump --split-3 $id;done命令，去解压缩sra文件，出现部分报错是什么原因？: 1025 浏览 3 关注 3 回答 0 评论

根据GFF和fasta等文件提取某一基因的ATG位点信息，比如具体的位置？: 984 浏览 2 关注 1 回答 0 评论

请问一下，R语言环境下，如何在dataframe数据中添加一个key列呢？具体描述见正文: 1063 浏览 4 关注 4 回答 1 评论

chipseq原始文件fastq如何截取DNA长度小于等于260bp进行后续mapping: 850 浏览 2 关注 2 回答 0 评论

请问RNA-seq采用poly A(+)策略建库，处理数据时若不去除rRNA会对后续分析有何影响？: 767 浏览 2 关注 1 回答 0 评论

请问hisat2+stringtie+deseq2的标准流程: 1147 浏览 2 关注 2 回答 0 评论

chip-seq 绘制bed 文件区域的 h3k27ac等修饰的信号强度图: 726 浏览 2 关注 1 回答 0 评论

请问在不考虑测序质量的情况下，如何根据bam文件获得比对到指定position的碱基？

请问在不考虑测序质量的情况下，如何根据bam文件获得比对到指定position的碱基？

0 评论

1 评论

关于作者

问题动态

相关问题

推荐内容