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

Adding a gene for a custom reference #273

Open
eigenP opened this issue Oct 23, 2024 · 1 comment
Open

Adding a gene for a custom reference #273

eigenP opened this issue Oct 23, 2024 · 1 comment

Comments

@eigenP
Copy link

eigenP commented Oct 23, 2024

Describe the issue
I'm making a custom reference (for the chicken genome) and I also want to add the GFP sequence. I have successfully done so following the instructions for cellranger (https://www.10xgenomics.com/support/software/cell-ranger/latest/tutorials/cr-tutorial-mr) (and the analysis works as expected with cellranger), but after inspecting the final outputs of kb ref and count, there is no GFP (also no trace of it in t2g.txt, etc.)

Any pointers (e.g. do I need to add more fields in the GFP gtf or FASTA than what the 10x tutorial) would be greatly appreciated! Thank you in advance.

What is the exact command that was run?

kb ref -i index.idx -g t2g.txt -f1 cdna.fasta --overwrite Gallus_gallus.bGalGal1.mat.broiler.GRCg7b.dna_GFP.toplevel.fa Gallus_gallus.bGalGal1.mat.broiler.GRCg7b.112.filtered.GFP.gtf
@Yenaled
Copy link
Collaborator

Yenaled commented Oct 24, 2024

You need to put more fields into the GTF. Take a look at an example for another gene and use that as a template for making your custom GTF entry.

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

2 participants