摘要: Samtools merge 时出错:弄了个软链接,如图2,使用链接过去的文件时出现图1的错误,主要是软链接的时候没有全部提供绝对路径。出现以下错误,是因为事先没提前对bam文件进行sort,genomeCoverageBed中的输入文件同样需要经过排序,如若是bed,则须用sort排序,如若是ba... 阅读全文
posted @ 2015-07-16 16:59 凉皮子 阅读(534) 评论(0) 推荐(0) 编辑
摘要: 有很多种格式来记录mapping的结果,大多数都收录在UCSC的帮助文档中。在mapping之后,如何对mapping的结果进行评估呢? 最简单的,就是通过samtools来评估mapping质量了。samtools idxstats aln.sorted.bam注意,这一步之前需要经过sort和i... 阅读全文
posted @ 2015-07-16 11:19 凉皮子 阅读(5418) 评论(0) 推荐(0) 编辑