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

Error in "cellranger mkref" #244

Closed
grimwoo opened this issue Mar 15, 2024 · 1 comment
Closed

Error in "cellranger mkref" #244

grimwoo opened this issue Mar 15, 2024 · 1 comment

Comments

@grimwoo
Copy link

grimwoo commented Mar 15, 2024

Thanks for the development of cellranger-7.1.0. I successfully created the count files of human, mouse, and etc. via "cellranger count" command. However, I encountered an error in making genome reference for axolotl (salamander).

The gtf file was downloaded from https://www.axolotl-omics.org/dl/AmexT_v47-AmexG_v6.0-DD.gtf.gz.
The fa file was downloaded from https://www.axolotl-omics.org/dl/AmexG_v6.0-DD.fa.gz.

I tried to fix the gtf file according to your recommendation:
https://kb.10xgenomics.com/hc/en-us/articles/4707448154381-Common-mkref-formatting-errors-when-building-custom-reference-from-NCBI-UCSC-or-RefSeq-genomes
https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/advanced/references#ref-pack

Using the filtered gtf file ("AmexT_v47-AmexG_v6.0-DD_fittered.gtf"), I still encountered error after the following second command.

cellranger mkgtf \
  AmexT_v47-AmexG_v6.0-DD_fittered.gtf \
  AmexT_v47-AmexG_v6.0-DD.filtered.again.gtf \
  --attribute=gene_biotype:protein_coding
cellranger mkref \
  --genome=Axolotl_genome \
  --fasta=AmexG_v6.0-DD.fa \
  --genes=AmexT_v47-AmexG_v6.0-DD.filtered.again.gtf
Indexing genome FASTA file...
...done

Writing genes GTF file into reference folder...
WARNING: The following transcripts appear on multiple chromosomes in the GTF:
TMEM11
SGMS1
HTR3E
SELENOO
SYNJ2BP
NXPE2
GNGT1
ARPC1A
MYH2
N/A
PRDM9

This can indicate a problem with the reference or annotations. Only the first chromosome will be counted.
...done

mkref has failed: error building reference package
Error detected in GTF file: Previous 'exon' entry for transcript_id ACER1 had strand == (-), but 'transcript' entry has strand == (+)

Could you please help me to produce a reference for axolotl (salamander).

Copy link

Thank you for your interest in Cell Ranger. GitHub Issues are not monitored by our support team. Please e-mail [email protected] for support with Cell Ranger or visit our documentation site at https://www.10xgenomics.com/support/software/cell-ranger.

@github-actions github-actions bot closed this as not planned Won't fix, can't repro, duplicate, stale Mar 15, 2024
@github-actions github-actions bot locked and limited conversation to collaborators Mar 15, 2024
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant