教程:
https://www.jianshu.com/p/e65964f39eb3
对于
10X
的数据,可以直接从 比对后的.bam
文件开始,需要的文件及过程如下此步骤在
Shell
中以命令行的方式进行,后台调用Python
,要求Python>=3.6
1. 需要的文件
possorted_genome_bam.bam
:Cellranger
结果中的outs
文件夹中的文件,排序好的bam
文件
genome.gtf
: 在运行Cellranger
时使用的.gtf
文件;genes/genes.gtf
expressed repeats annotation
: 此步骤是选做
的,为了去除重复序列 (repetitive elements
) 的影响;此文件可以通过如下方式得到:
image.png image.png3.2 如果
UCSC genome brower
没有你的reference
版本,则需要自己手动注释:
https://yanzhongsino.github.io/2021/08/02/omics_genome.annotation_repeat/
Repbase: https://www.girinst.org/server/RepBase/index.php
2. 运行 10X
数据
velocyto run10x -m repeat_msk.gtf mypath/sample01 somepath/refdata-cellranger-mm10-1.2.0/genes/genes.gtf
-m
: 此项为expressed repeats annotation
; 此项为选做