课前准备--基因组(WGS/WES)联合单细胞获取突变信息
作者,Evil Genius
还是希望大家在单细胞的基础上往前多走一步。
今天我们来实现目标
严格实现这个目标,需要2个组学的数据,WGS/WES和单细胞数据(示例是单细胞核转录组)。
当然了,都需要bam文件的内容。
之前介绍过cellsnplite的方法,今天我们来更新一个方法。
第一步,检测WES/WGS的突变信息,这个关于基因组的全流程代码,包括SNV + ANNOVAR + MSI + Fusion + CNV + TMB ,课上我会发给大家,大家根据自己的情况进行调整,各种数据库需要大家自己补充了,供大家参考。
主要获取突变信息文件
第二步,分析单细胞的突变信息,获取单细胞级别的突变数据。
代码也很简单
bin/scReadCounts -r 单细胞bam文件 -s 基因组的突变信息文件(vcf格式) -G STARsolo -b Barcode文件 -o 单细胞突变信息矩阵