利用picard 进行序列去重复

for id in $(ls /nfs_genome/home/zhaoshanzhong/zs/data/metagenome64/metacleandata/mapOnly/work_files/*.filtered.bam | sed "s/.filtered.bam//"); do base=`basename $id`;echo "picard -Xmx256g MarkDuplicates -I $id.filtered.bam -O $base.markdup.bam -M  $base.markdup -REMOVE_DUPLICATES true"; done > markdu.sh

nohup ParaFly -c markdu.sh -CPU 20 &