pp.import_data()

[pass-c]

Hello, I am newer to analyzing snATACseq data. I am in the initial preprocessing stages and have a very basic question. Some background I am importing fragment files that were generated in from a 10x Multiomic experiment and cellranger-arc.

When I run pp.import_data() the output: AnnData object with n_obs x n_vars = 50540 x 0. In scRNA-seq the n_obs is the nuclei count, does the same go for ATAC? If it is the case and n_obs is cell/nuclei count, could there be an error in the fragment file input (as I know the cell count and the output here does not align with that).

Thank you for your assistance.

[kaizhang]

import_data imports all barcodes/cells. Many of them may not be valid cells. You need to filter them based on number of reads or TSS enrichment.

[pass-c]

Thank you. I have come across your resources for generating snap files. Is this something that needs to be generated from my 10x Genomics cellranger output prior to running snapATAC2 in python? Additionally how do you get the .snap file to a .hdf5 file for downstream analysis? Thank you for your assistance.