Inconsisent shape of region_annotation vs. region_annotation_b

[fwzhao]

toolkit/ngs_toolkit/atacseq.py

Line 72 in 1dfa4d0

>>> a.get_peak_genomic_location()

Running analysis.get_peak_genomic_location() leads to inconsistent number of regions between background and real peak set.

Dataset: projects/pad-map (feel free to use)

See:

In [15]: [f'{i}: shape {getattr(analysis, i).shape}' for i in analysis.dict if isinstance(getattr(analysis, i), pd.core.fram
...: e.DataFrame)]
Out[15]:
['comparison_table: shape (425, 5)',
'gene_annotation: shape (211108, 6)',
'closest_tss_distances: shape (211468, 6)',
'region_annotation: shape (211105, 4)', <-----
'region_annotation_b: shape (211104, 4)', <----
'region_annotation_mapping: shape (395590, 4)',
'region_annotation_b_mapping: shape (339561, 4)',
'chrom_state_annotation: shape (211108, 4)',
'chrom_state_annotation_b: shape (211108, 4)',
'chrom_state_annotation_mapping: shape (250895, 4)',
'chrom_state_annotation_b_mapping: shape (223218, 4)']