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

Where has UMI count for AIR chains gone? #478

Closed
adamwoolfe opened this issue Jan 11, 2024 · 1 comment · Fixed by #487
Closed

Where has UMI count for AIR chains gone? #478

adamwoolfe opened this issue Jan 11, 2024 · 1 comment · Fixed by #487
Labels
bug Something isn't working

Comments

@adamwoolfe
Copy link

Describe the bug
I am looking to access UMI count for each chain of a set of antibodies read into the scirpy anndata frame. It appears that read count is present in the awkward array associated with the airr anndata obsm awkward array under the title 'consensus_count' like so:

adata["airr"].obsm['airr']['consensus_count']

but the UMI count is nowhere to be found. None of the available fields point to the UMI count.

where has this field gone (it is present in the original 10X filtered_contig_annotation.csv input file)?

System

  • OS: Linux
  • Python version: 3.10.3
  • Versions of libraries:
    anndata 0.10.4
    awkward 2.5.1
    matplotlib 3.8.2
    mudata 0.2.3
    muon 0.1.5
    numpy 1.26.2
    pandas 2.1.3
    scanpy 1.9.3
    scirpy 0.14.0

Additional context
Add any other context about the problem here.

@adamwoolfe adamwoolfe added the bug Something isn't working label Jan 11, 2024
@grst
Copy link
Collaborator

grst commented Jan 11, 2024

Hi Adam,

thanks for bringing this up!
The UMI count goes to adata["airr"].obsm["airr"]["duplicate_count"].

I always found this a bit confusing, and I've just seen that there has been an update to the AIRR Rearrangement standard that introduces a dedicated umi_count field which is a lot clearer. I'll update scirpy to use this field instead in the future.

Best,
Gregor

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
Status: Done
Development

Successfully merging a pull request may close this issue.

2 participants