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

Possible to use custom species (non-human/non-mouse)? #26

Closed
slukowski opened this issue Jun 27, 2022 · 1 comment
Closed

Possible to use custom species (non-human/non-mouse)? #26

slukowski opened this issue Jun 27, 2022 · 1 comment

Comments

@slukowski
Copy link

Hi Jamie

I'd like to be able to use mouse constant regions (TRAC/TRBC) with the human TRA/TRB genes, but am not able to do so. It wasn't clear whether I could do this by adding custom regions to the additional-genes.fasta in combination with '-xg', and editing the thimble file to list the mouse genes also didn't work.

I thought of creating a custom human/mouse hybrid reference, where I replaced human TRAC and TRBC genes with the mouse versions in the imgt, and also in the 'C-region-motifs.tsv' file. To try and use this, I called the folder 'HUMAN_mCONST' and specified this with the '-s' option. But no luck here either. I couldn't see in your code why this wouldn't work, as it should only be looking for the named folder in the Data directory, but I get the following error:

  File "/data/github/stitchr/Scripts/thimble.py", line 137, in <module>
    raise IOError("No data available for requested species: " + input_args['species'])
OSError: No data available for requested species: HUMAN_mCONST

Would you have any advice on how or if this might be possible?

Thanks!
Sam

@slukowski
Copy link
Author

slukowski commented Jun 27, 2022

Update: Changing the name to all upper case (HUMAN_mCONST --> HUMAN_MCONST) lets it pass without the error.

This works well now, and the key was to add the mouse TRAC/TRBC genes to the 'additional-genes.fasta' file and specify -'xg' in the command.

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