在做这个去除PCR重复reads时候必须要明白为什么要做这个呢？WGS？WES？RNA-SEQ?CHIP-SEQ？都需要吗？随机打断测序才需要？特异性捕获不需要？

搞明白了，我们就开始做，首先拿一个小的单端测序数据比对结果来做测试！

samtools rmdup -s tmp.sorted.bam tmp.rmdup.bam

[bam_rmdupse_core] 25 / 53 = 0.4717 in library

我们的测试数据里面有53条records根据软件算出了25条reads都是PCR的duplicate，所以去除了！