bam.py 文件源码

def get_number_reads(bam_filename):
    return reduce(lambda x, y: x + y, [ eval('+'.join(l.rstrip('\n').split('\t')[2:]) ) for l in pysam.idxstats(bam_filename) ])

# NOTE!! pgs argument used to be called pg and take a single pg item
# When resolving to tenkit:master, clients that inject PG tags must be updated to pass a list