Bug in call-ataqc module

[encode-dcc]

Bug in call-ataqc module

[leepc12]

It looks like your BAM file is empty. Please provide error logs (.tar.gz).

[nicolerg]

I'm not sure if either of these are the BAM file you're referring to:
call-ataqc/shard-0/inputs/789850838/${SAMPLE}.trim.merged.bam is 1.5G
call-ataqc/shard-0/inputs/-623782124/${SAMPLE}.trim.merged.nodup.bam is 667M.

debug_64.tar.gz

[leepc12]

This will be fixed in v1.1.4. But if you want to fix it yourself then modify function get_chr_m() in line 203-221 of /users/nicolerg/anaconda2/envs/encode-atac-seq-pipeline/bin/encode_ataqc.py like the following.

def get_chr_m(sorted_bam_file):
    '''
    Get fraction of reads that are mitochondrial (chr M).
    '''
    logging.info('Getting mitochondrial chromosome fraction...')
    chrom_list = pysam.idxstats(sorted_bam_file, split_lines=True)
    tot_reads = 0
    chr_m_reads = 0
    for chrom in chrom_list:
        chrom_stats = chrom.split('\t')
        if chrom_stats[0] == 'chrM':
            chr_m_reads = int(chrom_stats[2])
        tot_reads += int(chrom_stats[2])
    if tot_reads==0:
        fract_chr_m = 0
    else:
        fract_chr_m = float(chr_m_reads) / tot_reads

    return chr_m_reads, fract_chr_m

[nicolerg]

That function is not in encode_ataqc.py, but I modified it in /users/nicolerg/anaconda2/envs/encode-atac-seq-pipeline/bin/run_ataqc.py.