源碼閱讀網(wǎng)myBatis源碼分析師(精剪版)
java -jar $picard BuildBamIndex \
? ?INPUT=dedup_sample1.bam
#查看@RG
samtools view -H dedup_sample1.bam | grep '@RG'
如果同一樣本（同一文庫）多次上樣測序，需要將多次測序結(jié)果合并后，再進行sort,mark。
#merge bam files(利用picard可以保留@RG; samtools合并時會丟掉@RG）
time java -jar $picard MergeSamFiles\
I=1.bam I=2.bam\
O=12.bam
標簽：
