1，进入R，安装Rsamtools

BiocManager::install(“Rsamtools”)

2，载入Rsamtools包

        library(Rsamtools)

3，读取bam文件

        testPairedEndBam("test.bam")

4，返回结果

若返回TRUE，则test.bam为双端；返回FALSE，则test.bam为单端

转载本文请联系原作者获取授权，同时请注明本文来自陈明杰科学网博客。

链接地址：https://wap.sciencenet.cn/blog-707141-1375276.html?mobile=1

收藏

分享