Namespace(genes='/oak/stanford/groups/akundaje/soumyak/refs/gencode/hg38/hg38.gencode.protein_coding.tss.bed', list='/oak/stanford/groups/akundaje/projects/neuro-variants/variant_summary/rosmap/encode_2024/K562_bias/encode_2024.LV.cardiomyocyte.H.mean.variant_scores.tsv', out_prefix='/oak/stanford/groups/akundaje/projects/neuro-variants/variant_annotations/rosmap/encode_2024/K562_bias/encode_2024.LV.cardiomyocyte.H', peaks='/oak/stanford/groups/akundaje/projects/neuro-variants/data/processed/encode_2024/peaks/overlap/encode_2024.LV.cardiomyocyte.H.overlap.peaks.bed.gz', schema='chrombpnet') Namespace(genes='/oak/stanford/groups/akundaje/soumyak/refs/gencode/hg38/hg38.gencode.protein_coding.tss.bed', list='/oak/stanford/groups/akundaje/projects/neuro-variants/variant_summary/rosmap/encode_2024/K562_bias/encode_2024.LV.cardiomyocyte.H.mean.variant_scores.tsv', out_prefix='/oak/stanford/groups/akundaje/projects/neuro-variants/variant_annotations/rosmap/encode_2024/K562_bias/encode_2024.LV.cardiomyocyte.H', peaks='/oak/stanford/groups/akundaje/projects/neuro-variants/data/processed/encode_2024/peaks/overlap/encode_2024.LV.cardiomyocyte.H.overlap.peaks.bed.gz', schema='chrombpnet') chr pos end allele1 allele2 variant_id 837619 chr1 850506 850507 A T chr1:850507:A:T 837620 chr1 850562 850563 G C chr1:850563:G:C 837621 chr1 850855 850856 C T chr1:850856:C:T 837622 chr1 851095 851096 G A chr1:851096:G:A 837623 chr1 852046 852047 C T chr1:852047:C:T Variants table shape: (2531122, 6) annotating with closest genes 0 1 2 3 4 5 6 7 8 9 \ 0 chr1 850506 850507 A T chr1:850507:A:T chr1 923923 923924 SAMD11 1 chr1 850506 850507 A T chr1:850507:A:T chr1 959308 959309 NOC2L 2 chr1 850506 850507 A T chr1:850507:A:T chr1 960584 960585 KLHL17 3 chr1 850562 850563 G C chr1:850563:G:C chr1 923923 923924 SAMD11 4 chr1 850562 850563 G C chr1:850563:G:C chr1 959308 959309 NOC2L 10 11 12 13 14 0 0 + ENSG00000187634.13 protein_coding 73417 1 0 - ENSG00000188976.11 protein_coding 108802 2 0 + ENSG00000187961.15 protein_coding 110078 3 0 + ENSG00000187634.13 protein_coding 73361 4 0 - ENSG00000188976.11 protein_coding 108746 Closest genes table shape: (7593366, 15) annotating with peak overlap chr pos end allele1 allele2 variant_id 0 chr1 850506 850507 A T chr1:850507:A:T 1 chr1 850562 850563 G C chr1:850563:G:C 2 chr1 858636 858637 C T chr1:858637:C:T 3 chr1 858749 858750 G A chr1:858750:G:A 4 chr1 866125 866126 G A chr1:866126:G:A Peak overlap table shape: (156232, 6) chr pos allele1 allele2 variant_id logfc.mean \ 0 chr10 124524 T A chr10:124524:T:A 0.004236 1 chr10 124526 T A chr10:124526:T:A 0.008436 2 chr10 124528 T A chr10:124528:T:A 0.004839 3 chr10 124530 T A chr10:124530:T:A 0.003009 4 chr10 124542 T A chr10:124542:T:A 0.013604 logfc.mean.pval abs_logfc.mean abs_logfc.mean.pval jsd.mean ... \ 0 0.731995 0.012551 0.738737 0.020696 ... 1 0.746218 0.012213 0.755924 0.020241 ... 2 0.791019 0.009688 0.799129 0.020175 ... 3 0.813530 0.008368 0.819998 0.019508 ... 4 0.670151 0.013604 0.709242 0.016616 ... quantile_change.mean.pval abs_quantile_change.mean \ 0 0.652774 0.001085 1 0.663417 0.001112 2 0.676740 0.000940 3 0.729696 0.000728 4 0.565129 0.001560 abs_quantile_change.mean.pval closest_gene_1 gene_distance_1 \ 0 0.660752 ZMYND11 5565 1 0.673643 ZMYND11 5563 2 0.685470 ZMYND11 5561 3 0.736699 ZMYND11 5559 4 0.603081 ZMYND11 5547 closest_gene_2 gene_distance_2 closest_gene_3 gene_distance_3 \ 0 TUBB8 50361 DIP2C 565144 1 TUBB8 50363 DIP2C 565142 2 TUBB8 50365 DIP2C 565140 3 TUBB8 50367 DIP2C 565138 4 TUBB8 50379 DIP2C 565126 peak_overlap 0 False 1 False 2 False 3 False 4 False [5 rows x 38 columns] Annotation table shape: (2531122, 38) DONE Done Annotating