Linux: 如何使用samtools看你的 PE librar

test.png

ATAC-seq 的insert size 怎么搞，来干的，直接套底下的。

samtools view P12.oLT.R1.trim.PE2SE.nodup.bam | awk '$9>0' | cut -f 9 | sort | uniq -c | sort -b -k2,2n | sed -e 's/^[ \t]*//' > ./size_density/p12.oLT.sd.txt