czbiohub-sf / tabula-muris-senis

Tabula Muris Senis
http://tabula-muris-senis.ds.czbiohub.org
BSD 3-Clause "New" or "Revised" License
93 stars 26 forks source link

h5ad to BAM matching #52

Open davek44 opened 7 months ago

davek44 commented 7 months ago

Hi, perhaps related to this issue in human https://github.com/czbiohub-sf/tabula-sapiens/issues/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!