ChIA-PET2报错系列——第六步报错

最新推荐文章于 2022-04-29 16:32:22 发布

很酷的女超人

最新推荐文章于 2022-04-29 16:32:22 发布

阅读量365

点赞数

分类专栏：生信学习笔记 ChIA-PET2

本文链接：https://blog.youkuaiyun.com/viola_/article/details/117716690

版权

博客内容介绍了在处理ChIA-PET2数据时遇到的第六步错误，目标是获取四个特定结果。内容提到了已成功得到的前两个结果，包括组内11列数据和组间数据的详细情况。MICC文件包含11+2列，其中PP和FDR列是基于贝叶斯混合模型的评估指标。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

[bwa_aln_core] convert to sequence coordinate... 0.41 sec
[bwa_aln_core] refine gapped alignments... 0.07 sec
[bwa_aln_core] print alignments... 0.09 sec
[bwa_aln_core] 33353877 sequences have been processed.
[main] Version: 0.7.17-r1188
[main] CMD: bwa samse hg19.fa OUTdir15/index-1_2.valid.fastq.sai OUTdir15/index-1_2.valid.fastq
[main] Real time: 266.367 sec; CPU: 265.635 sec

[06-08 19:07:55] Running Step 3: Build Bedpe ...

buildBedpe OUTdir15/index-1_1.valid.sam OUTdir15/index-1_2.valid.sam OUTdir15/index-1.bedpe 30 1 1
Running buildbedpe...
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 2000000 PET
Processed 1353876 PET
All	33353876
Reads1 Low MAPQ	22004694	65.97%
Reads1 Unmapped	8067789	24.19%
Reads2 Low MAPQ	22752786	68.22%
Reads2 Unmapped