library(DropletUtils)
bc <- read10xCounts('/restricted/projectnb/waxmanlab/kkarri/tabula_muris/drop-seq/10X_P4_2_New_UseThis/outs/filtered_gene_bc_matrices_h5.h5')
bc$barcodenew <- gsub('-1','',bc$Barcode)
list <- bc$barcodenew
mol.info <- read10xMolInfo('/restricted/projectnb/waxmanlab/kkarri/tabula_muris/drop-seq/10X_P4_2_New_UseThis/outs/molecule_info.h5')
mol.df <- as.data.frame(mol.info$data[mol.info$data$cell %in% list,])
gene_reads <- mol.df[,c('gene','reads')]
gene_reads_df <- aggregate(gene_reads['reads'], by=gene_reads['gene'], sum)
bc_reads <- mol.df[,c('cell','reads')]
bc_reads_df <- aggregate(bc_reads['reads'], by=bc_reads['cell'], sum)
write.table(as.data.frame(gene_reads_df), "/restricted/projectnb/waxmanlab/kkarri/tabula_muris/drop-seq/10X_P4_2_New_UseThis/outs/gene_reads_df", )
write.table(as.data.frame(bc_reads_df), "/restricted/projectnb/waxmanlab/kkarri/tabula_muris/drop-seq/10X_P4_2_New_UseThis/outs/bc_reads_df", )
write.table(mol.info$genes,"/restricted/projectnb/waxmanlab/kkarri/tabula_muris/drop-seq/10X_P4_2_New_UseThis/outs/mol_info_read_gene", sep='\t')


library(DropletUtils)
library(dplyr)
library(stats)
bc <- read.table('/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G171/10X_premrna_G171B_mm10_ensembl_onlyPCG_transcript_UseThis/outs/barcodes.tsv')
bc$barcodenew <- gsub('-1','',bc$V1)
list <- bc$barcodenew
mol.info <- read10xMolInfo('/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G171/10X_premrna_G171B_mm10_ensembl_onlyPCG_transcript_UseThis/outs/molecule_info.h5')
mol.df <- as.data.frame(mol.info$data[mol.info$data$cell %in% list,])
gene_reads <- mol.df[,c('gene','reads')]
gene_reads_df <- aggregate(gene_reads['reads'], by=gene_reads['gene'], sum)
bc_reads <- mol.df[,c('cell','reads')]
bc_reads_df <- aggregate(bc_reads['reads'], by=bc_reads['cell'], sum)
write.table(as.data.frame(gene_reads_df), "/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G171/Analysis/gene_reads_df", sep='\t')
write.table(as.data.frame(bc_reads_df), "/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G171/Analysis/bc_reads_df", sep='\t' )
write.table(mol.info$genes,"/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G171/Analysis/molinfo_geneid", sep='\t')