源码阅读网myBatis源码分析师(精剪版)
java -jar $picard BuildBamIndex \
INPUT=dedup_sample1.bam
#查看@RG
samtools view -H dedup_sample1.bam | grep '@RG'
如果同一样本(同一文库)多次上样测序,需要将多次测序结果合并后,再进行sort,mark。
#merge bam files(利用picard可以保留@RG; samtools合并时会丢掉@RG)
time java -jar $picard MergeSamFiles\
I=1.bam I=2.bam\
O=12.bam