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
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