Hi Biostars,
I try to do velocity analysis and it needs a loom file so I need to run velocyto. To run velocyto, I need a gtf file and after reading the material below, I still get lost on how to get it. Would you please have a suggestion? Thank you so much!
https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/advanced/references
https://velocyto.org/velocyto.py/tutorial/cli.html#run10x-run-on-10x-chromium-samples
•
written
by
Chris •
You need the same GTF file that was used during mapping of your 10x data. If you used CellRanger and it was mouse (given you used recent CellRanger) you can find it in the CellRanger folder under refdata-gex-mm10-2020-A/genes/genes.gtf.
The input for the analysis is the BAM file from CellRanger, and the output is a loom file. loom can be read into R with LoomExperiment::import().
Thanks ATpoint! I think I can get the gtf file at: https://cf.10xgenomics.com/supp/cell-exp/refdata-gex-GRCh38-2020-A.tar.gz but velocyto also need repeat_msk.gtf which I don't know how to get after reading 10x website. Would you please have a suggestion?
velocyto run10x -m repeat_msk.gtf mypath/sample01 somepath/refdata-cellranger-mm10-1.2.0/genes/genes.gtf
Yes, thank you! I figured it out but I got an error:
velocyto run10x -m Hg38_rmsk.gtf WT/outs /cellranger/reference/refdata-cellranger-GRCh38-3.0.0/genes/genes.gtf
Inside WT/outs I have some folders such as filtered_feature_bc_matrix. Inside filtered_feature_bc_matrix, I have barcodes.tsv.gz
a possible suggestion: the masked repeats you were pointed to in this thread are from UCSC while your alignment was done using Ensembl reference?
Log in to answer this question.
