samtools flags 分解

　　bam文件的第二列的flag值，代表了当前read比对到基因组序列的基本情况。flag值由2的0次方到2的11次方这11个数字组成，即1，2，4，8，16，32，64，128，256，512，1024，2048。每个数字表示的比对情况这里不赘述，当获得当前read的flag值时，如何拆解成flag最初的值呢，这里提供一个脚本，权当记录备忘了，大家有更便捷的方法，欢迎交流。

程序主体如下：

#!/usr/bin/perl
use strict;
use autodie;
use warnings;

sub o2bn{
        my $n = shift;
        my $b = sprintf("%b", $n);
        my @str = split //, reverse $b;
        my @result;
        for(my $i=0; $i<@str; $i++){
                if($str[$i] eq '1'){
                        push @result, 2**$i;
                }
        }
        return @result;
}

print "please input flag number [quit:0]: ";
while(my $number = <STDIN>){
        last if $number == 0;
        my @ret = o2bn($number);
        print "result is : @ret\n";
        print "please input flag number [quit:0]: ";
}

结果：