Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

h5ad to BAM matching #52

Open
davek44 opened this issue Feb 12, 2024 · 0 comments
Open

h5ad to BAM matching #52

davek44 opened this issue Feb 12, 2024 · 0 comments

Comments

@davek44
Copy link

davek44 commented Feb 12, 2024

Hi, perhaps related to this issue in human czbiohub-sf/tabula-sapiens#42, I have also found apparently mismatched BAMs in the mouse data.

For example, cell row 96552 in the tabula-muris-senis-facs-official-raw-obj.h5ad has 'cell' column 'L22.MAA000593.3_8_M.1.1'. The table reports 2099 n_genes and 189281 n_counts.

That appears to match s3://czb-tabula-muris-senis/Plate_seq/3_month/170907_A00111_0051_BH2HWLDMXX/results_gencode_ercc/L22-MAA000593-3_8_M-1-1.gencode.vM19.ERCC.Aligned.out.sorted.bam. However, this BAM only has 486 aligned read rows, so it can't reach those reported gene expression values.

Could you help me understand where I'm going wrong here? Thanks!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant