在R里面write.csv或者write.table时,quote默认是True,如果需要保存的文件中字符串两边不显示双引号的话,要设置quote=F

    1. #本流程主要为了得到pbmc3k的细胞barcode信息,为了下一步过滤pbmc3k的BAM文件
    2. library(SeuratData)
    3. pbmc3k=pbmc3k
    4. pbmc3k_meta=pbmc3k@meta[,c("seurat_annotations")]
    5. pbmc3k_meta$cell_id=rownames(pbmc3k_meta)
    6. pbmc3k_meta$`grouping_level_1`=c(rep("Control",1350),rep("Case",1350))
    7. pbmc3k_meta=pbmc3k_meta[,c("cell_id", `grouping_level_1`, "seurat_annotations")
    8. colnames(pbmc3k_meta)=c("cell_id","grouping_level_1","grouping_level_2")
    9. pbmc3k_meta$cell_id=paste0('CB:Z:',pbmc3k_meta$cell_id,'-1')
    10. write.table(pbmc3k_meta,'pbmc3k_meta.tsv',row.names = F, col.names = F,quote = F)
    11. 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?