最近在跑程序的时候遇到一个问题，简直不要太坑，记录下来：
[W::bam_hdr_read] EOF marker is absent. The input is probably truncated.
em.......看到这个的时候，也就意味着重跑了。一般bam文件跑完出错的概率比较小，如果出错了的话，就得重新返工了。所以检查一下bam文件的完整性比较好！
如何检查呢？
方法一：
利用samtools quickcheck
samtools quickcheck 帮助文档
代码实战
方法二：来自于健明大神的帖子：https://www.jianshu.com/p/c6dd7edd6e80 这个方法一行代码就可以解决~（虽然上一个也是一行代码呀）
tail -c 28 SRR1286778_1.Cri.bam| xxd -p
输出结果
这个才是正常的bam的输出结果，如果出现了别的数字（乱码啥的。。。。）重跑吧。
最近在跑程序的时候遇到一个问题，简直不要太坑，记录下来： em.......看到这个的时候，也就意味着重跑...
一.数据完整性 数据检查：一般使用CRC -32方式检查数据完整性 1.对本地文件I/o的数据完整性检查 hdfs...
参考文章：1.如何统计BAM文件中的reads数2.Samtools常用命令的总结 当你有很多个bam文件时，想知...
bam文件的排序方式针对不同的应用场景，如果不能了解它将会极大影响计算的效率，本文简单介绍如何快速查看bam文件的...
IGV需要的文件：比对文件为 bam和bam.bai文件 参考基因组文件...
做生信分析的小伙伴们，相信大家对bam文件都不陌生，但具体到如何get到bam文件提供给我们的信息，却少有人真...
前言 如果对较大的bam文件进行操作时，往往相当耗时，所以将一个很大的bam文件分割为几个较小的bam文件，在并行...
常用samtools的命令： samtools view命令：查看sam，bam文件，进行sam及bam文件的转换...
1.bam文件读取 samtools view xxx.bam samtools view xxx.bam |le...
what is bam bam（或者是sam，cram）文件，是比对后拿到的文件，sam文件是纯文本，非常占用存储...
本文标题：如何检查bam文件的完整性？
本文链接：https://www.haomeiwen.com/subject/yblefctx.html
网友评论