2019-06-15 提取align后不成对的reads

samtools view LTR3-2_L4_P701501.bam | perl -alne '{print if $F[6] ne "="}' > unaligned.sam

包括未比对上的和比对到不同染色体的

cut -f 3,7 unpaired.sam |sort |uniq -c