用Linux命令行或软件对Bam文件进行操作
统计CIGAR，从多到少排序
1 | sambamba view -F "ref_name=='******' and not secondary_alignment and not supplementary" ******.bam | cut -f 6 | sort | uniq -c | sed -r 's/^\s+//' | sed -r 's/\s+/\t/' | sort -k 1 -nr > CIGAR.tsv |
提取Reads，转fasta
1 | sambamba view -F "ref_name=='******' and not secondary_alignment and not supplementary" ******.bam | cut -f 1,10 | sed 's/^/>/' | sed -r 's/\t/\n/' > Target.fa |