标签:sort mergeSAM 排序 Pro HiC 报错 bam xxxR1
报错信息
Forward and reverse reads not paired. Check that BAM files have the same read names and are sorted.
分析
提示是没有排序或者R1、R2没有配对,手动sort再分步运行HiC-Pro
samtools sort -@ 10 -n -m 10G -T temp1 xxxR1.bam > xxxR1.sort.bam
samtools sort -@ 10 -n -m 10G -T temp2 xxxR2.bam > xxxR2.sort.bam
mv xxxR1.sort.bam xxxR1.bam
mv xxxR1.sort.bam xxxR1.bam
注意
- 不能按染色体位置排序,因为HiC-Pro在处理HiC数据这一步需要配对关系,所以需要让R1、R2的bam按名称排序
- -m不能设置过低,HiC测序深度高的话会导致排序的过程文件过多而不能合并
- 最好设置-T,因为此时IO交互量很大
- 最后把排好序的文件按原来的名字修改命名
参考
标签:sort,mergeSAM,排序,Pro,HiC,报错,bam,xxxR1 来源: https://www.cnblogs.com/johnsonzzz/p/15716083.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。