在R里面write.csv或者write.table时,quote默认是True,如果需要保存的文件中字符串两边不显示双引号的话,要设置quote=F
#本流程主要为了得到pbmc3k的细胞barcode信息,为了下一步过滤pbmc3k的BAM文件
library(SeuratData)
pbmc3k=pbmc3k
pbmc3k_meta=pbmc3k@meta[,c("seurat_annotations")]
pbmc3k_meta$cell_id=rownames(pbmc3k_meta)
pbmc3k_meta$`grouping_level_1`=c(rep("Control",1350),rep("Case",1350))
pbmc3k_meta=pbmc3k_meta[,c("cell_id", `grouping_level_1`, "seurat_annotations")
colnames(pbmc3k_meta)=c("cell_id","grouping_level_1","grouping_level_2")
pbmc3k_meta$cell_id=paste0('CB:Z:',pbmc3k_meta$cell_id,'-1')
write.table(pbmc3k_meta,'pbmc3k_meta.tsv',row.names = F, col.names = F,quote = F)
write.table(filter,'filter.txt',row.names = F, col.names = F,quote = F)
后面根据filter文件过滤BAM文件的教程见:Is there way to filter the BAM file produced by 10x pipelines with a list of barcodes?