检测bam文件的完整度 samtools view T_recal.bam|head samtools view T_recal.bam|tail for i in *.bam ;do (samtools quickcheck $i && echo "ok" || echo $i error);done